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Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


TURBULENCE IN CIVIL ENGINEERING: 
AN INTRODUCTION TO THREE RESEARCH PAPERS 


J. M. Robertson,* M. ASCE 
for Hydraulics Research Committee 
(Proc. Paper 1391) 


The “Symposium on Turbulence” for which this paper serves as an intro- 

| duction is more than a mere collection of papers about turbulence. It is evi- 

dence of the hydraulic engineer’s interest in learning more about that funda- 
mental subject. Although as engineers we are interested in practical 
problems, we also desire to understand the basic mechanics of the phenomena 
that enter into our practice. Thus hydraulic engineers are interested in 
more than just coefficients. In many instances a clearer understanding of 
the flow fundamentals has lead to advances in our methods of analysis and 
design. For example, increased knowledge of turbulence has led to advances 


in our practical approach to pipe flow problems and sediment transport 
problems. 


Turbulence is a subject of concern to civil engineers in many areas be- 
side just those interested in flow in various channels. Thus, the wind loads 
on buildings and other structures are variable due to large-scale turbulent 
fluctuations. The 1952 paper by R. H. Sherlock(!) is an excellent illustration 
of the structural engineer’s concern with turbulence. The group of three re- 
search papers which follow are limited to cases of primary concern to hy- 
draulic engineers; however, they do contain information on turbulence funda- 
mental to all areas. Reflecting an interest in turbulence beyond that speci- 
fied by strictly hydraulic problems, the Committee on Research of the So- 
ciety’s Hydraulics Division has encouraged and sponsored studies and re- 
views of the nature of turbulent motion. The papers which follow are pre- 
sented as further evidence of our conviction that the civil engineer, and 
particularly the hydraulics engineer, desires to be cognizant of turbulence 
and appreciates this knowledge. 

The function of this introduction is to indicate where and how the three re- 
search papers conform to the general status of our knowledge of turbulence. 
| Since most of the readers are presumed to be ASCE members, this will be 
|. attempted through a review of previous coverages of turbulence in ASCE 


Note: Discussion open until March 1, 1958. Paper 1391 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No. HY 5, October, 1957. 


*Prof. of Theoretical and Applied Mechanics, University of “iinois, Urbana, 
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publications. The literature on turbulence is voluminous and far exceeds that 
in ASCE publications; an attempt is made therefore to assess the adequacy of 
such a coverage. 


ASCE Publications on Turbulence 


Presentation of the “modern” theory of turbulence to ASCE members 
dates back to the 1936 paper by Hunter Rouse(2) and the 1939 paper by A. A. | 
Kalinske. Rouse’s presentation was principally concerned with the effects 
on flow in pipes—head loss and velocity distribution relations. Analytically, 
the results presented were based on the phenomenological theories of 
Prandtl’s mixing length and von Karman’s similarity hypothesis both very 
valuable concepts today for practical problems. In concluding his paper, 
Rouse noted the need for a complete rational analysis of the inner mechanism 


of turbulence movement. Kalinske presented many of the basic concepts of i 
“the statistical theory of turbulence” initiated by G. I. Taylor—the correlation 

coefficient, length scales, spectrum of turbulence, etc. The knowledge of | 
turbulence was correlated with the solutions of certain practical hydraulic J 


problems; energy dissipation, suspended material transport, etc. Kalinske’s 
experimental data were primarily obtained by laborious photographic tech- 

niques; he indicated a great need for an instrument, like the hot-wire anemo- 
meter used in air, that could be used for turbulence measurements in water. 

A review of the ASCE Transactions and Proceedings since the papers by 
Rouse and Kalinske indicates papers on turbulence and suspended material 
by Dobbins (4) (1943), Vanoni (5) (1944), Camp (6) (1945), Ismail (7) (1951); 
papers on energy conversion by Kalinske'8) (1944) (based on work sponsored 
by the Hydraulic Research Committee) and Bakhmeteff and Allen'9) (1945) 

(also based on an investigation under the auspices of this committee). Con- 

cern with atmospheric turbulence is represented by the paper by Sherlock(1). 

More recently we have the 1953 paper by sence on the turbulent boundary 

layer on steep slopes and the 1955 paper by Ross(11) on the turbulent boundary aa 
layer in adverse pressure gradients. We have also papers by Bennett and 

Lee (12) (1955) on boundary layer transition and by Einstein and Li (13) (1956) 

on the laminar film or sublayer. 

These are interesting papers—all of them from a turbulence consideration = 
and many of them from the standpoint of practical hydraulics problems. It is | 
significant to note that they are primarily concerned with one or the other of 
the practical problems introduced by Rouse or Kalinske. 

Referring to the papers which follow it is quite appropriate that two of 
them involve solutions to the need expressed by Kalinske for a turbulence 
measuring instrument for use in water. The hot-wire anemometer (so satis- 
factory for use in air) is not suitable in water and two of its hydraulic cousins 
are described. A third cousin, the hot-film anemometer, has recently ap- 
peared 14) put is not discussed. The experimental results to be presented in 
the three papers also help us in our progress towards fulfilling the need ex- 
pressed by Rouse twenty years ago for an understanding of the turbulence 
mechanism. However, progress in this direction is quite slow and there is 
much yet to be learned. 


| 
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Status of Turbulence Knowledge 


To assess the extent to which the ASCE, through its publications, is keep- 
ing up on turbulence it is necessary to take a quick look at the overall 
progress in the field. 

Activity in the field of turbulence during the present century has been con- 
centrated in three periods: in 1925-30 when the phenomenological theories of 
Prandtl, von Karman, etc. were developed; again in 1935-37 when the statisti- 
cal theory of turbulence saw rapid development and the first real measure- 
ments of turbulence were obtained; and again in the 1947-53 post war period. 
In this last period a recent review\15) of the significant papers on the 
measurement of turbulence indicates the publication of twenty-nine such, with 
eight in the year 1951. Although theoretical development continued, this was 
principally a period of learning from experiment more and more about the 
nature of turbulence. 

It would obviously require a paper several times this in length to present 
an adequate review of the present status of our knowledge on turbulence. 
However, a few remarks will be attempted in the hopes of giving a glimpse of 
what we know or think we know. It has been noted that Kalinske introduced 
certain of the basic elements of “the statistical theory of turbulence.” In 
spite of this most of our applications have been in terms of the results of the 
phenomenological theories involving the mixing length, momentum transfer, 
etc. As of today, most of the leading turbulence researchers have discarded 
these concepts. Although it is true that actual measurements of turbulence 
discredit the mixing length approach, we are not in a position to replace it 
with a better theory for practical solutions. Therefore, we continue to use it, 
but with caution. The statistical theory of turbulence, which for a long time 
has been limited to the rather idealized case of homogeneous isotropic turbu- 
lence, 16) is just starting to be used in an attack on the more complicated 
case of turbulence in shear flows. 

The three research papers which follow faithfully reflect this introduction 
of the statistical theory in to shear flow problems. The study of correlation 
coefficients and frequency spectrum measurements in the first paper is a 
direct result of attempts to tie in with the statistical theory. The research 
reported did not progress as far with the statistical theory as some have, (17 
due primarily to limitations in experimental equipment. In two of the three 
papers this was a direct result of the fluid medium used, namely water. The 
fact that progress in water has been as good as it has is very encouraging. 

From the temporal- mean velocity standpoint our understanding has 
changed little since the summary paper by Rouse. We still conceive of a thin 
laminar film separating the wall from the turbulent flow. It is true there are 
various indications that our idealized picture of this is not correct and that 
even the laminar film contains turbulence fluctuations which, however, de- 
crease in magnitude as the wall is approached. But the laminar film concept 
is still valuable and useful. Outside this thin film the turbulent velocity pro- 
file is found to follow the logarithmic type relation due to Prandtl and von 
Karman. However, the 25 year old concept of a universal profile, deduced 
from Nikuradse’s experiments, is now known to be a poor approximation. A 
single logarithmic law does not adequately describe the velocity variation 
from the edge of the laminar film to the center of the pipe or outer edge of 
the boundary layer. Instead the logarithmic velocity distribution applies for 


| 
| 


1391-4 HY 5 October, 1957 


a distance of about ten or fifteen percent of the pipe radius from the wall. 
Outside of this region, in the central or core region of pipe flow, other re- 
lations apply. Description of the velocity distribution in these several regions 
is still based on phenomenological concepts. 

Progress towards our understanding of turbulent shear flows (i.e. pipes 
and boundary layers) is being made with the aid of energy balance consider- 
ations as well as the statistical theory. In this additional approach, one 
considers the production, convection, diffusion and dissipation of the turbulent 
energy. Basic discussions along these lines were introduced by von Karman 
(1937) and Prandtl (1945). The approach is very similar to the analysis util- 
ized by Bakhmeteff and Allen. Some very interesting results of this ap- 
proach are now apparent. It is found that, although turbulence production and 
dissipation are very high near the wall, there is significant transport and dif- 
fusion of turbulent energy to other regions. 


Future Plans of Hydraulics Research Committee 


This present symposium is not a summary of our knowledge of turbulence, 
with civil engineering applications and interests, but rather a report on some 
current researches into turbulence. Your Committee on Research of the 
ASCE Hydraulics Division feels that it should do more in this regard than just 
collect groups of papers when they become available. With this idea in mind, 
discussions have been held with several members of our profession who are 
active in the turbulence field. 

This has led to an indication of specific phases of turbulence that should 
and can be summarized in papers at this time. Some of these are as follows: 
application of basic ideas of turbulence to liquid mixing problems; application 
to micrometeorology or surface runoff, to flow through porous media, open 
channels, hydraulic structures and tidal estuaries; and effects of turbulence 
on sediment transport, river pollution problems and heat transfer in flumes 
and rivers. Along this line we have the suggestion that the sediment sus- 
pension problem should be amenable to analysis, following energy balance 
considerations. If so, such an analysis could give us a much better under- 
standing of suspension. 

It is impossible to predict the extent to which the above subjects will be 


reported in the near future. However, it is hoped that at least some of them 
will be. 
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HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


TURBULENCE IN CIVIL ENGINEERING: 
MEASUREMENTS IN FREE SURFACE STREAMS* 


Arthur T. Ippen** M. ASCE and Fredric Raichlen*** J.M. ASCE 
(Proc. Paper 1392) 


Various turbulence phenomena in free surface streams were investivated 4 


by means of a combination of a total head tube with a capacitance type pres- 
sure transducer. Design criteria for the instrument are given as well as a 
series of performances data derived from turbulence measurements. 


SYNOPSIS 


Certain turbulence characteristics of free surface streams were deter- 
mined in the laboratory by means of a combination of a total head tube with a 
capacitance type pressure transducer. The design criteria for such an 
instrument are given together with performance data. 

The experimental studies covered primarily boundary layer growth in a 
high velocity stream, the relative turbulence intensity, the properties of the 
auto-correlation coefficiént and of the mean intensity spectrum in relation to 
depth. The purpose of the studies was to check the instrument performance 
in relation to turbulence properties known from hot wire measurements in air 
streams. 

In general, satisfactory agreement was found and it may be concluded that 
an instrument of this type gives reliable values of mean velocities, root mean 
square values of the turbulent velocity fluctuations in the direction of the 
mean flow and the scales of macro turbulence in water. 


Note: Discussion open until March 1, 1958. Paper 1392 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 5, October, 1957. 

*Presented at a meeting of the A.S.C.E. in Pittsburgh, Pa., October 18, 
1956. 
**Prof. of Hydraulics, Hydrodynamics Lab., Massachusetts Inst. of Tech- 
nology, Cambridge, Mass. 

***2nd Lt., MeClellan Air Force Base, California, formerly Research Asst., 
Hydrodynamics Lab., Massachusetts Inst. of Technology, Cambridge, 

Mass. 


1392-1 


| 
| 
i 
ee HY 5 
} 
| 
| 


October, 1957 
INTRODUCTION 


Turbulent motion in fluids is of interest to hydraulic engineering in con- 
nection with many practical problems. These include in general the me- 
chanism of energy dissipation and the correlated mixing processes of waters 
with different temperatures, salinities, oxygen content and solid suspensions. 
More specifically the transport of sediment in streams and by wave action 
along beaches and the diffusion of waste waters in streams and reservoirs 
may be mentioned. 

Knowledge concerning the statistical properties of turbulent streams has 
been rapidly advanced in recent years through the development and extensive 
experimental application of the hot-wire anemometer for turbulent air flows. 
However, a suitable instrument for water has been lacking, although many 
methods have been tried. The hot wire is affected by the impurities in the 
water, which do not permit a constant calibration to be maintained.(1) The 
study of particle motion with oil droplets or other suspensions is very 
cumbersome and unsuited to the exploration of local conditions. 

The following study reports summarily on a program initiated in 1950 at 
the Hydrodynamics Laboratory of M.I.T. to combine the familiar Pitot tube, 
normally used for mean velocity determination, with a rapidly responding 
pressure pickup, so as to register the time history of the instantaneous 
pressures exerted by the turbulent flow on the tip of the tube. 

Of all the systems considered initially the capacitance type pressure cell 
proved to be the most suitable pickup with respect to possible sensitivity and 
frequency response. In the gage developed the dynamic head existing at a 
given point in the flow is converted to a pressure, which is transmitted 
through the tube to the pressure cell, where it causes a minute deflection of 
a clamped circular diaphragm. This deflection reduces the small airgap 
between the diaphragm and a fixed electrode and results in a change in 
capacitance, which is electronically amplified and recorded. 

For each point in the flow field, these records must permit an adequate 
definition of the state of turbulence by means of the following statistical 
parameters: 


1. The root mean Square of the velocity fluctuations in the direction of 
mean motion: u' 

2. The auto-correlation factor R(-+) , which compares these random 
fluctuations at one point with time. It is defined as: 


li T 
Treo u'lt)u'(t+7) dt 
Rr) = 
u 
wherein T indicates a delay time in seconds. 


. The mean intensity spectrum, which gives the distribution of mean 
power in the turbulence as a function of the frequency. The portion of 
kinetic energy for a frequency interval df is defined by u'2 W(f)df 
wherein according to the Wiener-Khinchine theorem the spectrum func- 
tion is: 


W (Ff) = af Rex) cos WT dT (2) 


| 
| 
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In general, the turbulent kinetic energy of the flow is given by the mean 
square of the fluctuations, while the distribution of the auto-correlation 
factors with respect to delay time results in definitions of certain turbulent 
eddy characteristics. The smallest eddies in the flow are defined by the 
shape of the correlation curves near their origin and are described by the 
scale of micro-turbulence A . By contrast, the entire area under the auto- 
correlation curve defines the scale of macro-turbulence or the mean eddy 
size. The total area under the mean intensity spectrum employed here is 
evidently only the portion of the total kinetic energy of turbulence defined by 
u'2. The spectrum function as defined by Equation (2) represents a com- 
parison of the auto-correlation function over its entirety with cosine functions 
of various frequencies w. 


Mechanics of Total Head Tube-Pressure Cell Combination 


Before discussing the gage actually used in these turbulence investigations, 
the basic mechanics of such an instrument along with the criteria to be em- 
ployed in design will be described. 

A schematic drawing of the gage is shown in Figure 1. The gage consists 
essentially of two parts; the total head tube which is inserted in the flow, and 
the pressure transducer mounted above the flow to which this tube is directly 
connected. At one end of the pressure chamber is a flexible diaphragm firm- 
ly clamped to the chamber walls around its periphery. A fixed electrode is 
placed a small distance away from the membrane separated by a thin air gap. 
A change in the instantaneous total head at the tube tip causes a change in 
pressure at the diaphragm face with a resulting deflection of the membrane. 
This deflection is recorded by electronically measuring the change in 
capacitance between the fixed electrode and the diaphragm. By proper static 
calibration this change in capacitance can be directly related to an instan- 
taneous change in total head. 

Considering the total head tube inserted in an open channel and aligned 
with the direction of mean tiow the instantaneous total head recorded by this 
gage is: 


wherein: U temporal mean velocity in the x direction 


u' = turbulent velocity fluctuation in the x direction 
p' = turbulent local pressure fluctuation 
hg = mean hydrostatic head on the tube tip 


Thus also: 


U 2 Vu ’ u' 2 p ‘ 
Ss 2g 9 29 
The values of H and H; are varying in random fashion. Assuming that 
the gage records the instantaneous values of Ht correctly it is seen that the 


difference between the instantaneous velocity u = y 2gH and the temporal 
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mean velocity U is primarily dependent upon the product = . While the 


2 
' ' 
— and ‘ terms are of the same order of magnitude, they are normally 
not in phase. Thus it is conceivable that for very small values of u' approach- 
' 
ing zero the turbulent pressure e exceeds in value the other —— 
terms. However, assuming for example u' max = .10 U and thus Me = 


' 
wl it is apparent that the maximum contribution of a to H 


when u' = 0 is only of the order of one percent and its influence on u' not 
more than 5%. The errors inherent in the instrument are given in more de- 
tail in reference (4). 

The method employed in the analysis of time histories of instantaneous 
head will be taken up later. Of importance is the fact that the dominant ran- 
dom quantities in the above equation must be recorded by an instrument, 
which will not distort these fluctuations felt at the tip of the tube. It is clear 
from dynamic considerations that the design criterion of prime importance 
for accurate measurement of these fluctuations is the natural frequency of 
the gage. 

In order to determine the natural frequency theoretically, the mechanical 
system of the gage is represented by a one-dimensional system composed of 
a spring (the diaphragm), a mass (the combined effective mass of the water 
in the impact tube and the pressure chamber), and a dashpot for viscous 
damping. By using this approach to the problem, many graphs and para- 
meters familiar to engineers can be used. This one-dimensional analysis 
has been shown satisfactory from both experimental and theoretical con- 
siderations.(4) 

For the system described above, the natural frequency can be obtained 
either by the use of Lagrange’s equation or by adding up the forces acting on 
the system; i.e., the sum of the inertial, spring and viscous forces being set 
equal to the sum of the externally applied forces. In both methods the follow- 
ing major assumptions are made in the evaluation of the resonant frequency. 


1. Velocity distributions in the tube and in the diaphragm well are uniform 
and may be related by the continuity equation: 
A2 
m4 = 
1 
where, m, = the velocity of water in the tube 
mg =the velocity of water in the diaphragm well 
A; =the area of the tube bore 
Ag =the area of the diaphragm well 
. The kinetic energies of the water in the diaphragm well and of the 
diaphragm are considered small compared to the kinetic energy of the 


A2 
water in the tube, since (2 ) >> 4. 


. The effect of viscous damping on the resonant frequency of the instru- 
ment is assumed to be small. 


| 
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The expression defining the natural frequency of the gage, in units of cps, 
derived by this type of analysis is thus: 


/ k 


eT \A; pA L, 


where, k = spring constant of the diaphragm 
Ly = tube length 
p =density of water 


The spring constant of the diaphragm treated as a circular membrane 
clamped at its edge is: 


a? (I- €?) 


where, a = radius of the diaphragm 
h = thickness of the diaphragm 
E = modulus of elasticity of the diaphragm 
£ = Poisson’s ratio = .30 


The instrument first developed at the Hydrodynamics Laboratory and used 
for the initial turbulence studies described herein is designated as Gage 1 and 
is shown in Figure 2. The “pressure pick-up” which was inserted in the flow 
was a Small tube, with a vertical leg 6 in. long and a 2 in. horizontal leg. The 
tube was of stainless steel with a 1/8 in. outside diameter and a 1/16 in. 
inside bore diameter. It was connected directly to the pressure transducer 
which was located above the water surface. The diaphragm was 0.0075 inches 
thick and 1 in. in diameter and was made of beryllium copper with a modulus 
of elasticity of 15 x 106 psi. The brass electrode, 0.938 in. in diameter, was 
insulated from the gage housing by a sheath of hard rubber and was separated 
from the diaphragm by an air gap of no less than 0.0025 inches. A small 
cylinder valve attached to the pressure cell allowed the system to be filled 

by flushing water through the diaphragm well and out through the impact tube. 

With these physical characteristics of Gage 1, a natural frequency of 15.4 
cps is obtained from Equations 4a and 4b. 

At first glance, this frequency response appears to be too low for an in- 
strument to be used in turbulence measurements. However, the results ob- 
tained with this instrument are in good agreement with those obtained using 
a modified version of Gage 1 (designated as Gage 2) having a ten times higher 
natural frequency. This would indicate that the major portion of the turbulent 
energy in the type of flow investigated is of low frequency. This conclusion 
will be supported by experimental data presented in connection with the re- 
sults of turbulence investigations with Gage 2. 


Turbulence Measurements with Gage 1 


The major portion of the measurements with Gage 1 and its successors 
was carried out in an open channel in the laboratory. This flume has 


| 
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aluminum walls, a glass bottom and is 45 ft. long and 4 ft. wide. The slope of 
the flume can be varied fully from zero to 0.1 ft. per ft. The mean velocities 
ranged from 4 to 10 fps and the depth of flow was varied from 1 to 3 inches. 
The circulation of water is accomplished by pumping from a large reservoir 
in the basement of the laboratory through a 10 x 6 in. Venturi for metering 
purposes and a 10 in. gate valve into a closed transition section. This sec- 
tion transforms the flow cross-section from a circular one of 10 in. diameter 
to a rectangle 10 in. x 4 ft. After this the flow passes upward to the elevation 
of the flume bottom and into the nozzle section, where several honeycomb 
baffles eliminate large scale disturbances. The flow is accelerated rapidly 
through a 4 ft. wide adjustable nozzle and enters the so-called high velocity 
flume with a free surface with practically negligible turbulence. 

E. M. Lowry, Jr.(2) in 1951 used Gage 1 to investigate the variation of the 
temporal mean velocity and of the root mean square of the turbulent velocity 
fluctuations depthwise in a uniform free surface flow. The results obtained 
(see Figure 3) were compared to measurements of Laufer(3) made in an air 
stream in a rectangular test section at a comparable Reynolds number and 
indicated good agreement. 

In obtaining the mean velocities and the root mean square values of the 
turbulent velocity fluctuations, it was assumed that each horizontal line on the 
Brush recorder paper that was used for signal display corresponded to a 
particular velocity head found by static calibration of the gage. Knowing the 
head of water H necessary to cause a given deflection on the recorder, the 


velocity was obtained from the relation u = V 2gH . Thus, in analyzing the 
record, it was assumed that every time the trace touched or crossed a parti- 
cular horizontal line on the recorder paper, a certain velocity associated with 
that line existed instantaneously. Designating the frequency of occurrance of 
the velocity as f, and the total velocity at that point as u, the expression 


(5) 
=F 


was used to determine the value of the temporal mean velocity U for a record 
of suitable length. For each record so analyzed, the value of the root mean 
square of the velocity fluctuations was computed from: 


Yur? = (6) 


These expressions actually evolve from Equation 3 after neglecting the 
pressure fluctuations, p'. 

Immediately downstream from the nozzle of the flume a turbulent boundary 
layer begins to develop which grows from the bottom and eventually extends 
over the full depth of the flow within a distance dependent upon the flow condi- 
tions established. For more detailed confirmation of the applicability of the 
gage Lowry investigated this boundary layer growth at stations every six 
inches downstream from the exit of the nozzle for a distance of four feet. The 
thickness of the turbulent boundary layer was determined by varying the tube 
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submergence slowly from the smooth glass bottom of the flume upwards in 
the flow until the velocity fluctuations existing in the boundary layer and ob- 
served in the record decreased to a few short intermittent bursts and even- 

tually disappeared. The values of the boundary layer thickness as shown for 
two runs in Figure 4 appear in good agreement with the well known expres- 

sion: 


3 0.38 


= 
x 


where, 6 = boundary layer thickness 
x = distance downstream 


(7) 


Rx Umax 
v 

Water surface elevations and theoretical boundary layer displacement 
thicknesses are also given. 

In addition the statistical distribution of the velocity fluctuations about the 
temporal mean velocity was found nearly Gaussian in form which is to be ex- 
pected for a random variable. This is shown in Figure 5 for numerous points 
at the same station. 

In 1952, D. L. Favin constructed an electrodynamometer at the laboratory 
and J. P. Lawrence(5) incorporated it into the circuit of Gage 1 to obtain, 
directly, values of the temporal mean velocity and the mean square of the 
fluctuations. This dynamometer consisted essentially of two sets of interact- 
ing coils, one set fixed and the other set suspended by a spring. When a cur- 
rent is imposed in both sets, an angular deflection of the novable set is ex- 
perienced which is proportional to the mean of the product of the applied 
currents. 

The actual values of interest can be obtained by applying the following 
currents to the coils: 


1. By imposing a constant current on one set of coils and a fluctuating 
signal from the gage on the other set, the value of the temporal mean velocity 
was determined. 


2. By removing the D.C. level of the signal and sending only the A.C. 
portion through both the fixed and the movable coils, a squaring and meaning 
effect is produced from which the mean square of the fluctuating velocity com- 
ponent is determined. 


A typical example of Lawrence’s results of the variation of the relative 
mean velocity and relative turbulent intensity with relative depth for Reynolds 
number value of R = 210,000 is shown in Figure 6. The different check points 
that were made both with a Prandtl tube and at a later time with Gage 1 are 
also given in this figure illustrating the reproducibility of the data obtained 
with this instrument. Figure 7 gives a summary of results for relative tur- 
bulence intensities obtained for a wide range of Reynolds numbers. 

The dynamometer is also capable of giving correlation factors directly; 
however, this possibility was not pursued, since a magnetic recording and 
playback of the record would have been required for auto-correlation or a 
duplicate instrument for space correlation. Instead the attention was directed 
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towards the development of a gage with much higher natural frequency. 


Turbulence Measurements with Gage 2 


Mechanical Features 


This new instrument denoted as Gage 2, was designed for a natural fre- 
quency of better than ten times that of the original capacitive element. Re- 
ferring to Equations 4a and 4b, it is seen that this can be accomplished in 
several ways by: 


1. increasing the spring constant k of the diaphragm 
2. increasing the area ratio, Aj/A9 


3. decreasing the impact tube length 


Of most practical significance is the increase of the spring constant of the 
diaphragm. Hence, in the first version of Gage 2 the diameter of the dia- 
phragm in Gage 1 was retained while the thickness was increased from 
0.0075 in. to 0.030 in. This change alone will result in a natural frequency 
of approximately eight times that of Gage 1; i.e., approximately 112 cps. 

The possibility of viscous damping to attain a relatively flat response curve 
was investigated by dividing the internal portion of the gage into essentially 
two parts; (1) a damping chamber that could be oil-filled for viscous damping 
and (2) an impact tube connection that would be water-filled and 

separated from the damping chamber by a thin secondary diaphragm. It was 
found that mechanical damping even with the damping well filled with an ex- 
tremely viscous oil was not appreciable. However, the secondary diaphragm 
was nevertheless retained, with the gage completely water-filled, in order to 
achieve a higher frequency response. Since it was found that a small amount 
of air trapped in the gage would substantially reduce the frequency response 
of the instrument, the gage was filled while submerged in a container of 
deaerated water. 

The features of this modified gage may be seen in Figure 8-a, a schematic 
drawing of Gage 2, along with Figures 8-b and c showing photographs of the 
instrument and of its associated electronics. 

A typical set of records obtained using Gage 1 and Gage 2 is shown in 
Figure 9. In these cases the gages were positioned in the high velocity flume 
approximately 2 ft. downstream from the nozzle exit. This station is up- 
stream from the point where the turbulent boundary layer reaches the water 
surface. The first record for each gage was obtained outside of the turbulent 
boundary layer where the small fluctuations evident in the signal may be con- 
sidered to be low intensity free stream turbulence. The second set of 
records for each gage was obtained near the channel bottom: i.e., within the 
boundary layer and thus the shear zone. 

It was found at a later date that the secondary diaphragm mentioned above 
introduced certain instabilities in Gage 2. For this reason, it was eliminated 
and the resultant decrease in the natural frequency was compensated for by 
increasing the thickness of the primary diaphragm still further by a factor of 
two (from 0.030 in. to 0.060 in.). This increased the calculated natural fre- 
quency of Gage 2 to approximately 255 cps. 

Due to the frequency limitations of the Brush Recorder used in conjunction 
with Gage 1, a Dumont type 322 dual beam oscilloscope was employed with 
Gage 2 with the signal being fed into one beam of the oscilloscope while 
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simultaneously a square wave of 60 cps was displayed on the other. The lat- 
ter afforded both a time and a reference trace for the signal as illustrated in 
the lower portion of Figure 9. Recordings were made by photographing the 
oscilloscope by means of an open-shutter, 35 mm. motion picture camera. 


Data Analysis and Auto-Correlations 


A micro-film viewer was used in the analysis of the data and a transparent 
grid was placed over this enlarged record with time intervals assigned to the 
vertical grid lines by means of the 60 cps square wave. Using this grid ar- 
rangement, values of the amplitude of the signal were obtained from the 
record at equally spaced time intervals. Approximately 1200 such points 
were then submitted to the Digital Computer Laboratory of M.I.T. where a 
program had been evolved for obtaining values of the auto-correlation factor 
from T=0 to Tt = 300 7) , where T, is the time spacing of the raw data 
points. 

Employing the method described above, auto-correlation curves were ob- 
tained at one depth in each of two flows having Reynolds numbers of 52,000 
and 54,300. These results are presented in Figures 10 and 11. 

It is apparent that these curves have both a high and a low frequency 
periodicity superimposed on the random portion. In examining the form of 
the numerator of the auto-correlation function (Equation 1), it is seen that 
for a periodic signal, a periodicity occurs in the correlation function with the 
same frequency as that in the signal but with no phase relation. The high fre- 
quency component is due to the natural frequency of the gage and affords a 
method of experimentally determining the frequency response of the instru- 
ment (240 cps in this case compared to 255 cps calculated theoretically). 

By the use of the correlation curve, two characteristic lengths have been 
defined by Taylor(6) which help to describe the properties of the turbulent 
flow field in more detail. The scale of micro-turbulence, associated with the 
smallest eddies in the flow is defined as the intercept, on the time or distance 
axis, of a parabola fitted to the correlation curve in the vicinity of the origin. 
In order for this parameter to be used in this investigation confidence must 
be placed in the portion of the correlation curve in the vicinity of T = 0. With 
this type of gage, where the mechanical system can be said to be composed of 
essentially a spring, a mass and a dashpot, the signal becomes distorted near 
the resonance and at forcing frequencies greater than the natural frequency 
(240 cps or t= 0.00415 sec). Therefore, reliable information cannot be ob- 
tained from the correlation curves in this region, and hence, the scale of 
micro-turbulence cannot be determined. 

The second length termed the scale of macro-turbulence, L, is associated 
with the entire turbulent field and is a measure of the average eddy size 
present. It is defined in units of time as: 


= dtr (8) 


A macro-scale larger in one case than another indicates that the auto- 
correlation curve falls off more gradually, and there is better correlation at 
greater delay times in the former as compared to the latter. 

As has been shown by Dryden(7) and Laufer,(3) a curve of the form 
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=e (9) 


can often be fitted to the correlation curve away from the origin for cases of 
a high turbulence Reynolds number. With this function fitting the data, then: 


Ly 
Lr=fe dt = (10) 


The macro-scale was determined from these curves first by eliminating 
the 240 cps component from the curves received from the Digital Computer 
Laboratory by a three point averaging process. In this manner a correlation 
curve consisting of 300 weighted mean values of the correlation factor was 
constructed. The curves so resulting were then averaged by eye, the final 
smoothed curves being shown in the above mentioned figures. This facilitated 
the subtraction of any low frequency periodicities from the correlation curve 
before proceeding with the determination of the macro-scale. A cosine func- 
tion having the same amplitude and period as the predominant low frequency 
component was subtracted from the smoothed correlation curves. The ordi- 
nates of the curve which resulted after the cosine function was removed com- 
prised the curve to which Equation (9) was fitted to determine L,. The first 
third of these points was then plotted on semi-logarithmic paper, ln R(T) vs rT. 
Since the function e- TU/Lx treated in this manner reduces to a straight line, 
this was a relatively simple method of determining the mean eddy size. The 
analytical expression fitted in this manner shows rather good agreement with 
the meaned experimental data away from the origin. From this curve fitting 
process the values of L,/Yo were obtained for the relative depths of .157 and 
-5910 as 2.95 and 3.66 respectively. 

The magnitude of these results is reasonable, since other investigators 
have found that the scale Ly in a depth-wise direction is approximately 1/3 to 
1/2 of that in the downstream direction. If this factor is applied to these re- 
sults with y, = 1.27 and 1.09 in., the depth-wise macro-scale would be of the 
same order as the depth of the flow. In addition to this, Townsend(8) has in- 
dicated that “eddies of downstream length . . . over three times the boundary 
layer thickness, contain much of the turbulent energy... .” 


Mean Intensity Spectra | 
The mean intensity spectrum, which is the cosine transform of the auto- | 

correlation function determines the distribution of turbulent energy with fre- if 

quency. The Digital Computer Laboratory which assisted in the evaluation of 

the mean intensity spectra used a finite length of data in determining the 

cosine transform of the correlation function. Thus instead of an upper limit 

of infinity in the integral in Equation (2), the limit is actually the effective ui 

time length of the correlation curve, T;- In addition, the cosine transform bal 

is convolved with the cosine transform of a unit step function of length 7, . | 


Hence, the mean intensity spectrum obtained from the Digital Computer 
Laboratory and denoted as W(f)" is actually: 


| 
| 


ASCE IPPEN - RAICHLEN 1392-11 


W(t) = 


If a periodicity exists in the signal, this periodicity will be apparent as a 
spike at that frequency in the mean intensity spectrum caused by the quantity 
W(f). On both sides of the spike will be a series of peaks and troughs of 


SIN WT, 


amplitude T;/ a , wherein n = 3, 5, 7, 9, . . . occurring with a period of 
0.8 7, due to the term = ‘ 


A typical example of the mean intensity spectra obtained from the Digital 
Computer Laboratory is shown in Figure 12 for three frequency ranges; 0-25 
cps, 0-75 cps and 0-250 cps. 

Apparent in this curve is a negative spike at zero frequency. This point on 
the mean intensity spectrum corresponds to the ordinate of the auto-correla- 
tion curve at T = oo; i.e., the mean or D.C. level of the auto-correlation 
curve. The primary data frém which these curves were obtained were taken 
from the record about an assumed mean velocity level. Therefore, if this 
assumed mean level is correct, the auto-correlation curves should approach 
a line of zero correlation in the results received from the Digital Computer 
Laboratory. However, if this value is incorrect and the Computer fails to 
remove all the D.C. level from the signal in computing the auto-correlation 
factor, the line of zero correlation obtained with its assistance will possess 
either a positive or a negative value. A negative value of this line of zero 
correlation, such as was actually experienced, resulted in the negative spike 
at the origin of the mean intensity spectrum. 

Periodicities below 5 cps are evident in this spectrum and are believed to 
be caused by periodic components in the flow. The apparent irregularity of 


this spectrum can be attributed to the term ———__— in Equation (11) which 


super-imposes curves of this form throughout the spectral range, its effect 
being most noticeable on either side of the predominant periodicities. 

A procedure was developed to smooth these curves by fitting the cyrve 
sin 


— to the peaks and troughs on either side of the predominant spikes 


caused by the periodic components in the flow. In this manner a mean inten- 
sity was obtained at each predominant periodicity which would be the actual 
relative intensity with the absence of the spike and the elimination of the ef- 
fect of the sine term in Equation (11). Minor protuberances were remgved 
from these curves by eye. The smoothed curve obtained in this manner from 
Figure 12 is shown in Figure 13. 

Using the exponential expressian e~ TU/Lx for the correlation function 
and the definition for the mean intensity spectrum; 4.e., Equation {2), the fol- 
lowing analytic expression is obtained for the spectrum function: 


4 
I+ 


= 


(12) 
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This expression was fitted to the smoothed mean intensity spectrum ob- 
tained in the manner described above. The value of the turbulence macro- 
scale employed in this reference curve was that obtained from the auto- 
correlation curve. The agreement between experimental function and analytic 
reference points in Figure 13 indicates that the low frequency periodicities 
were successfully removed in the first L, evaluation; i.e., using the auto- 


correlation curve of Figure 11. 
The mean intensity spectrum shows a rapid decrease in energy with in- 
creasing frequency, with over 50% of the turbulent energy lying below 5 cps. 
The natural frequency of the instrument is evident in Figure 12 as a band of 
low amplitude spikes from approximately 200 to 235 cps. Thus both the auto- 
correlation curve and the mean intensity spectra can be employed to evaluate | 
the natural frequency of such an instrument. i 
A loss of sensitivity of the electronic circuit was experienced in subse- 
quent experiments which could not be corrected electronically without com- 
pletely redesigning this circuit. However, this was improved by decreasing | 
the thickness of the primary diaphragm from 0.060 to 0.045 in. Though this | 
lowered the natural frequency of the instrument to approximately 150 cps, it 
is evident from an examination of the mean intensity spectrum that this is high 
enough to record the turbulent fluctuations with little distortion. 


Turbulence in the Wake of a Cylinder 


The primary objective of this phase of the investigation was to measure 
variations in the temporal mean velocity, the root mean square of the turbu- 
lent velocity fluctuations and the scale of macro-turbulence in the wake of a 
right circular cylinder in super-critical open channel flow, thereby demon- i} 
strating the applicability of a total head tube-pressure cell combination in 
such investigations. A brass cylinder with a diameter d = 1 in. was used for 
testing in the high velocity flume in the laboratory. The cylinder was placed 
with its vertical axis 18 ft. from the nozzle exit, on the centerline of the 


flume, and with the following uniform flow characteristics established: | 
Yo = 1.34 in. U = 6.28 fps 
R = 58,500 (Reynolds No.) F = 3.31 (Foude No.) 


Measurements were made using the improved instrument, Gage 2, in the wake 
of the cylinder within the region enclosed by the standing waves emanating 
from the body and at three stations downstream from the centerline of the 
cylinder, at x/d = 40, 50, and 70. At each station the velocity traverses were 
made near the centerline of the flume and at a fixed distance of y = 0.84 in. 
from the channel bottom. 

It must be remembered that the wake investigated is not that of a cylinder 
in two-dimensional flow, but is also subject to shear conditions imposed by 
the channel bottom. 


Mean Velocity Distribution 


Figure 14 represents graphically the results obtained in the temporal | 
mean velocity phase of the wake study and is plotted as the ratio of the local 4 
temporal mean velocity U to the upstream mean velocity, Up computed at a 
relative depth y/yo = 0.63 vs. the ratio of the perpendicular distance from if 
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the centerline of the wake to the cylinder diameter, z/d. As a basis of com- 
parison at different stations downstream from the test cylinder, two equations 


were employed; the first was a function of the form of the Gaussian normal 
distribution: 


Us - K (2/d)* 
U, 


where uz = 1.06 U,- U 


Equation (13) showed good agreement with experimental data at the three 
downstream stations. The second expression was one derived by 
Schlichting.(9) 


« 
where: 


2b = the width of the wake x = distance measured downstream 
L/b from the cylinder in the direction 
of flow 


4 = Prandtl mixing length z = distance measured perpendicular 


Cy = coefficient of drag to and from the centerline of the 


d = diameter of cylinder wale 


U, = upstream mean velocity at distance 


U = temporal mean velocity y from bottom 


in wake 


The quantity 8 was evaluated experimentally; one set based on assump- t 
tions relating the mixing length to the macro-scale in a longitudinal direction 


and an additional set listed as Bexp. modified to agree with the data z/d = 0. 
These are: 


x/d B Bexp. 


40 0.144 0.174 
50 0.219 0.242 
70 0.245 0.245 8 


These values are quite close to the constant found experimentally by 
Schlichting; 8 = 0.207. Reference (4) contains a detailed discussion of this é 
phase of the investigation. 
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Figure 15 illustrates graphically the variation of the ratio of the root mean 
square of the velocity fluctuations to the mean stream velocity, Yur /Uo ‘ 


with distance z/d. A decrease of turbulence intensity with distance down- 
stream along with the transverse variations is indicated. The general trend 
of these results is in agreement with experiments of Townsend(11) and of 
Rdshko(12) performed with negligible boundary layer effects in air streams 
at subsonic velocities. The locations of intensity maxima observed by these 
investigators are also presented for comparison. At x/d = 70 the turbulence 
intensity has practically returned to the value of the undisturbed flow up- 
stream at the equivalent relative depth. 


Auto-Correlation Studies 


The auto-correlation curves shown in Figures 16, 17, and 18 for the 
velocity fluctuations are for points on or near the centerline of the wake at 
x/d = 40, 50, 70, respectively. They were obtained in accordance with the 
method described previously after the effect of the hydrostatic head and the 
mean velocity had been removed from the signal. 

The exponential function, Equation 9, fitted these curves quite well away 
from the origin as can be seefi from the figures. In the case of the auto- 
correlation function at x/d = 70, it was necessary to first subtract a cosine 
function of about 3.6 cps from the curve to obtain a better fit. This fre- 
quency is of the same order of magnitude as that experienced in the pre- 
viously menfioned uniform flow imvestigations. The periodicities at the 
natural frequency of the gage (140 cps) are barely perceptible indicating that 
the turbulence gage is capable of measuring accurately the major portion of 
the spectrum of turbulent energy in the wake. 

The scale of macro-turbulence, Ly, was determined by the procedure 


described previously and is presented in Figure 19 as the variation of (Ly/d)? 
with distance downstream on or near the centerline of the wake. 

A marked increase in Ly with distance downstream is seen from this 
figure, which may be explained as follows: Near the cylinder the body tends 
to break up the eddy pattern in the approaching flow and substantially reduces 
the size of the macro-scale below its free stream value. With distance down- 
stream the smaller eddies produced by the obstacle would tend to dissipate 
more rapidly than the larger eddies, thereby producing an increase in L, with 
x. Since the value of L, at the farthest station downstream is of the same 


order of magnitude as values obtained with Gage 2 in fully developed turbulent 
flow at comparable Reynolds numbers, a general return to upstream turbu- 
lence conditions is indicated. This same trend has been observed by other 
investigators in the flow behind grids. A curve of the form: 


(16) 


was fitted to the experimental data, showing the dependence of L, on x)/ 2 in 
the core of the wake. The dependence of lateral macro-scale on x1/2 has 
been found to hold downstream from a grid where the flow is known to be 
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isotropic.(13,14) Thus the longitudinal scale would also vary in the same 
manner. In all these considerations the vorticity generated by the interaction 
of the cylinder with the non-uniform flow of approach in the channel has been 
ignored. It would be of interest to investigate the entire problem in more de- 
tail. The present aim was confined however to establishing the general 
validity of results obtained with the instrument developed. 


CONC LUSIONS 


The studies presented indicate that an instrument consisting of a total head 
tube-pressure cell combination may be designed with sufficient frequency 
response to measure essential turbulence characteristics. The character- 
istics of a number of turbulent flows have been determined and have been 
found in essential agreement with comparable results obtained with the 
generally accepted hot-wire technique in airstreams. 

In general the following advantages may be stated: 


1. The instrument is suited for use in water, maintains its calibration and 
is mechanically simple. 


2. It may be built for laboratory as well as field use and for cltan water 
as well as for liquids with suspensions. 


3. Its mechanical features in relation to its natural frequency response 
are readily determined by a one-dimensional vibrational analysis. 


4. The instrument records with sufficient accuracy the instantaneous 
total head at a given point in the flow in the direction of the mean motion. 


5. From these records the turbulent fluctuations of the velocity in that 
direction can be determined with satisfactory accuracy. The error in the 
root mean square value of the turbulent component being normally less than 
5%. Auto and cross-correlations, the mean intensity spectra and scales of 
the macro-turbulence may be derived also. 


The disadvantages are: 


1. Turbulent pressure fluctuations are included in the record and cannot be 
extracted. 


2. The temporal mean velocity must be large as compared to the fluctuat- 


ing component. The error due to pressure fluctuations increases as the mean 
velocity decreases. 


3. Transverse velocity fluctuations cannot be determined at present with 
the instrument. 


4. The micro-scale of the turbulence apparently cannot be determined with 
sufficient accuracy. 


The claims made for the instrument are obviously not in any way tied to 
the present mechanical and electronic assembly, but will apply to any similar 
mechanical system. On balance it may be stated that the instrument provides 
a simple but restricted approach to turbulence measurements in liquid flow 
and that with adequate development of the electronic components a great 
many problems of liquid flow may be explored with respect to their most 


\ 
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important internal flow mechanics. In the attempt to verify the performance 
of the instrument a beginning has been made in that direction with the turbu- 
lence studies presented above. 


ACKNOWLEDG MENT 


The studies summarized above were conducted at the Hydrodynamics 
Laboratory of the Massachusetts Institute of Technology during the period of 
1949-1955 under the sponsorship of the David Taylor Model Basin and of the 
Office of Naval Research, Mechanics Branch, under Contract No. N5-ori- 
07874. Special assistance was given by the Digital Computer Laboratory of 
M.I.T. in the determination of auto-correlation factors and of the mean 
intensity spectra. 


REFERENCES 


. Macovsky, M. S., “The Measurement of Turbulence in Water,” David 
Taylor Model Basin Report No. 670, October, 1948. 


. Lowry, E. M., Jr., “Turbulence Studies from Recordings of Instantaneous 
Pitot Tube Pressures,” Master of Science Thesis, Massachusetts 
Institute of Technology, 1951. 


. Laufer, J., “Investigation of Turbulent Flow in a Two-Dimensional 
Channel,” NACA Report 1053, 1951. 


. Ippen, A. T.; Tankin, R. S. and Raichlen, F., “Turbulence Measurements 
in Free Surface Flow with an Impact Tube - Pressure Transducer Com- 
bination,” Technical Report No. 20, Hydrodynamics Laboratory, Depart- 
ment of Civil and Sanitary Engineering, Massachusetts Institute of 
Technology, July, 1955. 


. Lawrence, J. P., “Reduction of Turbulence Data with an Electrodyna- 
mometer,” Master of Science Thesis, Massachusetts Institute of 
Technology, 1953. 


. Taylor, G. I., “Statistical Theory of Turbulence,” Proc. Roy. Soc. of 
London, Series A, Vol. 156, 1935. 


. Dryden, H. L., “A Review of the Statistical Theory of Turbulence,” 
Quart, Appl. Math., Vol. 1, 1943. 

. Townsend, A. A., “The Structure of the Turbulent Boundary Layer,” 
Cambridge Philosophical Society Proceedings, Vol. 47, 1951. 

. Schlichting, H., “Lecture Series ‘Boundary Layer Theory’ Part I Laminar 


Flows and Part II Turbulent Flows,” NACA Tech. Mem. 1217 and 1218, 
1949. 


. Rouse, H., Engineering Hydraulics, John Wiley and Sons, Inc., New York, 
1950. 


. Townsend, A. A., “Measurements in the Turbulent Wake of a Cylinder,” 
Proc. Roy. Soc. of London, Series A, Vol. 190, 1947. 


| 

| 


ASCE IPPEN - RAICHLEN 1392-17 


12. 


13. 


14. 


Roshko, A., “On the Development of Turbulent Wakes from Vortex 
Streets,” NACA Tech. Note 2913, 1953. 


Batchelor, G. K. and Townsend, A. A., “The Decay of Turbulence in the 
Initial Period,” Proc. Roy. Soc. of London, Series A, Vol. 193, 1948. 


Batchelor, G. K. and Townsend, A. A., “The Decay of Turbulence in the 
Final Period,” Proc. Roy. Soc. of London, Series A, Vol. 194, 1949. 


| 
‘| 
| 
| \ 
| 
| 
| i 
| 
| 
4 


T e8ep Jo ydeusozoud 


(taqueyo 
einss add ) 
TT 9m 
w6erydetp 


——  aporzoata 


| 
1392-18 HY 5 October, 1957 | 

| 

| 

N 
| 


—— Pasnsoaw ajijoud > uny o 
{s uny x 
80 90 v0 


800 vo'o 


009 
4ajn07 


< 

=) 


Hid3Q 3AILV 134 
H1ld30 3AILV 134 


© 
° 


wo 


=u 


“kK /K 


ASCE 1392-19 
| 
| 
| 
: 
| | 
| 
° 


*@TZZON JO UT 


90 vO 20 O 90 v0 20 O 90 O 8090 ZO 20S 4ad 
31ZZON 


v 44 - JONVLSIO 


S600'0*S 
NNY 


Olbved 
8 NNY 


| 
1392-20 HY 5 October, 1957 ‘| 
| 
| 
| 
I] | 
| qa)" | | 
| | 
| | | | 
| | | | 
| F | 
| 00-0} | | 
i it 
|| 
| 
| \ | 
| 
| 
| 


= 
< 

z 


*srequny sptousey sno 


*yzdeg Jo 
SB SATPETSY 


[BUION 03 Aouenb 


ooo 012 —--—¢I 
----vol spucces o2 50 
| 000‘o¢ 6 | uo 
‘ON SOTONAZ nny 


| — — 


(2606) yuod yaeyo 
* 
2606 o 


Weed 
$4) 
| 94200*S | 


| 


| 


2 


Seen 


on 


> 


| 
| 1392-21 
| 
| 
xc | 
le 
| | 
418 
|" 
| a 
| | ~ 1° 
= Oo o ° ° 1 
| 4 —+ 3 + + - 
| = 
| | 
| | > | 2 
Je 
| 
| 
| 


abbey 
abey patquessy jo 
ainsserg jo werbetq 


*z pats Tpow 


October, 1957 


2 39Vv9 
1739 
3YuNSS3ud 


puy 
veemieg 
409 100° 


wBosydojg Asopuocses 
sequoyd bujdwog 
Bujy 
seqqny 
sbujy joes 


VILE 


I 
1392-22 ars 
| 
4 | 
Re: 
| 
> | 


IPPEN - RAICHLEN 


1392-23 


VI 
Gage 1 


Free stream signal 


it 
=== == 


Signal within 
boundary layer 


Gage 2 


Free stream and 
time signal 


Figure 9. Samples of signal records 
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Figure 1h. Mean Velocity Distribution in Wake 
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Figure 15. Turbulent Intensity in Wake 
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ABSTRACT 


Measurements of turbulence intensity at five stations and shear at one 
station in the boundary layer on the inside of a7 1/2 degree conical diffuser 
are reported. Comparison of the results with those reported for non- 
pressure gradient boundary layers indicates significant divergence from the 
universal inner region concept. The microscale of turbulence appears to be 
remarkably constant although other turbulence quantities are greatly in- 
creased as the flow develops in the adverse pressure gradient. 


Turbulence in a Diffuser Boundary Layer 


The state of knowledge regarding the exact nature of the flow in a turbulent 
boundary layer developed under the influence of an adverse pressure gradient 
is best described as incomplete. Information is required, both on the mean 
flow characteristics and on the nature of turbulence and turbulence processes 
in the developing boundary layer. Mean flow data have been reported by vari- 
ous investigators; nevertheless, a recent comprehensive study of the situation 
for plane two dimensional flows by D. Ross(1) has indicated the need for even 
more information. If, in spite of the relative wealth of literature available on 
this phase of the subject, it can be described as “incomplete,” then the state 
of knowledge regarding the turbulent aspects of the situation can only be 
described as “fragmentary.” This interpretation has been stated in slightly 


Note: Discussion open until March 1, 1958. Paper 1393 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 5, October, 1957. 

*Prof. of Theoretical and Applied Mechanics, Univ. of Illinois, Urbana, IJ. 

**Project Leader, Hydraulics Group, Covington Research Center, West 

Virginia Pulp and Paper Company, Covington, Va. Formerly Asst. Prof. 
of Eng. Research, Ordnance Research Lab., Pennsylvania State Univ., 
University Park, Pa. 


1393-1 


HY 5 
| 
| 
| 
| 
| 
| 
| — 
| 


1393-2 HY 5 October, 1957 


stronger terms by G. B. Schubauer and P. S. Klebanoff(2) who commented, 
“The real problem is the understanding of the mechanics of turbulent shear 
flow under the action of pressure gradient. The solution of this problem de- 
pends upon the understanding of the mechanics of turbulence and in this only 
rudimentary beginnings have been made.” 

The random velocity fluctuations known as turbulence must be character- 
ized with the aid of certain statistical values, the first of which was introduced 
by O. Reynolds. Thus the flow velocity at any time is defined by the relation, 
velocity = U +u inthe x direction. Here U is the temporal mean velocity 
which varies with location and u is the fluctuating component. This random- 
ly fluctuating quantity is characterized by u', its root-mean-square (rms) 
value, known as the turbulence intensity. From a mean velocity consideration, 
turbulent flows with small (as in pipes and channels) or zero (as in boundary 
layers) pressure gradients are similar, particularly near the wall. Within a 
wall region 10 or 15 percent of the pipe radius or boundary-layer thickness, 
the “Law of the Wall” applies and the velocity distribution is given by a 
logarithmic relation. Of course, there is a laminar film at the wall if it is 
smooth. In the outer portion of the flow, away from the wall, other velocity 
relations are more appropriate. 

The effect of adverse pressure gradients on turbulent boundary-layer type 
flows have been clearly indicated by Ross(1) from a mean flow standpoint. 
The “Law of the Wall” associated with the wall shear stress still applies in 
the inner region. The flow in this region is singularly unaware of the pres- 
sure gradient. In the outer region the flow is greatly modified through the 
cumulative effects of the pressure gradient. The velocity profiles become 
sharper with inflection points sometimes appearing as the wall shear stress 
becomes smaller. Eventually, this becomes zero and the flow separates from 
the wall. 

The status of our knowledge about the turbulence in flows with negligible 
pressure gradients (including pipes and channels) has been well summarized 
by G. B. Schubauer.(3) He shows that, except for an intermittency which oc- 
curs in the outer part of boundary-layer flows, the turbulence in pipes, (4) 
channels,(5) and boundary -layers(6) is in many respects the same. Thus if 


the intensities are expressed in terms of the shear velocity U.= TH /p 


(where Ty = wall shear stress and p = fluid density) as a function of the 
friction wall distance paramenter y Ur /v (where y = wall distance 


and v = fluid kinematic viscosity) unique relations result. Understanding of 
the turbulence in these flows is elucidated through a detailed study of the 
energy processes. Schubauer also shows the similarity for the various flow 
cases of the several terms in the turbulent energy balance equation, which 
may be written in simplified form 


*D+C+W=0 


where 


P,. indicates the production of turbulent energy from the mean motion 
D, the kinetic energy and pressure energy diffusion across the flow 
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C, the convection of turbulent energy by the mean motion, and W, the 
dissipation of the turbulent energy. 


The first of these terms may be expressed in dimensionless form as 


OU 
P y 


which is simply calculated. Part of the mean flow energy is directly dissi- 
pated; the dimensionless direct-viscous-dissipation rate is given by 


2 2 
Wy 


When integrated across the boundary-layer or conduit (P, + Wy) account for 


the loss of kinetic energy of the mean flow - some of it going directly to heat 
through the action of viscosity and the remainder going into the production of 
turbulence energy. Thus we have a connection between mean flow conditions 
and turbulence. Some of the turbulence terms involved in the convection and 
dissipation have been measured so that the trends for C and W are known with 
some reliability. The diffusion terms are found by balancing the equation and, 
as summarized in Schubauer’s paper, the diffusion is significant. 

It has already been remarked that information on the turbulence in adverse- 
pressure gradient flows is “fragmentary.” G. B. Schubauer and P. S. 
Klebanoff(2) at the National Bureau of Standards, B. G. Newman(7) in 
Australia and V. A. Sanborn and R. J. Slogar(8) at the NACA have contributed 
the major portion of the literature available on turbulence in adverse pres- 
sure gradients. All of these references describe studies of the plane two- 
dimensional boundary layer and the first two were concerned primarily with 
the separation phenomena. 

A study of boundary layers in axisymmetric flow was conducted at the 
Pennsylvania State University from 1946 to 1954. Recently E. M. Uram(9) 
along with J. W. Holl and J. M. Robertson(10) have investigated the mean 
flow parameters for an axismmetric boundary layer encountered in a 
circular diffuser. Near the conclusion of these studies, a two channel, 
constant-temperature hot-wire anemometer became available for a brief 
series of turbulence measurements in this facility. Time commitments did 
not permit a thorough, detailed study of the turbulent processes in the 
developing boundary layer. Nevertheless, it is hoped that the relatively 
modest amount of information obtained during the course of this investigation 
and presented in this paper will contribute toward the advancement of the 
knowledge regarding turbulent boundary-layer flow phenomena under the ac- 
tion of an adverse pressure gradient. 


Apparatus and Test Procedure 


Flow Circuit 


The turbulence measurements were made in the developing boundary layer 
of a7 1/2 degree (total angle) conical diffuser which was also the test facility 
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for the mean velocity studies(9,10,11) at Penn State. The particular configura- 
tion corresponds to the Phase II arrangement reported in the recent papers 
and was essentially an Eiffel type air tunnel. The diffuser had a 6-inch en- 
trance diameter and increased to a 11 1/4 inch diameter where it discharged 
freely. The test set up consisted of a 24-inch axial flow fan followed by re- 
ducing section, screens and a cast aluminum nozzle. Six diameters of 

straight conduit separated the nozzle from the start of the diffuser. This 
promoted the growth of a reasonably thick boundary layer before the start of 
the adverse pressure gradient. 

Holes suitable for the insertion and traversing of instrumentation were 
provided at a number of stations along the diffuser. The location of those 
stations at which turbulence measurements were made is shown in Figure 1 
along with a dimensional description of the diffuser. A plot of the pressure 
distribution obtained from the previous studies(9,10) is included. As in other 
diffuser flows, the initial pressure gradient is comparatively steep. It is 
seen from this figure that the first measurement station was at or slightly 
ahead of the start of the strong pressure gradient. 

The physical condition of the boundary layer entering the diffuser is de- 
fined by the velocity outside the boundary layer U, = 166 fps and the 
boundary-layer Reynolds number based on momentum thickness, Rg = 6900. 
This was determined at Station 1 in the mean flow investigation.(9,10) 

Numerous accurately machined, flanged construction joints in the diffuser 
section permitted disassembly of the facility to insert the hot-wire probe in 
the traversing apparatus. For each location, the zero reading of the traverser 
was adjusted to place the plane of the hot wire at the wall of the diffuser. In 
this manner, the reading of the traverser was actually the distance of the 
probe from the diffuser wall. It is estimated that the cumulative radial posi- 
tional error arising from the initial installation of the probe in the traverser, 
and the error in the traverser screw, amounted to no more than 0.02 inches. 
A Pitot tube was positioned in the same measurement plane as the hot wire 
through an instrumentation hole in the opposite wall. The dynamic pressure, 
derived from the difference in the Pitot impact pressure and the static pres- 
sure obtained from wall piezometers, was used to check the calibration of the 
hot wire probe and make suitable corrections for any temperature induced 
drift in the anemometer. The mean flow rate was monitered during the tests 
by measuring the pressure drop across the nozzle. It remained steady and 


corresponded very closely to the flow rates developed during the other 
studies (9,10). 


Hot Wire Anemometer 


This paper is not primarily concerned with the development and design of 
hot wire anemometers; therefore, the fine points of the subject will not be 
discussed here. To those who wish to review this field, we refer to the ex- 
cellent survey of hot wire anemometry by R. D. Cooper and M. P. Tulin.(12) 
The particular instrument used in this investigation was constructed by 
the Hubbard - Ling Company of Iowa City, Iowa. It is a constant temperature 
system designed to overcome the two major inherent defects of the hot wire 
as a velocity measuring device, namely the thermal capacity of the wire and 
its non-linear conversion of velocity to electrical signals. The frequency 
response of this hot wire anemometer is + 1 percent, 0-10 KC, down 10 per- 
cent at 25 KC, and the equivalent noise level (relative to mean velocity) is 
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0.1 percent or less for velocities over 20 fps. 
The two channel feature of this instrument along with the signal manipula- 
tion possible in the metering circuit enabled the determination of the follow- 
ing quantities with relative ease: U, the mean velocity in the axial direction, * 
u', and longitudinal component (rms value) of turbulence and A ., the micro- 
scale of turbulence, all using a single wire placed normal to the mean flow 
direction. The use of an X or crossed wire probe enabled the determination 


of v', the rms orthogonal component of turbulence and uv the Reynolds shear 
stress in the xy plane. 
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Test Procedure 


Calibration of the hot wire probe was conducted prior to each boundary 
layer traverse by placing the wire outside the boundary layer, and establish- 
ing a mean velocity calibration curve. In progressing down the diffuser, the 
maximum velocity, Uj at any station dropped from 166 fps to 98 fps. Lower 
velocities were obtained for the calibration process by placing a cloth over 
the entrance grill of the fan. This throttled the flow very effectively without 
disturbance. 

Following calibration of the probe, measurements were taken at frequent 
intervals across the boundary layer. As the wall was approached, the inter- 
vals between adjacent measurement stations were systematically reduced 
from 1/4 to 1/100 inches. The probe was then returned to the initial position 
outside the boundary layer; during the return, a number of points were re- 
peated as a check on the results. In several instances, the traverses were 
repeated on different days as a further check on reliability. 


Results 


The change in the mean velocity profile as the boundary layer grows along 
the diffuser is illustrated in Figure 2. Here the local mean velocities have 
been reduced to dimensionless numbers in terms of the maximum velocity 
measured at each station. All of the data were derived from hot wire 
measurements. To indicate agreement between hot wire and impact tube 
measurements, a curve from Uram’s data (at station 6) is also presented. 

It can be seen that the hot wire has followed the impact tube measurements 
quite faithfully except for the region near the wall. Here the difference be- 
tween the two sources of data are recognizable but are in the wrong direction. 
Normally, one would expect the impact tube to read high because of the effect 
of turbulence on the impact pressure. The multiple sets of data, obtained on 
different days for Stations 3 and 12 show some shift. However, it amounts to 
only a few of percent, which is within the repeatability limits of the hot wire 
anemometer. In general, the validity of the hot wire mean flow measure- 
ments is established. 

The mean velocity profiles have also been analyzed to determine the dis- 
placement and momentum thicknesses as given by the following definitions: 


*The mean velocity in the y direction, V is sensibly zero, since the maxi- 
mum angle of the streamlines relative to the diffuser axis is 3.75 degrees. 
In the absence of secondary motion, the circumferential component W is 

zero. 
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These are actually two-dimensional definitions rather than the three- 
dimensional definitions which should apply.(10) However, for comparative 
purposes, the differences are irrelevant. Comparison of the hot-wire deter- 
mined momentum and displacement thicknesses with those reported by 

E. M. Uram(9) are presented in Figure 3. In general, the agreement is good. 
Also included in this figure is the common shape parameter H = 6*/@ . 

As in the case of other boundary layers, H is seen to increase to a value of 
about 2.8 as the end of the diffuser (and separation) are approached.(10) 
These parameters along with the pressure distribution serve to describe the 
mean flow conditions under which the turbulent boundary layer developed. 

The distribution of the relative turbulence intensity at the five measure- 
ment stations along the diffuser is illustrated in Figure 4. Turbulence inten- 
sity is normally nondimensionallized through division by a reference velocity, 
in this case the local velocity U. Measurements at Stations 3 and 12 made on 
different days are indicated by modifications of the basic symbols. Variation 
from day to day does not seem particularly large. The greatest scatter 
coincides with the highest turbulence levels, particularly in the case of the 
single traverse at Station 17. Other investigators have found considerable 
scatter during measurements of high intensity turbulence.(12) 

Curves derived from Klebanoff’s(6) investigation of the zero pressure 
gradient boundary layer are included in both Figures 2 and 4, for comparison 
with the results found at Station 1 which was essentially in a zero pressure 
gradient and at the same Reynolds number. These comparisons serve to 
show that the velocity distribution and turbulence in the boundary layer under 


study was quite normal prior to its exposure to an adverse pressure gradient. 


They also indicate a fair degree of reliability in the hot-wire measurements. 

The curves of Figure 4 indicate relative intensities which do not permit 
comparison between sections or even across the boundary layer at the same 
section. For such comparisons, the absolute intensity u' is more suitable. 
Schubauer and Klebanoff(2) refer intensities to the velocity outside of the 
boundary layer in the flow just before the start of the adverse pressure 
gradient. The intensities measured in this present study were similarly 
converted to dimensionless rations u'/U49 where the initial reference 
velocity Uj9 was taken at Station 1. From smooth curves drawn through the 
data, the contour plot of Figure 5 was obtained, showing the variation in 
turbulence intensity in the boundary layer. This shows a definite peak in the 
u'/ Uyo values; this peak shifting outward as the flow progressed. The re- 
sults are compatible with data presented by Schubauer and Klebanoff except 
the present data seem more internally consistent. 

Shear stress and v' measurements were made only at Station 12. The 
shear stress results are shown in Figure 6. Considerable scatter is appar- 
ent in the data, especially in the region of greatest shear. Again, qualitative 
similarity with the observations of Schubauer and Klebanoff is found. 

The longitudinal microscale or dissipation length can be expressed by the 
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following relation: 


u'U 


Experimentally, the square root term was found from the hot wire signal with 
an electronic differentiating current and A, calculated with the aid of the 
other measured quantities. 

If the turbulence were isotropic, the turbulent-energy dissipation rate 
would be 


Various researchers have found that when this expression is integrated across 
the flow, the result may be low by a factor of 2.5 in the simpler flow case of 
negligible pressure gradient. However, even for such non-isotropic cases, 

Ax does give an indication of the size of the smallest eddies. The results of 
the measurements of this length are given in Figure 7. Smooth curves rather 
than individual data points are presented for the sake of clarity. In the case 
of Station 3 individual data points are given to indicate the scatter. 


Discussion 


The most noticeable feature of the relative intensity curves of Figure 4, 
other than the rapid increase in this quantity as the flow proceeds, is the 
occurance of a peak for Stations 6 and 12. Zero pressure gradient boundary 
layers also display a peak; however, this occurs at wall distances corres- 
ponding to the thickness of the laminar film or sub-layer. In the case of 
Stations 6 and 12, the distance from the peak to the wall is at least several 
orders of magnitude greater. If this peak is a characteristic feature of ad- 
verse pressure gradient boundary layers, its absence at Station 17 is quite 
unexplainable. For the plane two-dimensional boundary layer with an ad- 
verse pressure gradient, Schubauer and Kelbanoff(2) do not seem to have 
found such a peak. 

In the introduction it is pointed out that an inner region exists in boundary 
layer flows which is independent of external conditions. This inner region is 
about 10 percent of the boundary layer thickness; this was verified by the 
mean flow measurements(9,10) in this boundary layer. The inner concept 
would lead one to expect essential identity of the turbulence in this region. 

The measured intensities presented as u'/Uy and v'/U- in comparison 
with the pipe results of Lauffer(4) are shown in Figure 8. In the case of the 
u' measurements, there is considerable divergence between the values ob- 
tained with the X probe and the single wire*. However, this uncertainty is 
insufficient to invalidate the general results shown. It is obvious from this 


*The single wire points shown in Figure 8 do not represent actual measure- 
ments as in Figure 4, but average values from the measurements. 
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figure that in the inner region (out to about 46/10) the turbulence intensities 
are much greater than found in the pipe flow. In the laminar film (thickness 
= 64), these two quantities are probably identical as suggested by the dashed 
lines. In the outer region, the turbulence intensities reach quite large magni- 
tudes compared to those found in other boundary-layer-type flows. This 
outer effect is not unexpected since the flow there is the result of upstream 
generated turbulence.(13) 

Additional evidence concerning the high turbulence level in the inner region 
was obtained through calculation of the turbulence production rate P;. This 
quantity was determined from mean curves drawn through the shear stress 
and velocity measurements at Station 12. The direct viscous dissipation rate 
WH was also obtainable. The results of these calculations are shown in 
Figure 9. Comparison with results for flat plates and pipes(3) show that in 
the inner region the turbulence production is 5 to 30 times greater while the 
direct dissipation rate seems to be about the same (or at least insignificant). 
Again, in the laminar film, it appears that had measurements been made, the 
values would have been identical. Just as in the case of the intensity, there 
is a large increase in the production of turbulence in the outer boundary layer 
region. 

Proper understanding of the occurances in boundary layers such as that 
studied naturally would require evaluation of the other terms in the turbulence 
energy equation. The necessary measurements to evaluate these should have 
been made, but the extent of this study was limited by other factors. The high 
turbulence production rate in the inner region for the adverse pressure 
gradient case does not appear to have been commented on by other investi- 
gators. In retrospect, it seems hardly surprising that this was found since 
the shear stress in these flows rapidly increases to many times the wall 
value as the wall distance increases.(13) 

The measured shear stress distribution at Station 12 and shown in Figure 6 
differs greatly from that found in pipes and zero-pressure gradient boundary 
layers. In such cases, the shear has a maximum at the wall and decreases 
monotonically to zero at the outer edge. Peaked shear distribution curves in 
adverse pressure gradients have been noted previously.(2,7) An explanation 
for this in terms of the “history” of the flow has been presented elsewhere.(13) 
At the wall, Or/@y equals dp/dx from the momentum equation, written for 
the wall streamline. As the pressure gradient is positive, the shear has to 
increase as one departs from the wall. A theory was presented in 1950 ex- 
plaining the variation in the wall shear stress, as affected by the local ad- 
verse pressure gradient, wall shear stress and boundary layer thickness, 
and the boundary layer thickness and wall shear stress at the start of the 
pressure gradient.(13) Figure 6 also presents the calculated shear stress 
distribution based on this theory. The wall shear stress coefficients used 
for this calculation were based on pitot tube measurements.(10) The agree- 
ment between theory and experiment is seen to be good. The theory gives the 
total shear, viscous plus turbulent, while the hot wire measurements yield 
only turbulent shear. This consideration does not affect the comparison of 
Figure 6 for at a wall distance of 0.1 inches, the viscous shear has been esti- 
mated to be 0.2 percent of the measured shear. 

The dissipation length A, results are surprising in that they indicate little 
variation between stations. It would appear that except near the wall, the 
smallest eddies had a size of about 0.01 feet. Other recent boundary layer 
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measurements also indicate little variation.(14,15) Lauffer’s results shown 
in the Figure 7 range from 0.01 to 0.02 feet. As already noted, if the turbu- 


lence were isotropic, the dissipation of turbulence into heat would be given 
by the relation 


' 2 
= 2 
Wo dE = 


With the essential constancy in A the isotropic dissipation is proportional to 
u'2, Inspection of Figure 5 indicates that turbulent energy dissipation in- 

creases rapidly and has its maximum about the middle of the boundary layer 
at the downstream stations. One might conjecture that the rate of actual dis- 
sipation increases in a similar fashion. 


CONCLUSIONS 


Although the turbulence studies were not as extensive as desired and the 
precision of the data could be improved on, the results are felt to be suffi- 
ciently adequate to justify the following conclusions: 


1. The turbulence, its rate of production, and its rate of dissipation in an 
adverse pressure gradient are greatly in excess of similar quantities for 
zero-pressure-gradient boundary layers. 


2. The existence of an inner wall region having universal characteristics, 
although well established for mean flow conditions, has not been found for 
turbulence. Thus, the turbulence and its rate of production are far in excess 
of the zero-pressure-gradient case, even near the wall. 


3. The shear stress variation across the boundary layer is reasonably 
well predicted by the history theory. 


4. The longitudinal microscale of turbulence A, remains remarkably 
constant across and along the developing boundary layer. 
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Fig.5 Contour Plot of Turbulence Intensity 
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Fig.7 Turbulence Micro-Scale or Dissipation Length 
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INTRODUCTION 


The experimental investigation of the component velocity fluctuations in a 
turbulent fluid requires a method of velocity measurement of high sensitivity, 
accurate frequency response, and one which will preferably distinguish 
between motion in various directions. The hot-wire anemometer adequately 
satisfies these criteria in the case of gas flows but has proven unsatisfactory 
in the study of the details of turbulent fluctuations in liquid media. It is the 
purpose of this paper to describe a method of component velocity fluctuation 
measurement in turbulent liquid flows based upon the induction of an electri- 
cal potential field in a continuous medium moving relative to a stationary uni- 
form magnetic field. An analysis of the theoretical basis of the proposed 
method will be presented, as well as a description of the experimental appli- 
cation to the case of the measurement of the component velocity fluctuations 
in the fully developed turbulent flow of water in a circular pipe. Attention will 
be focused upon the physical phenomena underlying the method of electro- 
magnetic induction which will serve to point up the general advantages and 
limitations of the principle in the experimental investigation of turbulent 
flows. The directions in which further development appear to be fruitful, will 
be indicated in a concluding section. 


Electromagnetic Induction in a Fluid in Motion 


The physical principle upon which the experimental method is based is 
Faraday’s law of electromagnetic induction which indicates that electromo- 
tive forces are induced in a liquid conductor moving relative to a magnetic 
field, the electric field strength at any point being simply related to the mag- 
netic field strength and the velocity. This fact may be used for investigating 
the distribution of velocities in a moving liquid, for this distribution can be de- 


Note: Discussion open until March 1, 1958. Paper 1394 is part of the copyrighted 
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duced from observations on the e.m.f.’s induced at different points in the 
liquid by a known magnetic field. 
Let: 


> 

V =the velocity vector at any point in a fluid 
> 

B = the magnetic field intensity 


> 
E = the electric field intensity or induced e.m.f. per unit length 
= the induced current 


¢ = the electrical conductivity of the fluid 


9 = the electrical potential 


Faraday’s Law of Induction can be formulated as follows. Every changing 
magnetic field has associated with it an electric field E; the line integral of 
E around any closed circuit is related to the time rate of change of the flux 
across the circuit by the basic relation, 


= \\B-aA 
bE. ai dt |\ (1) 


wheredA is the element of area of a surface capping the circuit. 
The total magnetic flux through a given circuit appearing as the integral 
in the right hand side of equation (1) can change with time for two reasons: 
a) It can change because of changes in the external field with time. 
(“transformer action”) 
b) It can change because of motion of the circuit itself or parts of the 
circuit. (“motional induction”) 


In the case of a flow steady in time relative to a constant external magnetic 
field B, only the latter form of induction is operative. If one develops the 
form of the right hand side of (1) for a surface itself in motion with the fluid 
and applies Ohms Law to relate the current to the electric field, there results 
a basic relation for the e.m.f. and the motional induction, namely, 


= 
BxV + @) 
The product of magnetic field and velocity is a vector product following the 


familiar “Flemings right hand rule”. 
Since E is the gradient of the potential P ,» one may also write 


= BxV + Ve 


In the case of steady flows the continuity equation for electric charge is 


(3) 


div = O 


so that taking the divergence of both sides of (3) yields 


(4) 


= div (VxB) 


__ 

= 
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This Poisson equation indicates that for a given magnetic field the potential 
gradient is uniquely a function of the velocity field and the boundary con- 
ditions of the flow. 

Under certain conditions the induced current flow is small or identically 
zero, so that from equation (3) we see that a measurement of the potential 
gradient at a point where the magnetic field is known is a direct measure of 
the local velocity. Such a method of velocity measurement has the advantages 
that: a) it is based upon a universal and fundamental physical principle, b) 
it is possible to distinguish between velocity components due to the vector 
nature of the basic relation, c) it yields an instantaneous measure of veloci- 


ties fluctuating in time, and d) it is insensitive to the physical properties of 
the fluid medium. 


Application of the Method of Electromagnetic Induction to 
Local Velocity Measurement 


The discussion of Faraday’s principle of electromagnetic induction applied 
to continuous media in motion relative to a fixed magnetic field, led to the 
fundamental equation for the local potential gradient in terms of the local 
values of the velocity, magnetic field strength and current density, namely 
equation (3). If the velocity distribution is known, it is always possible to 
compute the resultant electrical field, including the potential and current 
distributions by solving a Poisson equation with suitable boundary conditions. 
The central question for experimental application, however, is whether 
one can deduce the velocity field from a measurement of the distribution of 
potential gradient. Here one faces the difficulty that the current density in 
equation (3) is not known a priori, but is rather a function of the velocity field 
that one wishes to deduce. The local value of the current density appearing 
in equation (3) is moreover determined by the total velocity vector field, so 
that the local potential gradient is dependent upon the entire velocity pattern. 
It might be supposed, at first view, that the influence of induced currents 
could be reduced by increasing the conductivity of the fluid; but unfortunately 


it is the quantity > that appears in the fundamental equation, and this quanti- 


ty is determined solely by the irrotationality of the induction vector Bx V 

In the case of the measurement of mean velocity it is always possible to 
apply the method of electromagnetic induction to determine the velocity 
distribution experimentally. In such cases one can compute the velocity 
pattern from measured values of the potential gradient assuming no current 
flow and then with this assumed velocity distribution compute the induced 
current flux by the methods previously outlined. A process of successive 
approximation would then yield the correct velocity distribution. For the 
case of two-dimensional flow with the vorticity always in the direction of the 
applied magnetic field the induced current flow is automatically zero. The 
induced e.m.f. sets up space and surface changes whose electrostatic field 
exactly balances the induced e.m.f. 

The application of the method to the measurement of local turbulent ve- 
locity components raises more serious questions. Here there is no possibil- 
ity of utilizing the simplification of the induction equation for two-dimension- 
al flow, for the velocity fluctuations are conceded to be three dimensional 
regardless of the mean flow characteristics. In the case of turbulent fluctu- 
ations, moreover, one measures the mean square of the potertial gradients 
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and the expression for the latter deduced from the induction equation by 


Reynolds’ rules of averaging involve product correlations whose magnitude 
can not easily be estimated a priori. 


Consider, for example, the case of pipe flow with the axis as shown and 


B uniform in the direction 
of the 4 axis. Then from 


we have 


- b/g 


dy = -BU ~ by/¢ 


29 _ 
(7) 


Splitting the velocity components into mean and fluctuating parts, and the 
same for the current, i.e. 
U = u + u' 
. 
s + Ly 


and applying the Reynolds rules of averaging we have 


\2y) ) + 28 (at, su) 
= B(pau") - 
+u") B (vt + ui )+ ( 


The expression for the mean square of the fluctuating potential gradients are 
then, 


| 


| 


= 

grad = (VxB) - Ye (6) 

| 

‘| 

— 
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(10) 


It is clear from these relations that provided the presence of induced 
currents could be neglected that the r.m.s. values of the u and v components 
could be computed directly from the second and third equations from 
measured values of the r.m.s. potential gradients. But in the actual case the 
relations are complicated not only by the unknown current intensities but by 
the presence of the correlations between velocity and current fluctuations re- 
sulting from squaring the induction equation. 

Fortunately it is possible to obtain a direct experimental estimate of the 
intensity of the current fluctuations by utilizing the first of the above 
equations. This shows that the r.m.s. potential gradient measured in the 
direction of the magnetic field gives the r.m.s. current in the same direction. 
The results of this measurement and its interpretation will be discussed in 
subsequent sections. 


Experimental Program 


Flow System 


The experiments of this research were performed on water flowing in a 
2" 1.D. lucite pipe. As shown in a schematic diagram, Fig. 1, water was 
pumped from a reservoir to an entrance head tank I, then passed through a 
gravel-filled passage to another entrance head tank II. The passage was 
introduced specifically for the purpose of eliminating large scaleflow dis- 
turbance. In order to obtain a smooth flow in the 26' - 6" long pipe of 2" I.D., 
the pipe entrance was bell-shaped. The bell-mouth was connected to the 
head tank by a rubber diaphragm. The flexible connection was designed to 
eliminate the transmission of mechanical vibration between the entrance tank 
and the experimental pipe. The pipe was arranged horizontally and well 
supported by nine shock mounts. A heavy masonry support was also placed 
near the pipe entrance to prevent horizontal movement. The flow discharged 
from the lucite pipe to an exit tank then drained back to the reservoir, the 
flow rate being regulated by varying the field resistance of a D.C. motor 
which was driving the pump. The flow rate was determined gravimetrically. 


Magnet and Measuring Probe 


Magnet—The magnetic field was provided by a direct current electro- 
magnet with 5.75 inch diameter circular pole pieces and a 2-1/2 inch gap. 
The maximum field intensity was about 9000 gauss. The field strength was 
practically uniform over the measuring region. D.C. current was supplied 
by a 110V, 50 KWD.C. generator of the hydro-laboratory of the University 
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of California, Berkeley. In order to help eliminate the ripple voltage in the 
D.C. supply, a 60-micro-farad condenser was connected across the coils of 
the magnet. An ammeter was used to indicate the current passing through 
the coils of the magnet. The readings of the ammeter were calibrated 
against the readings of a standard gauss- meter, which was also used to 
check the uniformity of the field. 


Measuring Probe 


One electrode, two-electrode, and three-electrode probes were used in 
this investigation, the arrangement of the probe depending on the type of 
measurement. A micrometer support for the electrode probe was provided 
for traversing the pipe cross-section. A photograph of the measuring probe 
assembly is shown in Fig. 2. 

The probe was constructed from 0.0125 inch diameter insulated copper 
wires, enclosed in a 0.054 inch O.D. copper tube (the assembly was manu- 
factured by the Precision Tube Co., Phila., Pa.), bent into an L shape so that 
the wire tips were 3/4" upstream from the stem. The tips of the wires were 
separated by the proper distance and insulated by being painted with glyptal 
except at the tips. At the other end of the probe the electrode wires were 


soldered to a microphone connector which was fastened to the micrometer 
assembly. 


Electronic Circuits 


The signal from the probe was fed to the differential input of a Tektronix 
type 122 low level preamplier; the output signal of the preamplier was applied 
directly to a Dumont 304-A Cathode Ray Oscillograph. In addition, the instan- 
taneous signal of the preamplier was amplified once again by a single stage 
amplifier and fed to a squaring amplifier to obtain the time average square 
g" of the signal. 

Two complete channels were provided so that the two input signals from 
the probe could be measured at the same time. 
The block diagram of the electronic equipment is shown below. 


122 t 
INPUT | CROSS | * DU-MONT 
| PRODUCT | CATHODE RAY 
| OSCILLOGR APH 


Results of Measurement 


A sample of the oscilloscope record is shown in Fig. 3. The upper trace 
represents the basic amplifier and flow noise without the magnetic field. 
The lower one represents the measurement signals with the magnetic field. 
This sample was obtained by the input of a two electrode probe with 0.10" 
separation at a Reynolds number of 1.0 x 109. 

The results are presented as the square root of the direct readings of the 


time average signal by means of the cross-product instrument. The follow- 
ing measurements were made: 


| 


| 
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(1) RMS potential gradient in the direction normal to the mean flow as a 
function of radial position in the pipe. A two electrode probe was used 
for this type of measurement. This arrangement of the probe is shown 
in Fig. 4. The experimental results of a two electrode probe with 
0.025" separation (center to center) are shown in Fig. 4. The abscissa 
is the fraction of the radius, while the ordinate is the ratio of the read- 
ing of the potential gradient to a fractional velocity an.s. 79 /s Uy 
The u, is defined as 


U, = | Tun 
where To is the shear stress at the wall 
is the density 


$ is the friction factor 


Uw. is the average velocity of the pipe flow. 


(2) RMS potential gradient in the axial direction as a function of radial ? 
position in the pipe—two electrode probe with the arrangement shown 
in Fig. 5. 
(3) RMS potential gradient in the direction of the magnetic field as a 
function of a radial position in the pipe. Fig. 6 represents the 
measurement results of this type. 


Discussion of Results 


With the orientation of the uniform magnetic field perpendicular to the 
pipe axis as employed in these experiments it is only possible to obtain in- 
formation on the component velocities in the direction of the mean flow and 
perpendicular to the mean flow. The experimental information relevant to 
these velocity component measurements is summarized in the curves of 
figures 4 through 6. Due to the presence of fluctuating induced currents 
whose magnitudes can not be estimated a priori, it was thought preferable to 
present the data in terms of the directly measured root mean square values 
of the component electrical potential gradients rather than estimated root 
mean square values of the velocity components. For purposes of comparison 
the data of Laufer* on velocity component fluctuations in the flow of air 
through a circular pipe obtained by means of the hot wire anemometer have 
been plotted in Figs. 4 and 5. 

As has been previously discussed (see equation 10), provided that no in- 
duced current flow is present, the potential gradient normal, to the flow isa . 
direct measure of the axial velocity component and the potential gradient in ,; 
the direction of the mean flow represents the normal velocity component. 

Comparing the measurements of rms potential gradient divided by the uni- 
form and constant magnetic field intensity, to the hot wire velocity measure- 
ments of Laufer in air flow, we find first that the distribution of potential 
gradient as a function of radial position for both the axial and transverse 
cases follows very closely the corresponding component velocity fluctuation 
measurements obtained by Laufer and second that the magnitudes compared 


*J. Laufer, NACA Technical Note 2954, June 1953. 
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on a dimensionless basis by division by u, show a marked discrepancy. If 
the potential gradients per unit magnetic field intensity be interpreted as 
velocity components, the present set of measurements are consistently lower 
than those obtained by means of the hot-wire anemometer. We have conclud- 
ed that the difference is due to the presence of current flow such that the in- 
ductional e.m.f. given by the vector product of velocity and magnetic field is 
greater than the measured potential gradient by the magnitude of the induced 
current intensity. 

Fortunately it is possible to measure at least one component of the current 
vector, for by the basic induction equation we conclude that the potential 
gradient in the direction parallel to the magnetic field measures the current 
flow in the same direction. In Fig. 6 measurements of the r.m.s. current 
density in the direction of the magnetic field are plotted as a function of the 
radial position. From equation (10) the magnitude of r.m.s. B42 is equal 


to r.m.s. Lele It may be noted that this quantity is of the correct order of 


magnitude to bring the potential gradient measurements into agreement with 
the data of Laufer if one assumes that the remaining components of current 
fluctuation are of the same order of magnitude as that measured, and if the 
correlation between velocity and current fluctuations is of negligible order 
compared to the other terms of equation (10). Since these hypotheses could 
not be verified experimentally in the present investigation, the potential 
gradient data have been reported directly. 

From the foregoing discussions of the basic theory of the method of 
electromagnetic induction and the results of the experimental measurements 
it will be clear that in order to render the method quantitative for the 
measurement of fluctuating velocity components it is necessary to measure 
the three components of the current density. This can be accomplished pro- 
vided that the orientation of the uniform magnetic field can be varied at will, 
since the potential gradient parallel to the magnetic field measures directly 
the current flux in the same direction. For this reason and another to be 
mentioned shortly, attention is presently being directed to the development of 
an electrode probe assembly of which a small permanent magnet forms an 
integral part. Alnico magnets of small size and high flux are currently being 
investigated for use in a self-contained probe. In addition to the advantage of 
changing the field direction such a probe assembly offers the possibility of 
considerably extending the range of flow geometries investigated. These are 
now limited by the restriction of gap sizes permissible for an electromagnet 
of high flux density completely enclosing the fluid field. 

There is moreover a further application of the method of electromagnetic 
induction which appears attractive, and is presently being investigated. This 
is the use of the method as a means of measuring the components of vorticity 
in a fluid flow. The basic equation for the potential gradient being, 


grad = (VxB) - Ye 
if one takes the divergence of both sides of the equation one has 
= dw (Vx B) 


V corl B 


(6) 


(11) 
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In a uniform magnetic field therefore 


> 

V g = B-tL (12) 

where ©. is the vorticity vector. An electrode probe in the form of a cross 
enables an experimental measurement of the Laplacian of the potential to be 
obtained which is a direct measure of the vorticity component in the direction 
of the known magnetic field. In this application the disturbing influence of 
induced current patterns is eliminated and with an orientable field the com- 
ponents of the vorticity could be measured separately. Here again the ad- 
vantages of a self-contained probe- magnet assembly are clear, as is the 


desirability of an experimental determination of the vorticity in turbulent 
shear flows. 
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FIG.2 PROBE ASSEMBLY, SHOWING SHIELDED 
ELECTRODES, OUTPUT CONNECTOR AND 
MICROMETER TRAVERSING MECHANISM 
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Journal of the 


HYDRAULICS DIVISION 


100 FREQUENCY CURVES OF NORTH AMERICAN RIVERS 


E. Kuiper,* M. ASCE 
(Proc. Paper 1395) 


This paper presents frequency curves of maximum annual flood flows on 
important rivers in North America. The frequency curves are drawn as 
straight lines. It was found that the inclination of the lines is a function of 
drainage area, climate and soil conditions. 


SYNOPSIS 


The present study was undertaken in the belief that the frequency curves of 
maximum annual flood peaks for one particular river station can be drawn 
with greater accuracy when due attention is paid to other frequency curves on 
the same or comparable streams. 

After having investigated one hundred frequency curves of major streams 
on the North American Continent, it was concluded: 


1. That frequency curves can generally be drawn as straight lines; 

2. That topography and lakes have little bearing on the inclination of a 
frequency line; 

3. That a dry climate and porous soils tend to steepen the frequency line; 

4. That an increase in size of drainage area tends to flatten the frequency 
line. 


It is believed that the results of this study can be of help in drawing fre- 
quency curves of flood flows of rivers that resemble the ones that are ana- 
lysed in this study. 


INTRODUCTION 


In most cases a frequency curve of maximum annual flood peaks at a given 
river station is based on the flood peaks recorded at that particular river 


Note: Discussion open until March 1, 1958. Paper 1395 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, Vol. 
83, No HY 5, October, 1957. 


* Chf. Engr., Manitoba Water Resources Investigation, Manitoba, Canada. 
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station. Little or no attention is paid to frequency curves at river stations 
upstream or downstream, or in similar drainage basins; however, it would 
seem reasonable that some relationship exists between such frequency 
curves. In order to investigate this possibility, a study was made of the fre- 
quency curves of maximum annual flood flows at one hundred river stations 
on the more important rivers in North America. 

In order to facilitate the comparison between the different frequency 
curves, every series of flood peaks was converted into ratio-to-the-mean 
values. For plotting position and graph paper, the Hazen method and the 
Hazen log-probability paper were chosen. After the points were plotted, a 
method of drawing the frequency curves had to be selected. It was desirable 
to adopt the straight-line-by-eye method since this would facilitate the com- 
parison of the different frequency curves. In order to determine if this would 
be justified, special attention was given to the plotted points of river stations 
with the longest period of record, and also some experiments were carried 
out that will be discussed in the following paragraphs. It was concluded that 
for the present study the use of the straight-line method was permissable. 

The main part of the study was then devoted to the questions “if” and ‘how” 
the inclination of the straight frequency line is related to various drainage 
basin characteristics. Plain reasoning suggested that certain characteristics 
would have a certain bearing upon the inclination of the lines. These suppo- 
sitions were then checked with the available data. 


Frequency Curves as Straight Lines 


While studying the one hundred frequency curves that are presented in this 
paper, it appeared that the flood peaks from the longest periods of record had 
a tendency to fall on a straight line, when plotted on logarithmic probability 
paper. In Table 1 are assembled the stations with a period of record of sixty 
years and longer. In the last column is shown to what extent the plotted 
points fall on a straight line. It can be seen from the Table, and the corres- 
ponding Figures, that twelve out of twenty-two series of records can be rep- 
resented by a straight line. Nine out of the remaining ten series approach a 
straight line fairly well. In only one case, at St. Paul on the Mississippi 
River, would a straight line fit poorly, and then, only with respect to the 
lower third of the points. 

Demonstration of the misleading appearance of a number of points on a 
seemingly curved line is provided by the following experiment. A straight 
line was drawn on log-probability paper. From this line, one thousand dis- 
charge values were read at the frequency values: 0.05, 0.15, 0.25, 0.35, etc., 
up to 99.95. The one thousand items that were thus obtained represented a 
true straight frequency line. From this group of one thousand items, fifty 
items were selected, one by one, and by pure chance. This was repeated un- 
til twenty groups of fifty items each were obtained. These twenty series were 
plotted on frequency paper. Five series fitted the original straight line well; 
nine series fitted the line fairly well; six series fitted the line poorly; two of 
them would not fit a straight line at all! 

It is recognized that such an experiment is no conclusive proof that fre- 
quency lines have to be straight. The basic series in this example had a 
straight-line distribution to start with, which is an unknown quantity in an ac- 
tual case. However, it is demonstrated that a curved appearance of the 
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TABLE 1 


Stations with a period of record of 60 years and longer 


Dr. Area Period of Plotted Points 


Fig.* No.** River Station in SqeMie Record on Str. Line 


Cumberland Nashville 12,900 12 
Connecticut Hartford 10,500 107 
Columbia The Dalles 237 ,000 
Arkansas Little Rock 157,900 
Mississippi St. Louis 701,000 
Mississippi Memphis 932,800 
Mississippi Vicksburg 1, 
Tennessee Knoxville 

Tennessee Chattanooga 

Tennessee Florence 

Ohio Pittsburgh 

Ohio Cincinnati 

Ohio Metropolis 

Mississippi St. Paul 

Mississippi Lé Claire 

Mississippi Keokuk 

Savannah Augusta 

Kanawah Kanuchin Falls 

Missouri Fort Benton 

Missouri Kansas City 

Chippewa Chippewa Falls 

Hudson Mechanicville 


* 


Figure refers to Figure number. See Figure 1 for meaning of symbols. 
** Number refers to number in Table 7. 
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plotted points can easily be a chance deviation from a true straight line 
distribution. 

From the observations presented in Table 1, and from the above experi- 
ment it is tentatively concluded for the present studies that frequency curves 
can be drawn as straight lines unless there is conclusive evidence that this 
would not be justified. 

In about twenty-five out of the one hundred cases, mostly containing the 
stations with a short period of record, it was necessary to give some special 
thought to the drawing of a straight line through the plotted points. Difficul- 
ties were caused mostly by irregular scattered points from low to high, or by 
a break in the general alignment of the points, or by one or a few of the high- 
est points that were out of the general direction. In these cases, the follow- 
ing considerations were made before the line was drawn. 

Widely scattered points are usually caused by a short period of record 
that is not representative for the real distribution of the flood peaks. In such 
a case, it is difficult to draw a line in which confidence can be placed. The 
only thing to do is to draw a line through the general direction of the points. 
In the final appraisal of frequency curves, not too much weight should be 
given to the curves that have to be derived in this way. 

A break in the general alignment of the points is either caused by a non- 
representative period of record or by special circumstances. The last pos- 
sibility was investigated first and if no abnormalities were found, the line 
was drawn by giving most of the weight to the upper part of the points, since 
in this study, most of the interest in the frequency curves was focussed on 
the region with the low frequencies. 

One or a few of the highest points that were either too low or much too 
high with respect to the general direction of the other points caused some 
trouble. There were cases where the highest point was so much higher than 
the other points that there could be no question as to this point not belonging 
in the series. There were other cases, however, where the highest point was 
rather high but where it was difficult, if not impossible, to determine whether 
or not this point properly belonged in the records. 

There are some considerations that were kept in mind when dealing with 
these troublesome high points. First of all, the total number of station years 
that was available for the present study amounts to about five thousand. They 
are not all independent items because the stations are sometimes close to- 
gether, or on one river, so that correlation is inevitable. However, it can be 
assumed that they are equivalent to at least, say, a few thousand independent 
streamflow data. In such an extensive series, it can be expected that some of 
the relatively highest peak flows would fit nicely in a, say, 2,000-year period 
of record. However, these high flows are bound to fall out of place in a 50- 
year period of record. 

Furthermore, there is the fact that in the present collection of streamflow 
data, a few river stations have been especially selected because of the exces- 
sive high floods that they contained; for instance, seven stations on the Kansas 
River and tributaries and one on the Neosha River, all affected by the 1951 
floods; four stations on the Arkansas River, affected by the 1921 flood; and the 
Miami River at Dayton, affected by the 1913 flood. These extremes may not 
even belong in a record period of a few thousand years and would have to be 
studied separately. 

Finally, there is the probability which is almost a certainty, that due to 
chance, some periods contain records that are relatively too high and that 
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other periods contain records that are relatively too low. Even if it is as- 
sumed that all the records together represent a period of a few thousand 
years, then it must still be expected that in the breakdown into one hundred 
Small periods, some of these periods receive more and others less than their 
proper share of large floods. 


As a result of the above considerations, the frequency curves in this study 
are drawn by giving more weight to the general direction, indicated by the 
majority of the plotted points, than to the position of one or a few of the high- 
est points. 


The Influence of Drainage Basin Characteristics 


The peak flow of a flood depends upon many variables. First of all, there 
is the factor of chance. In one year, the flood peak may be large; in the next 
year, it may be small. This factor of chance is the principal cause of the 
variability of peak flows that can be found at any stream gauging station. 
However, chance is not the only factor that determines the variation between 
the extreme low and high peaks. The characteristics of the drainage basin 
also play a role. 

In the following paragraphs, the various factors that influence the magni- 
tude of the peak flow will be discussed. The purpose of the discussion is to 
find out to what extent these factors have a bearing on the inclination of the 
frequency curves of the flood peaks. This discussion will be in general terms 
only. The following section will deal with specific cases, and it will be seen 
how the findings of this section can be applied to the frequency curves of one 
hundred river stations on the North American Continent. 


Effect of the Size of the Drainage Area 


Consider two river stations a and b. The respective drainage areas of a 
and b, and also the climatological conditions, are identical. Hence, the fre- 
quency curve of flood flows at a must be the same as the one at b. What type 
of frequency curve can be expected at c, just below the confluence of the two 
rivers? If there was complete correlation in magnitude and complete co- 
incidence in time between the flood peaks at a and b, the peaks at c would be 
the sum of the peaks at a and b. When expressed in ratio to the mean, the 
curve c will coincide with curves a and b. 

However, it is very seldom that complete correlation and complete co- 
incidence occur in nature. For practical purposes, it can be assumed that 
they exist only partly or not at all. In order to determine the effect of cor- 
relation and co-incidence separately, it will be assumed that first of all there 
is a complete correlation and no co-incidence. This means that the peak 
flows in one year are the same at a and b, but that they occur at different 
dates. When the peaks at a and b are far enough apart, the peaks at c will be 
practically the same as the ones at a or at b. When the peaks are fairly close 
together, that is, when there is part-coincidence, the peaks at c will consist 
of the peaks at a (or b), plus a part of the peaks at b (or a). When a long 
enough period of time is considered, the years with some co-incidence will 
cancel the years with almost complete co-incidence and the practical result 
is a certain degree of co-incidence. The curve c can then be obtained by add- 
sing to the peak flows of a, a certain percentage of the peak flows at b. When 
expressed in ratio to the mean flood, curve c again coincides with the curves 
a and b. 
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Secondly, it will be assumed that there is no correlation at all, and com- 
plete co-incidence. This means that the peak flows at a and b occur every 
year on the same date but that there is no relationship between the magni- 
tudes. The one may be extremely large while the other is average or even 
extremely small. This condition results in a curve c, obtained by combining 
items of curve a with arbitrary items of curve b, after representative series 
have been taken from each curve. When expressed in ratio to the mean flood, 
curve c is flatter than curves a and b. 

From the above discussion, it follows that frequency curves have a ten- 
dency to become flatter when the drainage area increases, provided that the 
drainage basin characteristics remain the same. This tendency towards a 
decrease in ratio between the extreme floods and the mean flood is caused by 


the non-existence of complete correlation between the peak flows on the tribu- 


taries and main stem of a river system. One part of a drainage basin may 
produce a large flood but that does not necessarily imply that the other parts 
produce large floods also. The larger the drainage basin, the less correla- 
tion there is between the production of floods in the various sub-basins, and 
the smaller the ratio between the extreme floods and the mean. 


Effect of Topography 


Consider two drainage basins, identical except for the fact that one has a 


topography that is twice as steep as the other. What effect will this have upon 


the frequency curve of the flood flows? Steep slopes mean a quick runoff of 
the precipitation and therefore, a unit-hydrograph with a relatively sharp 
peak. On account of the principle of the unit-hydrograph, it can be expected 
that a change in its shape will affect all floods to the same extent. Hence, all 
flood peaks in the drainage basin with the steeper topography are a certain 
ratio larger than the flood peaks in the other basin, provided that the runoff 
co-efficient remains the same. In other words, topography has no direct 
bearing on the steepness of the frequency curve. It was assumed in the above 
that the runoff co-efficient remained the same in both cases and at all flows. 
However, it can be expected that for low flows, the runoff coefficient is 
smaller than for high flows, and that the difference is smaller with steeper 
topography. This secondary effect of a changing runoff co-efficient will be 
discussed and included under the section “Effect of Soils.” 


Effect of Climate 


Consider two drainage basins, identical except for the fact that one is 
located in a region with a wet climate, like the Atlantic Coast of the United 
States; and the other in a region with a dry climate, like the Canadian 
Prairies. The wet climate will result in a relatively large figure of average 
annual rainfall and it can be expected that the magnitude of the mean flood is 
also relatively large. However, it has been learned from experience that the 
extreme rainstorms, and therefore, the extreme floods, may very well be of 
the same order of magnitude in both places. There are even dry regions 
where the extreme rainfalls are larger than in wet regions. In terms of 
flood magnitudes, this means that the dry regions can be expected to produce 
a relatively large ratio between extreme floods and the mean flood. In terms 


of the frequency curve, it means that the drainage basin in the dry region will 
show a relatively steep curve. 
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Effect of Soils 


Consider two drainage basins, identical except for the fact that one has an 
impervious soil, resulting in a high runoff coefficient, and the other has a 
loose porous soil, resulting in a low runoff coefficient. What effect will this 
have upon the inclination of the frequency curve? In the case of the impervi- 
ous soil, the runoff coefficient may range from a low of, say, sixty per cent 
during moderate rainfalls to a high of ninety per cent during heavy rainfalls 
of short duration. This means that the peak flows of the floods are approxi- 
mately proportional to the size of the rainstorm. 

In the drainage basin with the pervious soils, it can be expected that during 
moderate rainfalls, the infiltration capacity of the ground is so large that 
most or all of the precipitation will be absorbed in the soil and hence, that 
very little water will reach the streams. The runoff coefficient may be as 
low as a few per cent. During heavy rainfalls, however, especially if they are 
of a short duration, the water that infiltrates in the ground constitutes only a 
fraction of the total amount of precipitation and as a result, a relatively large 
proportion will reach the streams. The runoff coefficient may become as 
high as, say, fifty per cent. This means that in the case of the pervious soils, 
the peak flows increase much more rapidly than the size of the rainstorm and 
as a result, the ratio between the extreme floods and the mean flood is larger 
than in the case of the impervious soils. In other words, impervious soils 
tend to produce flatter frequency curves. 


Effect of Natural Storage 


Natural lakes have a tendency to reduce small and large flood peaks in ap- 
proximately the same proportions. Consider, for instance, two river stations 
—a and b, one just upstream and one just downstream of a natural lake. It 
will be interesting to determine the frequency curve at b. For that purpose 
various surface areas of the lake were assumed and different types of floods 
were routed through the lake in order to determine the peaks at b. In each 
case, it was found that the ratio of the original peak to the reduced peak was 
almost independent of the size of the flood. In other words, natural lakes have 
little bearing on the inclination of a frequency curve. This conclusion does 
not hold for situations where the natural storage is carried over from one 
flood to another. In that case small flood peaks may be increased by a carry- 
over from preceding large peaks, with the result that the frequency curve will 
tend to flatten. In the extreme case of infinite storage the frequency line be- 
comes horizontal. 


Analysis of 100 Frequency Curves 


It will be shown in this section that the effect of various drainage basin 
characteristics, as discussed in the previous section, can be found in the fre- 
quency curves of actual river flows. The purpose of this part of the study is 
to find out if relationship between drainage basin characteristics and fre- 
quency curve inclination is so well-defined that it can be used as guidance in 
determining the frequency curve of flood peaks of a particular river station. 

In order to facilitate the analysis, it would be desirable if the inclination of 
the frequency line could be expressed in some numerical parameter. One of 
the simplest ways of doing this would be to quote the point where the frequency 
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line intersects the vertical line, representing the frequency value of one per 
cent. This parameter will be called R, which is then equivalent to the ratio 
of the one per cent flood to the mean flood. The problem is now to determine 
how R is influenced by the size of drainage basin, the climate, and the soils. 
It would be desirable if these independent variables could also be expressed 
in numerical parameters. 

The first variable, the size of the drainage basin, can be expressed easily 
in numerical terms; namely, in square miles. The second and third variables 
are more difficult to express. Fortunately, they both produce similar results 
with regard to R and the mean flood. When the climate is dry, the mean flood 
is bound to be low and R can be expected to be high. When the soils are 
porous, the mean flood will also be low and R also be high. Therefore, it 
would seem logical to relate R to the mean flood instead of taking the climate 
and the soils both into consideration. The mean flood in itself is not a desir- 
able indication because it also depends upon the size of the drainage basin. 

In order to overcome that difficulty the mean flood will be divided by the area 
of the drainage basin to the 0.8 power. Thus, a parameter F: “Flood Coef- 
ficient” is obtained that has from experience proven to be fairly constant 
when the characteristics of the drainage basin are constant, and when only 
the size of the basin varies. It can then be expected that, other variables be- 
ing constant, R will be low when F is high, and vice versa. 

Thus, there are two independent variables: size of drainage are D. A. and 
flood coefficient F; and the one dependent variable: the ratio of the one per 
cent flood to the mean flood R. Since there are two independent variables, it 
can be expected that the effect of one variable may obscure or even dominate 
the effect of the other one. For this reason the analysis will first of all be 
concerned with examples where either F or D. A. is kept more or less con- 
stant, and afterwards with examples where both vary at the same time. 

The first example concerns the Tennessee Basin, where F is fairly con- 
stant and where the effect of D. A. upon R can be observed. It can be seen 
that R decreases when D. A. increases. The data are assembled in Table 2. 

The second example concerns the Red River of the North Basin where F is 
also fairly constant. It can again be seen from Table 3 that R decreases 
when D. A. increases. 

The third example concerns a variety of river stations, all with approxi- 
mately the same D. A., but with different values of F. It can be seen from 
Table 4 that R decreases with an increase of F. 

The fourth example shown in Table 5 concerns the Arkansas Basin, where 
F decreases with increasing drainage area. Although R is supposed to de- 
crease with increasing drainage area, the variable F apparently dominates to 
such an extent that R increases instead of decreases. 

The last example shown in Table 6 concerns the Mississippi Basin, where 
F increases with increasing drainage area. Here, a rapid decrease of R with 
increasing D. A. can be noted because the increasing F and the increasing 
D. A. have the same effect upon R. 

In order to demonstrate the variation of R as a function of F and D. A. for 
all 100 frequency curves, Table 7 and Figure 25 have been prepared. In 
Table 7 all pertinent data of the 100 frequency curves have been listed. Al- 
though it was decided to determine R by drawing a straight frequency line by 
eye, the values of Cy, Cg, and R, (the computed one per cent ratio, by means 
of the Foster III table) have nevertheless been listed, in case they may prove 
to be of interest for future investigations. In Figure 25 the data have been 


| 
| 


ASCE KUIPER 1395-9 


TABLE 2 
| Figs No, River Station F, 
1 62 French Broad Ashville 950 72.6 2.6 
7 63 Tennessee Knoxville 8,900 65.7 2h 
64, Tennessee Chattanooga 21,400 69.0 1.8 
16 65 Tennessee Florence 30,800 64.4 1.7 
TABLE 3 
4 6 Pembina Neche, N. D. 3,200 3.0 5.0 
| 1 R.R. of the N. Fargo, N.D. 6,800 2.9 heb 
2 2 8 " " Grand Forks, N.D. 15,600 4.6 3.7 
| © © Emerson, Man. 36,200 36 
3 5 " Winnipeg, Man. 104,500 363 3.0 
TABLE 4 
5 nu Qu' Appelle Tantallon, Sask. 17,500 Oh 58 
69 Arkansas La Junta, Colo. 12,200 
29 Minnesota Mankato, Minn. 14,900 6.2 3.46 
95 Snake Neeley, Idaho 14,100 10.9 2.8 
39 Des Moines Keosanqua, Iowa 13,900 20.0 245 
15 61 Cumber land Nashville, Tenn 12,900 63.0 1.7 
FABLE 
66 Arkansas Salida 1,200 10.7 2.0 
17 67 bd Canon City 3,100 7.7 2.8 
68 n Pueblo 4,700 Tol 302 
69 La Junta 12,200 
18 70 Holly 25,100 2.69 
Wichita 40,200 1.9 5.0 
TABLE 6 
28 Minnesota Mont evideo 6,200 3el 
29 Mankato 14, 5900 662 346 
7 25 Mississippi St. Paul 36,800 8.3 2.5 
8 26 " Le Claire 88,600 15.1 2.1 
27 " Keokuk 119,000 154 2.0 
20 76 Vicksburg 1,1,,500 181 1.7 
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Period of record used. 

5 highest items in 
descending order. 
Drainage Area, 

Mean of annual peak flow 
Coefficient of Variation 
Coefficient of Skew 


= Ratio of 1% flood to M.F. 
= Greatest historical cen 


ACCORDING TO 
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Figure 2 [No. 2 in Table 7] 


RED RIVER OF THE NORTH 
Grand Forks, N. D. 


Records 1902 - 1951 
Max. 1950, 1947, 1948, 1904, 


Hist. 1897: 80,000 cfs. = 5.1 M.F 
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Figure 3 (No. 5 in Table 7] 


RED RIVER OF THE NORTH 
Winnipeg - Redwood, Man. 


Records 1902-1951 

Max. 1950, 1916, 1948, 1904, 
D. A. 104,500 sq. mi. | 
M. F. 34,100 cfs. 
Cy = 0.54 

F = 3.3 R = 3.0 
Hist. 1826: 225,000 cfs = 6.6 M.E 


| 


FREQUENCY IN PERCENT 
(No. 6 in Table 7] | 


PEMBINA RIVER, 
Neche, N. D. 


Max. 1950, 1949, 1913, 1948, 
D. As 3,200 sae mi. 
M. F. 1,910 cfs. 


Cy s 1.00 


F s 3.0 


Records 1913 - 1951 | 1 
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Figure 5 [No. 11 in Table 7] 


QU'APPELLE RIVER, 
Tantallon, Sask. 


Records 1913-1951 
Max. 1925, 1948, 1927, 1923, 
D. A. = 17,500 sq. mi. 


igure 6 [No. 22 in Table 7] 


SASKATCHEWAN RIVER, 
Nipawin, Sask. 


Records 1911-1951 
Max. 1915, 1932, 1948, 1929, 
D. Aw = 97 ,O00 mi. 
M. F. = 90,000 cfs. 
s 0.50 


= 9.6 
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Figure 7 [No. 25 in Table 7] 


MISSISSIPPI RIVER, 
St. Paul, Minn. 


Records 1878-1950 

Max. 1881, 1897, 1916, 1908, 
D. A. 36,800 sq. mi. 

M. F. 37,100 cfs. 

Cy = 0.51 
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Figure 8 (No. 26 in Table 7] 


MISSISSIPPI RIVER, | 
Le Claire, Iowa. 


| 
Max. 1880, 1888, 1892, 


De. Ae 88,600 SQe mi. | 
M. F. 136,600 cfs. 
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| Figure 9 [No. 28 in Table 7] 
MINNESOTA RIVER, 
Montevideo, Minn. 
| Records 1909-1950 
| 


Max. 1919, 1917, 1943, 1920, 
De As 6,200 mi. 
| M. F. 3,400 cfs 


Figure 10 (No. 32 in Table 7] 


| CHIPPEWA RIVER, 
Chippewa Falls, Wis. 


\ Records 1888-1950 
Max. 1941, 1943, 1920, 1938, 1905 
D. A. 5,600 sq. mi. 
M. F. 39,400 cfs. 
= 1.20 


F 40.0 = 24 


| 


98 95 90 80 70 60 50 40 30 20 10 5 2 ' O05 0201 


| ‘ 
| ° 
1.0 
| 9 
| 98 95 90 80 70 60 50 40. 30 20 10 5 2 | 05 0.2 Ol 50 
| 
: 
= 
— 


Figure 11 [No. 42 in Table 7] 


MISSOURI RIVER, 
Kansas City, Mo. 


Records 1881-1951 

Max. 1951, 1903, 1908, 1943, 
D. A. 489,200 sq. mi. 

M. F. 193,500 cfs. 

Cc s O44 

v 

2 5h R 
Hist. 1844: 625,000 cfs = 3. 


Figure 12 [No. 52 in Table 7] 


KANSAS RIVER, 
Topeka, Kans. 


Records 1905-1951 


Max. 1951, 1935, 1908, 1941, 1943 


D. A. 56,700 sq. mi. 
M. F. 60,700 cfs. 
Cy s 0.64 


F 


s 1.24 
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Figure 13 (No. 58 in Table 7] 


OHIO RIVER, 
Metropolis, Ill. 


Records 1875-1949 
Max. 1937, 1884, = 1883, 1897 


0.79 


Figure 1, (No. 59 in Table 7] 


KANAWHA RIVER, 
Kanawha Falls, W. Va. 


Records 1877-1949 


Max. 1878, 
D. A. 8,400 sq. mi 
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Figure 15 (No. 61 in Table 7] 


CUMBERLAND RIVER, 
Nashville, Tenn. 


Records 1838-1949 

Max. 1927, 1882, 1848, 1937, 
D. A. 12,900 sq. mi. 

M. F. 122,200 cfs. 

Cy = 0.23 


F 63.0 


—. 16 a 65 in Table 7] 


TENNESSEE RIVER, 
Florence, Ala. 


Records 1874-1949 

Max. 1897, 1882, 1875, 1948, 
D. Ae 30,800 mi. 

M. 251,000 cfs. 

s 0.28 


F s 
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Figure 17 [No. 67 in Table 7] 


ARKANSAS RIVER, 
Canon City, Colo. 


Records 1895-1950 
ax. 1921, 1941, 1937, 1944, 1942 
is 
= 
= 2.8 


FREQUENCY PERCENT 


Figure 18 [No. 70 in Table 7] 


ARKANSAS RIVER, 
Holly, Colo. 


Records 1913-1950 

Max. 1921, 1942, 1927, 1925, 1936 
D. Aw 25,100 sq. mi. 

M. F. 10,900 cfs. 

= = 0.90 Cc = 1.40 


F s 2.9 
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Figure 19 [No. 72 in Table 7] 


ARKANSAS RIVER, 
Little Rock, Ark. 


Records 1872-1951 
D. A. 157,900 sq. mi. 
M. F. 291,300 cfs. fe, 
=0.43 = 1.32 

F «20.2 R = 2.3 


02 0.1 


16.0 
Figure 20 [No. 76 in Table 7] 


MISSISSIPPI RIVER, 
Vicksburg, Miss. 


Records 1872-1951 
Max. 1937, 1945, 1950, 1922, 
D. Aw 1,144,500 sq. mi. 

M. F. 1,271,200 cfs. 
F 


s 0.23 
= 18.1 R =#1.7 
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Figure 21 (No. 77 in Table 7] 


NIAGARA RIVER, 
Buffalo, N. Y. 


Records 1905-1950 

Max. 1929, 1947, 1943, 1913, 
D. A. 263,000 sq. mi. 

M. F. 208,600 cfs. 

s 0.09 


= 9.6 


Figure 22 [No. 81B in Table 7] 


CONNECTICUT RIVER, 
Hartford, Conn. 


Records 1843 - 
Max. ime, 1928, 1862 
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Figure 23 [No. 90B in Table 7] 


COLUMBIA RIVER, 
The Dalles, Oreg. 


Records 1858-1949 
Max. 1894, 1876, 1948, 
D. A. 237,000 sq. mi. 
M. 609,900 cfs. 

Cy =0.29 


Figure 24 [No. 100 in Table 7] 


FRASER RIVER, 
Hope, B.C. 


Records 1912-1950 
Max. 1948, 1950, 1921, 1920, 1936 
De Aw 85,600 mi. 
M, F. 307 ,500 efs. 
# 0.22 
v 
F 34.0 
Hist. 1894: 618,000 cfs = 2.0 M.P 


FREQUENCY IN PERCENT 
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plotted. The size of the drainage area is plotted on the horizontal scale; the 
magnitude of F is expressed in the size of the plotted points; and the value of 
R is plotted on the vertical scale. It can be seen that R tends to decrease 
with an increase in drainage area and an increase in flood coefficient. 

However, there is still an appreciable amount of scatter left in the dia- 
gram. There are several possible reasons for this. 


1. It is likely that not all variable factors that affect the inclination of a fre- 
quency line are properly expressed in D.A. and F. 

2. It is likely that a certain amount of scatter is due to non-representative 
records. 

3. A certain amount of error is involved in drawing the straight frequency 
line. 

4. It could be that a straight line is not the proper representation of the true 
frequency curve. 


In view of the considerable amount of work involved, a further analysis | 
has not been pursued. 


Possible Application 


If an engineer is faced with the problem of drawing a frequency curve of 
maximum annual flood flows for a river station, and if he has decided that : 
the frequency curve can be drawn as.a straight line, and if he has access to 
flood frequency data at other places in the drainage basin or in similar drain- 
age basins, he would be well advised to study all these available flood fre- 
quency data. There is a good possibility that he will find definite trends in the 
inclication of the frequency lines under consideration, in which case he could 
use these trends as guidance for his estimates. 

For instance, suppose that a frequency curve has to be drawn for the an- : 
nual peak flows at Saskatoon, Saskatchewan (Figure 25, Number 19) and let it ra 
be assumed that records were only available from 1932-1951, instead of the 
actual available period of 1911-51. It can be seen from Figure 27, where the 
twenty year period is plotted, that without further information the frequency | 
line would probably be drawn as shown by the dashed line, corresponding to 
an R value of 2.5. However, in view of the short period of record, some fur- 
ther guidance in drawing the line would be desirable. In order to provide this 
guidance, all available F, R and D.A. values for stations in the Nelson River 
Drainage Basin have been plotted in Figure 28. Moreover, lines of equal F 
values have been estimated and drawn in the same figure. These lines are 
based not only on the points in Figure 28, but also on the general trend of F 
values in Figure 25. 

From an appraisal of Figure 28 and using the computed F value of 8.4 from 
Figure 27 the R value of Saskatoon is selected at 3.1 and the frequency line in 
Figure 27 is corrected accordingly. From a comparison of Figures 26 and 
27, it can be seen that the corrected frequency line in Figure 27 is now in 
close agreement with the originally established line in Figure 26. It should 
be noted that a discrepancy remains in the value of the mean flood, due to the 
short period of record. 


Figure 26 [No. 19 in Table 7] 


SOUTH SASKATCHEWAN RIVER, 
Saskatoon, Sask. 


Records 1911-1951 
D.A. 50,900 sq. mi. 
M.F. 54,000 cfs. 
s 0.50 


SOUTH SASKATCHEWAN RIVER, 
Saskatoon, Sask. 


Records 1932-1951 

Max. 1932, 1948, 19475 1942, 
D. A. 50,900 sq. mi. 

M. F. 48,400 cfs, 


F = 8.4 R = 3.1 


ASCE KUIPER 1395-29 
40 30 20 10 5 2 | 0.5 0.2 O1 
| | || titi, 
Pg z 
A 
3 
98 95 90 80 70 60 50 40 30 20 10 5 2 1 os 0.2 0.1 
| 
— 
Sea 
209 
TTT 
9a 95 90 B80 70 60 80 40 30 20 10 5 2 1 OS 0.20.1 


NISVG SOVNIVYC NOSTSN pesieid 


October, 1957 


000‘000'! 


|: 
1395-30 HY 5 
@o - 
N 
Cor | 
N 
| 
| 2 
| 
ai. 
= 
a 
o pe 
= 


KUIPER 
CONC LUSIONS 


1. When maximum annual flood flows, expressed in ratio to the mean, are 
plotted on logarithmic probability paper, they tend to fall on a straight line. 


2. The inclination of the frequency line tends to decrease with an increase of 
drainage area, all other drainage basin characteristics remaining the same. 


3. The inclination of the frequency line tends to decrease with an increase of 
the flood coefficient F (equal to the mean of the annual peak flows in cfs, 
divided by the drainage area in square miles, to the 0.8 power), all other 
drainage basin characteristics remaining the same. 


4. The personal judgment that is involved in drawing or extending the fre- 
quency line of one particular river station can be improved by appraising fre- 
quency lines in the same or in similar drainage basins. 
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Journal of the 
HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 


THE HYDRAULIC DESIGN OF STILLING BASINS: 
HYDRAULIC JUMPS ON A HORIZONTAL APRON (BASIN I) 


J. N. Bradley! M. ASCE and A. J. Peterka2 M. ASCE 
(Proc. Paper 1401) 


FOREWORD 


This is one of a group of six papers on the hydraulic design of stilling ba- 
sins and their associated appurtenances. Although original data from hydrau- 
lic models are presented in tabular form, the papers emphasize practical de- 
sign procedures. Basin sizes and dimensions are given in easy-to-use 
dimensionless forms so that dependable stilling basins may be designed with- 
out the need for exceptional judgment or extensive experience. Sample prob- 
lems are included. 

The “Introduction” and “Experimental Equipment” are included in Paper 
1401 but apply to all papers in this group. Paper 1401 is an academic study 
of the hydraulic jump on a flat floor. Useful information is presented con- 
cerning energy losses, applicability of the hydraulic jump formula, length of 
jump, and a new classification of the various types of jumps. 

Paper 1402 covers the design of a short hydraulic jump stilling basin hav- 
ing an end sill and chute blocks. Paper 1403 describes a shorter stilling ba- 
sin which utilizes baffle piers. In both papers information is given to deter- 
mine the critical dimensions of the stilling basin, the tail water range, and 
the water surface profile in the basin. 

Paper 1404 describes a special type of hydraulic jump basin for use when 
the Froude number of the incoming flow is low, and the jump produces waves 
in the downstream channel. An alternate design and two types of wave sup- 
pressors are also developed. 

Paper 1405 describes the design of a stilling basin having a sloping apron. 
The extra tail water depth required and the merits of a sloping apron are 
evaluated, as is the basin length for a range of apron slopes. 

Paper 1406 develops the design of impact basins for discharges up to 340 
second feet and incoming velocities up to 30 feet per second. No tail water is 
required. The bibliography in this paper applies to all papers in this group. 


Note: Discussion open until March 1, 1958. Paper 1401 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 5, October, 1957. 


1. Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly of 
U. S. Bureau of Reclamation, Denver, Colo. 
2. Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 
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ABSTRACT 


A study of the hydraulic jump from a practical viewpoint is presented and 
the results obtained are correlated with those of others. Included in the study 
are the applicability of the hydraulic jump formula, the length of jump over 
the entire practical range, and a new classification for the various forms of 
the hydraulic jump. 


INTRODUCTION 


Although hundreds of stilling basins and energy dissipating devices have 
been designed in conjunction with spillways, outlet works, and canal struc- 
tures, it is still necessary in many cases to make model Studies of individual 
structures to be certain that these will operate as anticipated. The reason for 
these repetitive tests, in many cases, is that a factor of uncertainty exists, 
which in retrospect is related to an incomplete understanding of the overall 
characteristics of energy dissipators. 

Many laboratory studies have been made on individual structures over a 
period of years, by different personnel, for different groups of designers, 
with each structure having different allowable design limitations. Attempts 
to analyze the assembled data resulted in sketchy and, at times, inconsistent 
results having only vague connecting links. Extensive library research into 
the works of others revealed the fact that the necessary correlation factors 
are non-existent. 

To fill the need for up-to-date hydraulic design information on stilling 
basins and energy dissipators, a research program on this general subject 
was begun with a rather academic study of the hydraulic jump, observing all 
phases as it occurs in open channel flow. With a broader understanding of 
this phenomenon it was then possible to proceed to the more practical aspects 
of stilling basin design. 

Existing knowledge, including laboratory and field tests collected from 
Bureau of Reclamation records and experiences over a 23-year period, was 
used to establish a direct approach to the practical problems encountered in 
hydraulic design. Hundreds of tests were also performed on both available 
and specially constructed equipment to obtain a fuller understanding of the 
data at hand and to close the many loopholes. Testing and analysis were 
synchronized to establish valid curves in critical regimes, providing suffi- 
cient understanding of energy dissipators in their many forms to establish 
workable design criteria. Since all the test points were obtained by the same 
personnel, using standards established before testing began, and since results 
and conclusions were evaluated from the same datum, the data presented are 
believed to be consistent and reliable. 

These papers, therefore, are the result of an integrated attempt to genera- 
lize the design of stilling basins, energy dissipators, and associated appur- 
tenances. In most instances design rules and procedures are clearly stated 
in simple terms with limits fixed in a definite range. In other cases, however, 
it is necessary to state procedures and limits in broader terms, making it 
necessary to carefully read the accompanying text. 

Proper use of the material in these papers will eliminate the need for hy- 
draulic model tests on many individual structures, particularly the smaller 
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ones. Structures obtained by following the recommendations in this report 
will be conservative in that they will contain a desirable factor of safety. 

However, to further reduce structure sizes, to account for unsymmetrical 
conditions of approach or getaway, or to evaluate other unusual conditions not 
covered in this discussion, model studies will still prove beneficial. 


Experimental Equipment 


Five test flumes were used at one time or another to obtain the experimen- 
tal data required on the hydraulic jump basins, Flumes A and B, Figure 1; 
Flumes C and D, Figure 2; and Flume F, Figure 3. The arrangement shown 
as Flume E, Figure 3, actually occupied a portion of Flume D during one 
stage of the testing, but it will be designated as a separate flume for ease of 
reference. Each flume served a useful purpose either in verifying similarity 
or extending the range of the experiments. Flumes A, B, C, D, and E con- 
tained overilow sections so that the jet entered the stilling basin at an angle 
to the horizontal. The degree of the angle varied in each case. In Flume F, 
the entering jet was horizontal since it emerged from under a vertical slide 
gate. 

The experiments were started in an existing model of a flat-chute spillway, 
Figure 1A, having a small discharge and low velocity. This was not an ideal 
piece of equipment for general experiments as the training walls on the chute 
were diverging. The rapid expansion caused the distribution of flow entering 
the stilling basin to shift with each change in discharge; nonetheless, this 
piece of equipment served a purpose in that it aided in getting the research 
program underway. 

Tests were then continued in a glass-sided laboratory flume 2 feet wide and 
40 feet long in which an overflow section was installed, Flume B, Figure 1. 
The crest of the overflow section was 5.5 feet above the floor, while the down- 
stream face was on a slope of 0.7:1. The capacity was about 10 cfs. 

Later, the work was carried on at the base of a chute 18 inches wide having 
a Slope of 2 horizontal to 1 vertical and a drop of approximately 10 feet, 
Flume C, Figure 2. The stilling basin had a glass wall on one side. The dis- 
charge capacity was 5 cfs. 

The largest scale experiments were made on a glass-sided laboratory 
flume 4 feet wide and 80 feet long, in which an overfall crest with a slope of 
0.8:1 was installed, Flume D, Figure 2. The drop from headwater to tail water 
in this case was approximately 12 feet, and the maximum discharge was 
28 cfs. 

The downstream end of the above flume was also utilized for testing small 
overflow sections 0.5 to 1.5 feet in height. The maximum discharge used was 
10 cfs. As stated above, this piece of equipment will be designated as Flume 
E, and is shown on Figure 3. 

The sixth piece of equipment was a tilting flume which could be adjusted 
for slopes up to 12°, Flume F, Figure 3. This flume was 1 foot wide by 20 
feet long; the head available was 2.5 feet, and the flow was controlled by a 
Slide gate. The discharge capacity was about 3 cfs. 

Each flume contained a head gage, a tail gage, a scale for measuring the 
length of the jump, a point gage for measuring the average depth of flow enter- 
ing the jump, and a means of regulating the tail water depth. The discharge in 
all cases was measured through the laboratory venturi meters or portable 
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FIGURE 1 


TEST FLUME A 
Width of basin 5 feet, drop 3 feet, discharge 6 cfs 


TEST FLUME B 
Width 2 feet, drop 5.5 feet, discharge 10 cfs 
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FIGURE 3 


TEST FLUME (E) 


Width 4 feet, Drop 0.5-1.5 feet, Discharge 10 cfs 


TEST FLUME (F) 
Adjustable tilting type, maximum slope 12 degrees, 
width one foot, discharge 5cfs 
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venturi-orifice meters. The tail water depth was measured by a point gage 
operating in a stilling well. The tail water depth was regulated by an adjust- 
able weir at the end of each flume. 

Other test setups were used to develop the impact basin, wave suppressor, 
etc. These are described in the appropriate sections. 


General Investigation of the Hydraulic Jump on Horizontal Apron (Basin I) 


Introduction 


A tremendous amount of experimental, as well as theoretical, work has 
been performed in connection with the hydraulic jump on a horizontal apron. 
To mention a few of the experimenters who contributed basic information, 
there are: Bahkmeteff and Matske,! 9 Safranez,? Woycicki,* Chertonosov,!® 
| Einwachter,!1 Elms,!2 Hinds,!4 Forcheimer, Kennison, 22 Kozeny,23 

Rehbock,24 Schoklitsch,29 Woodward,26 and others. There is probably no 
phase of hydraulics that has received more attention, yet, from a practical 
viewpoint, there is still much to be learned. 

The first phase of the present study (Section 1) was academic in nature 
consisting of observing the hydraulic jump in its various forms and correlat- 
ing the results with those of others; the primary purpose being to become bet- 
ter acquainted with the overall jump phenomenon. The objectives were: (1) to 
determine the applicability of the hydraulic jump formula for the entire range 
of conditions experienced in design; (2) to determine the length of the jump 
over the entire practical range and to correlate the findings with results of 
other experimenters where possible; and (3) observe, catalog, and evaluate 
the various forms of the jump. 


Current Experimentation 


To satisfactorily observe the hydraulic jump throughout its entire range 
required tests in all six test flumes. This involved about 125 tests, Table 1, 
for discharges of 1 to 28 cfs. The number of flumes used enhanced the value 
of the results in that it was possible to observe the degree of similitude ob- 
tained for the various sizes of jumps. Greatest reliance was naturally placed 
on the results from the larger flumes, since the action in small jumps is too 
rapid for the eye to follow. Friction and viscosity also become a measurable 
factor. This was demonstrated by the fact that the length of jump obtained 
from the two smaller flumes, A and F, was consistently shorter than that ob- 
served for the larger flumes. Out-of-scale frictional resistance on the floor 
and side walls produced a short jump. As testing advanced and this deficiency 


became better understood, some allowance was made for this effect in the 
observations. 


Experimental Results 


Definitions of the symbols used in connection with the hydraulic jump on a 
horizontal floor are shown on Figure 4. The procedure followed in each test 
of this series was to first establish a flow. The tail water depth was then 
gradually increased until the front of the jump moved upstream to Section 1, 
indicated on Figure 4. The tail water depth was then measured, the length of 

| | the jump recorded, and the depth of flow entering the jump, Dj, was obtained 
: by averaging a generous number of point gage measurements taken 
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immediately upstream from Section 1. The results of the measurements and 
succeeding computations are tabulated in Table 1. The measured quantities 
are tabulated as follows: total discharge, (Column 3); tail water depth, 
(Column 6); length of jump, (Column 11), and depth of flow entering jump, 
(Column 8). 

Column 1 indicates the test flumes in which the experiments were per- 
formed, and Column 4 shows the width of each flume. All computations are 
based on discharge per foot width of flume; unit discharges (q) are shown in 
Column 5. 

The velocity entering the jump V1, Column 7, was computed by dividing q 
(Column 5) by Dy (Column 8). 


The Froude Number 


The Froude number, Column 10, Table 1, is: 


(1) 
"Vai 
where Fy is a dimensionless parameter, V1 and Dj are velocity and depth of 
flow, respectively, entering the jump, and g is the acceleration of gravity. 
The law of similitude states that where gravitational forces predominate, as 
they do in open channel phenomena, the Froude number should have the same 
value in model and prototype. Therefore, a model-shaped jump in a test flume 
will have the identical characteristics of a prototype jump in a stilling basin, 
if the Froude numbers of the incoming flows are the same. Although energy 
conversions in a hydraulic jump bear some relation to the Reynolds number, 
gravity forces predominate, and the Froude number becomes most useful in 
plotting stilling basin characteristics. Bakhmeteff and Matzke! demonstrated 
this application in 1936 when they related stilling basin characteristics to the 
2 
square of the Froude number, 2D,’ which they termed the kinetic flow factor. 


The Froude number (1) will be used throughout this presentation. As the 
acceleration of gravity is a constant, the term g could be omitted. Its inclu- 


sion makes the expression dimensionless, however, and the form shown as (1) 
is preferred. 


Applicability of Hydraulic Jump Formula 


The theory of the hydraulic jump in horizontal channels has been treated 
thoroughly by others (see Bibliography), and will not be repeated here. The 


expression for the hydraulic jump, based on pressure-momentum may be 
written; 15 
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Transposing Dj to the left side of the equation and substituting F,? for 
1 


D, 


2 

D 

1/2 + - 1) (3) 


Expression (3) shows that the ratio of conjugate depths is strictly a func- 
D 
tion of the Froude number. The ratio = is plotted with respect to the Froude 
1 
number on Figure 5. The line, which is virtually straight except for the lower 
end, represents the above expression for the hydraulic jump; while the points, 
which are experimental, are from Columns 9 and 10, Table 1. The agreement 
is quite good for the entire range. 

There is an unsuspected characteristic in the curve, however, which is 
mentioned here but will be enlarged on later. Although the tail water depth, 
recorded in Column 6 of Table 1, was sufficient to bring the front of the jump 
to Section 1 (Figure 4) in each test, the ability of the jump to remain at Sec- 
tion 1 for a slight lowering of tail water depth became more difficult for the 
higher and lower values of the Froude number. The jump was least sensitive 


to variation in tail water depth in the middle range, or values of Fj from 4.5 
to 9. 


Length of Jump 


The length of the jump measurement, Column 11, Table 1, was the most 
difficult to determine. Special care was therefore given to this measurement. 
Where chutes or overfalls were used, the front of the jump was held at the in- 
tersection of the chute and the horizontal floor, as shown on Figure 4. The 
length of jump was measured from this point to a point downstream where 
either the high-velocity jet began to leave the floor or to a point on the surface 
immediately downstream from the roller, whichever was the longer. In the 
case of Flume F, where the flow discharged from a gate onto a horizontal 
floor, the front of the jump was maintained just downstream from the com- 
pleted contraction of the entering jet. In both cases the point at which the 
high-velocity jet begins to rise from the floor is not fixed, but tends to shift 
upstream and downstream. This is also true of the roller on the surface. It 
was at first difficult to repeat length observations within 5 percent by either 
criterion, but with practice satisfactory measurements became possible. It 
was the intention to judge the length of the jump from a practical standpoint; 
in other words, the end of the jump, as chosen, would represent the end of the 
concrete floor and side walls of a conventional stilling basin. 

A system devised to measure velocities on the bottom, to aid in determin- 
ing the length of jump, proved inadequate and too laborious to allow comple- 
tion of the program planned. Visual observations, therefore, proved to be the 


most satisfactory as well as the most rapid method for determining the length 
measurement. 
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The length of jump has been plotted in two ways. The first is perhaps the 
better method while the second is the more common and useful. The first 
method is shown in Figure 6 where the ratio, length of jump to Dy (Column 13, 
Table 1), is plotted with respect to the Froude number (Column 10) for results 
from the six test flumes. The resulting curve is of fairly uniform curvature, 
which is the principal advantage of these coordinates. The second method of 
plotting, where the ratio, length of jump to the conjugate tail water depth D2 
(Column 12) is plotted with respect to the Froude number, is presented on 
Figure 7. This latter method of plotting will be used throughout the study. 
The points represent the experimental values. 

In addition to the curve established by the test points, curves representing 
the results of three other experimenters are shown on Figure 7. The best 
known and most widely accepted curve for length of jump is that of 
Bakhmeteff and Matzke! which was determined from experiments made at 
Columbia University. The greater portion of this curve, labeled 1, is at 
variance with the present experimental results. Because of the wide use this 
curve has experienced, a rather complete explanation is presented regarding 
the disagreement. 

The experiments of Bakhmeteff and Matzke were performed in a flume 6 
inches wide, with limited head. The depth of flow entering the jump was ad- 
justed by a vertical slide gate. The maximum discharge was approximately 
0.7 cfs, and the thickness of the jet entering the jump, Dj, was 0.25 foot fora 
Froude number of 1.94. The results up to a Froude number of 2.5 are in 
agreement with the present experiments. To increase the Froude number, it 
was necessary for Bakhmeteff and Matzke to decrease the gate opening. The 
extreme case involved a discharge of 0.14 cfs and a value of Dj of 0.032 foot, 
for F; = 8.9, which is much smaller than any discharge or value of D1 used in 
the present experiments. Thus, it is reasoned that as the gate opening de- 
creased, in the 6-inch-wide flume, frictional resistance in the channel down- 
stream increased out of proportion to that which would have occurred in a 
larger flume or a prototype structure. Thus, the jump formed in a shorter 
length than it should. In laboratory language, this is known as “scale effect,” 
and is construed to mean that prototype action is not faithfully reproduced. 

It is quite certain that this was the case for the major portion of curve 1. In 
fact, Bahkmeteff and Matzke were somewhat dubious concerning the small 
scale experiments. 

To confirm the above conclusion, it was found that results from Flume F, 
which was 1 foot wide, became erratic when the value of Dj approached 0.10. 
Figures 6 and 7 show three points obtained with a value of Dj of approximate- 
ly 0.085. The three points are given the symbol and fall short of the 
recommended curve. 

The two remaining curves, labeled 3 and 4, on Figure 7, portray the same 
trend as the curve obtained from the current experiments. The criterion used 
by each experimenter for judging the length of the jump is undoubtedly respon- 
sible for the displacement. The curve labeled 3 was obtained at the Technical 
University of Berlin on a flume 1/2 meter wide by 10 meters long. The curve 
labeled 4 was determined from experiments performed at the Federal Insti- 
tute of Technology, Zurich, Switzerland, on a flume 0.6 of a meter wide and 7 
meters long. The curve numbers are the same as the reference numbers in 
the Bibliography which refer to the work. 

As can be observed from Figure 7, the test results from Flumes B, C, D, 
E, and F plot sufficiently well to establish a single curve. The five points 
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from Flume A, denoted by squares, appear somewhat erratic and plot to the 
right of the general curve. Henceforth, reference to Figure 7 will concern 
only the recommended curve which is considered applicable for general use. 


Energy Absorption in Jump 


With the experimental information available the energy absorbed in the 
jump may be computed. Columns 14 through 18, Table 1, list the computa- 
tions, and the symbols may be defined by consulting the specific energy dia- 
gram in Figure 4. Column 14 lists the total energy, E1, entering the jump at 
Section 1 for each test. This is simply the depth of flow, Dj, plus the velocity 
head computed at the point of measurement. The energy leaving the jump, 
which is the depth of flow plus the velocity head at Section 2, is tabulated in 
Column 15. The differences in the values of Columns 14 and 15 constitute the 
loss of energy, in feet of water, attributed to the conversion, Column 16. 
Column 18 lists the percentage of energy lost in the jump, EL, to the total 
energy entering the jump, Ey. This percentage is plotted with respect to the 
Froude number and is shown as the curve to the left on Figure 8. Fora 
Froude number of 2.0, which would correspond to a relatively thick jet enter- 
ing the jump at low velocity, the curve shows the energy absorbed in the jump 
to be about 7 percent of the total energy entering. Considering the other ex- 
treme, for a Froude number of 19, which would be produced by a relatively 
thin jet entering the jump at very high velocity, the absorption by the jump 
would amount to 85 percent of the energy entering. Thus, the hydraulic jump 
can perform over a wide range of conditions. There are poor jumps and good 
jumps, with the most satisfactory occurring over the center portion of the 
curve. 

Another method of expressing the energy absorption in a jump is to express 
the loss, Ey, in terms of Dj. The curve to the right on Figure 8 shows the 


ratio = (Column 17, Table 1) plotted against the Froude number. Losses in 
1 


feet of head are obtained from this method. 


Forms of the Hydraulic Jump 


The hydraulic jump may occur in at least four different distinct forms ona 
horizontal apron, as shown on Figure 9. All of these forms are encountered in 
practice. The internal characteristics of the jump and the energy absorption 
in the jump vary with each form. Fortunately these forms, some of which are 
desirable and some undesirable, can be cataloged conveniently with respect to 
the Froude number, as shown in Figure 9. 

When the Froude number is unity, the water is flowing at critical depth; 
thus a jump cannot form. This corresponds to Point 0 on the specific energy 
diagram of Figure 4. For values of the Froude number between 1.0 and 1.7 
there is only a slight difference in the conjugate depths Dj and Dg. A slight 
ruffle on the water surface is the only apparent feature that differentiates this 
flow from flow at critical depth. As the Froude number approaches 1.7, a 
series of small rollers develop on the surface as indicated in Figure 9, and 
this action remains much the same but with further intensification up to a 
value of about 2.5. In this range there is no particular stilling basin problem 
involved; the water surface is quite smooth, the velocity throughout is fairly 
uniform, and the energy loss is low. 
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Figure 9 indicates the type of jump that may be encountered at values of 
the Froude number from 2.5 to 4.5. This type has a pulsating action and is 
usually seen in low head structures. The entering jet oscillates from bottom 
to surface with no regular period. Turbulence occurs near the bottom at one 
instant and entirely on the surface the next. Each oscillation produces a 
large wave of irregular period which in prototype structures has been ob- 
served to travel for miles causing damage to earth banks and riprap. This 
case is of sufficient importance that a separate section, Section 4, has been 
devoted to the practical aspects of design. 

A well stabilized jump can be expected for the range of Froude numbers 
between 4.5 and 9, Figure 9. In this range, the downstream extremity of the 
surface roller and the point at which the high-velocity jet tends to leave the 
floor occur in practically the same vertical plane. The jump is well balanced 
and the action is thus at its best. The energy absorption in the jump for 
Froude numbers from 4.5 to 9 ranges from 45 to 70 percent, Figure 8. 

As the Froude number increases above 9, the form of the jump gradually 
changes to that shown in Figure 9. This is the case where Vj is very high, 
Dj is comparatively small, and the difference in conjugate depths is large. 
The high-velocity jet no longer carries through for the full length of the jump. 
In other words, the downstream extremity of the surface roller now becomes 
the determining factor in judging the length of the jump. Slugs of water rolling 
down the front face of the jump intermittently fall into the high-velocity jet 
generating additional waves downstream, and a rough surface can prevail. 
Figure 8 shows that the energy dissipation for this case is high and may reach 
85 percent. 

The limits of the Froude number given above for the various forms of jump 
are not definite values, but overlap somewhat depending on local factors. Re- 
turning to Figure 7, it is found that the length curve catalogs the various 
phases of the jump quite well. The flat portion of the curve indicates the 
range of best operation. The steep portion of the curve to the left definitely 
indicates an internal change in the form of the jump. In fact, two changes are 
manifest, the form shown in Figure 9A and the form, which might better be 
called a transition stage, shown in Figure 9B. The right end of the curve on 
Figure 7 also indicates a change in form, but to less extent. 


Practical Consideration 


Although it was the intention to stress the academic rather than the prac- 
tical viewpoint in this paper, this is the logical place to point out a few of the 
practical aspects of stilling basin design using horizontal aprons. Consider- 
ing the four forms of jump just discussed, the following are pertinent: 


1. All jump forms are encountered in stilling basin design. 

2. The form in Figure 9A requires no baffles or special devices in the 
basin. The only requirement is to provide the proper length of pool, which is 
relatively short. This can be obtained from Figure 7. 

3. The form in Figure 9B is one of the most difficult to handle and is fre- 
quently encountered in the design of canal structures, diversion or low dam 
spillways, and even outlet works. Baffle piers or appurtenances in the basin 
are of little value. Waves are the main source of difficulty and methods for 
coping with them are discussed in Paper 1404. The present information may 
prove valuable in that it will help to restrict the use of jumps in the 2.5 to 4.5 
Froude number range. In many cases the critical range cannot be avoided, 
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but in others the jump may be brought into the desirable range by altering the 
dimensions of the structure. 

4. No particular difficulty is encountered in the form valuable as a means 
of shortening the length of basin. This is discussed in Papers 1402 and 1403. 

5. As the Froude number increases, the jump becomes more sensitive to 
tail water depth. For numbers as low as 8, a tail water depth greater than the 
conjugate depth is advisable to be certain that the jump will stay on the apron. 
This phase will be discussed in more detail in following Papers. 

6. When the Froude number is greater than 10 the difference in conjugate 
depths is great, and, generally speaking, a very deep basin with high training 
walls is required. On high spillways the cost of the stilling basin may not be 
commensurate with the results obtained. A bucket type dissipator may give 
comparable results at less cost. On lower head structures the action in the 
basin will be rugged in appearance with surface disturbances being of great- 
est concern. 

7. High Froude numbers will always occur for flow through extremely 
small gate openings on even the smallest structures. Unless the discharge 
for these conditions represents an appreciable percentage of the design flow, 
the high Froude numbers have no significance. 


Water-surface Profiles and Pressures 


Water-surface profiles for the jump on a horizontal floor were not mea- 
sured as these have already been determined by Bakhmeteff and Matzke,1 
Newman and LaBoon,!9 and Moore. 27,18 It has been shown by several ex- 
perimenters that the vertical pressures on the floor of the stilling basin are 
virtually the same as the static head indicated by the water-surface profile. 


Conclusions 


The foregoing experiments and discussion serve to associate the Froude 
number with the hydraulic jump and stilling basin design. The ratio of con- 
jugate depths, the length of jump, the type of jump to be expected, and the 
losses involved have all been related to this number. The principal advantage 
of this form of presentation is that one may see the entire picture at a glance 
with a minimum of calculation. 


Application of Results (Example 1) 


Water flowing under a Sluice gate discharges into a rectangular stiiling 
basin the same width as the gate. The average velocity and the depth of flow 
after contraction of the jet is complete are: Vy = 85 ft/sec and Dy = 5.6 feet. 
Determine the conjugate tail water depth, the length of basin required to con- 
fine the jump, the effectiveness of the basin to dissipate energy, and the type 
of jump to be expected. 


Jed, x 5.6 


Entering Figure 5 with this value 
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The conjugate tail water depth 


Do = 8.5 x 5.6 = 47.6 feet 


Entering the recommended curve on Figure 7 with a Froude number of 
6.34 


= = 6.13 
2 


Length of basin necessary to confine the jump 


L = 6.13 x 47.6 = 292 feet 


Entering Figure 8 with the above value of the Froude number, it is found 
that the energy absorbed in the jump is 58 percent of the energy entering. 

By consulting Figure 9 it is apparent that a very satisfactory jump can be 
expected. 

The following Papers deal with the more practical aspects of stilling 


basin design, such as modifying the jump by baffles and sills to increase sta- 
bility and shorten the length. 
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FOREWORD 


This is one of a group of six papers on the hydraulic design of stilling 
basins and their associated appurtenances. Although original data from hy- 
draulic models are presented in tabular form, the papers emphasize practical 
design procedures. Basin sizes and dimensions are given in easy-to-use 
dimensionless forms so that dependable stilling basins may be designed with- 
out the need for exceptional judgment or extensive experience. Sample prob- 
lems are included. 

The “Introduction” and “Experimental Equipment” are included in Paper 
1401 but apply to all papers in this group. Paper 1401 is an academic study of 
the hydraulic jump on a flat floor. Useful information is presented concerning 
energy losses, applicability of the hydraulic jump formula, length of jump, and 
a new classification of the various types of jumps. 

Paper 1402 covers the design of a short hydraulic jump stilling basin hav- 
ing an end sill and chute blocks. Paper 1403 describes a shorter stilling 
basin which utilizes baffle piers. In both papers information is given to de- 
termine the critical dimensions of the stilling basin, the tail water range, and 
the water surface profile in the basin. 

Paper 1404 describes a special type of hydraulic jump basin for use when 
the Froude number of the incoming flow is low, and the jump produces waves 
in the downstream channel. An alternate design and two types of wave sup- 
pressors are also developed. 

Paper 1405 describes the design of a stilling basin having a sloping apron. 
The extra tail water depth required and the merits of a sloping apron are 
evaluated, as is the basin length for a range of apron slopes. 

Paper 1406 develops the design of impact basins for discharges up to 340 
second feet and incoming velocities up to 30 feet per second. No tail water is 
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required. The bibliography in this paper applies to all papers in this group. 


ABSTRACT 


An analysis of the critical dimensions of 36 existing hydraulic jump still- 
ing basins and hydraulic model verification tests are used to generalize the 
design of a stilling basin, utilizing chute blocks and an end sill. Tail water 


ranges, jump lengths, water surface profiles, etc., values are provided in 
dimensionless forms. 


INTRODUCTION 


Stilling basins are seldom designed to confine the entire length of the 
hydraulic jump on the paved apron as was assumed in Paper 1401; first, for 
economic reasons and secondly, because there are means for modifying the 
jump characteristics to obtain comparable or better performance in shorter 
lengths. It is possible to reduce the jump length by the installation of acces- 
sories such as blocks, baffles, and sills in the stilling basin. In addition to 
shortening the jump, the accessories exert a stabilizing effect and in some 
cases increase the factor of safety. 

This study concerns stilling basins of the type which have been used on 
high dam and earth dam spillways, and large canal structures, and will be 
denoted as Basin II, Figure 10. The basin contains chute blocks at the up- 
stream end and a dentated sill near the downstream end. No baffle piers are 
used in Basin II because of the relatively high velocities entering the jump. 
The principal aim in these tests was to (1) generalize the design, and (2) de- 
termine the range of operating conditions for which this basin is best suited. 
Since many of these basins have been designed, constructed, and operated, 
some of which were checked with models, the principal task in accomplishing 
the first objective was to tabulate and analyze the dimensions of existing 


structures. Only structures on which firsthand information was available 
were used. 


Results of Compilation 


With the aid of Figure 10, most of the symbols used in Table 2 are self- 
explanatory. The use of baffle piers is limited to Basin III. Column 1 lists 
the reference material used in compiling the table. Column 2 lists the maxi- 
mum reservoir elevation, Column 3 the maximum tail water elevation, 
Column 5 the elevation of the stilling basin floor, and Column 6 the maximum 
discharge for each spillway. Column 4 indicates the height of the structure 
studied, showing a maximum fall from headwater to tail water of 179 feet, a 
minimum of 14 feet, and an average of 85 feet. Column 7 shows that the width 
of the stilling basins varied from 1,197.5 to 20 feet. The discharge per foot 
of basin width, Column 8, varied from 760 to 52 cfs, with 265 as an average. 
The computed velocity, Vj (hydraulic losses estimated in some cases), enter- 
ing the stilling basin (Column 9) varied from 108 to 38 feet per second, and 
the depth of flow, Dj, entering the basin (Column 10) varied from 8.80 to 0.60 
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feet. The value of the Froude number (Column 11) varied from 22.00 to 4.31. 
Column 12 shows the actual depth of tail water above the stilling basin floor, 
which varied from 60 to 12 feet, while Column 14 lists the computed, or con- 
jugate, tail water depth for each stilling basin. The conjugate depths, Do, 
were obtained from Figure 5, (paper 1401). The ratio of the actual tail water 
depth to the conjugate depth is listed for each basin in Column 15. 


Tail Water Depth 

The ratio of actual tail water depth to conjugate depth shows a maximum of 
1.67, a minimum of 0.73, and an average of 0.99. In other words, on the aver- 
age, the basin floor was set to provide a tail water depth equal to the conju- 
gate or necessary depth. 


Chute Blocks 

The chute blocks used at the entrance to the stilling basin varied in size 
and spacing. Some basins contained nothing at this point, others a solid step, 
but in the majority of cases a serrated device, known as chute blocks, was 
utilized. The chute blocks at the upstream end of the basin tend to corrugate 
the jet, lifting a portion of it from the floor, resulting in a shorter length of 
jump than would be possible without them. These blocks also tend to improve 
the action in the jump. The proportioning of chute blocks has been the subject 
of much discussion. The tabulation in Columns 19 through 24 of Table 2 
shows the sizes which have been used. Column 20 shows the height of the 
chute blocks, while Column 21 gives the ratio of height of block to the depth, 
Dy. The ratios of height of block to Dj indicate a maximum of 2.72, a mini- 
mum of 0.81, and an average of 1.35. This is somewhat higher than was shown 
to be necessary by the verification tests discussed later; a block equal to Dj 
in height is sufficient. 

The width of the blocks is shown in Column 22. Column 23 gives the ratio 
of width of the block to height, with a maximum of 1.67, a minimum of 0.44, 
and an average of 0.97. The ratio of width of block to spacing, tabulated in 
Column 24, shows a maximum of 1.91, a minimum of 0.95, and an average of 
1.15. The three ratios indicate that the proportion: height equals width, 
equals spacing, equals Dj should be a satisfactory standard for chute block 
design. The wide variation shows that these dimensions are not critical. 


Dentated Sill 

The sill in or at the end of the basin was either solid or had some form of 
dentated arrangement, as designated in Column 25. A dentated sill located at 
the end of the apron is recommended. The shape of the dentates and the angle 
of the sills varied considerably in the spillways tested, Columns 26 through 
31. The position of the dentated sill also varied and this is indicated by the 


ratio o in Column 26. The distance, X, is measured to the downstream edge 
II 
of the sill, as illustrated in Figure 10. The ratio es varied from 1 to 0.65, 


with an average of 0.97. 

The heights of the dentates are given in Column 27. The ratio of height of 
block to the conjugate tail water depth is shown in Column 28. These ratios 
show a maximum of 0.37, a minimum of 0.08, and an average of 0.20. The 
width to height ratio, Column 30, shows a maximum of 1.25, a minimum of 
0.33, and average of 0.76. The ratio of width of block to spacing, Column 31, 
shows a maximum of 1.91, a minimum of 1.0, and an average of 1.13. For the 
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purpose of generalization, the following proportions are recommended: 

(1) height of dentated sill = 0.2D9, (2) width of blocks = 0.15Dg9, and (3) spacing 

of blocks = 0.15D9, where Dg is the conjugate tail water depth. It is recom- 

mended that the dentated sill be placed at the downstream end of the apron. 
Columns 32 through 38 show the proportions of additional baffle blocks 

used on three of the stilling basins. These are not necessary and are not 

recommended for this type of basin. 


Additional Details 

Column 18 indicates the angle, with the horizontal, at which the high- 
velocity jet enters the stilling basin for each of the spillways. The maximum 
angle was 34° and the minimum 14°. The effect of the vertical angle of the 
chute on the action of the hydraulic jump could not be evaluated from the in- 
formation available. This factor will be considered, however, in Section 5 in 
connection with sloping apron design. 

Column 39 designates the cross section of the basin. In all but three cases 
the basins were rectangular. The three cross sections that were trapezoidal 
had side slopes varying from 1/4:1 to 1/2:1. The generalized designs present- 
ed in this report are for stilling basins with rectangular cross sections. 
Where trapezoidal basins are contemplated a model study is strongly 
recommended. 

Column 40, Table 2, indicates that in the majority of basins constructed for 
earth dam spillways the wing walls were normal to the training walls. Five 
basins were constructed without wing walls using a rock blanket for protec- 
tion. The remainder utilized angling wing walls or warped transitions down- 
stream from the basin. The latter are common on canal structures. The ob- 
ject, of course, is to build the cheapest wing wall that will afford the 
necessary protection. The type of wing wall is usually dictated by local con- 
ditions such as width of the channel downstream, depth to foundation rock, 


degree of protection needed, etc., thus wing walls are not amenable to 
generalization. 


Verification Tests 


An inspection of the data shows that the structures listed in Table 2 do not 
cover the desired range of operating conditions. There is insufficient infor- 
mation to determine the length of basin for the larger values of the Froude 
number, there is little or no information on the tail water depth at which 
sweepout occurs, and the information available is of little value for generaliz- 
ing the problem of water-surface profiles. Laboratory tests were therefore 
performed to extend the range and to supply the missing data. The experi- 
ments were made on 17 Type II basins, proportioned according to the above 
rules, and installed in Flumes B, C, D, and E (see Columns 1 and 2, Table 3). 
Each basin was judged at the discharge for which it was designed, the length 
was adjusted to the minimum that would produce satisfactory operation, and 
the absolute minimum tail water depth for acceptable operation was measured. 
The basin operation was also observed for flows less than the designed dis- 
charge and found to be satisfactory in each case. 

Table 3 is quite similar to Table 2 with the exception that the length of 
basin Ly (Column 11) was determined by experiment, and the tail water depth 


at which the jump just began to sweep out of the basin was recorded (Column 
13). 
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Tail Water Depth 


The solid line on Figure 11 was obtained from the hydraulic jump formula 
D 


= = 1/2 (V 1 + BF? - 1) and represents conjugate tail water depth. It is the 
1 
Same as the line shown on Figure 5, (paper 1401). The dashed lines on 
Figure 11 are merely guides drawn for tail water depths other than conjugate 
depth. The points shown as dots were obtained from Column 13 of Table 2 
and constitute the ratio of actual tail water depth to Dj for each basin listed. 
It can be observed that the majority of the basins were designed for conjugate 
tail water depth or less. The minimum tail water depth for Basin II, obtained 
from Column 14 of Table 3, is shown on Figure 11. The curve labeled “Mini- 
mum TW Depth Basin II” indicates the point at which the front of the jump 
moves away from the chute blocks. In other words, any additional lowering of 
the tail water would cause the jump to leave the basin. Consulting Figure 11 
it can be observed that the margin of safety for a Froude number of 2 is 0 
percent; while for a number of 6 it increases to 6 percent, for a number of 10 
it diminishes to 4 percent, and for a number of 16 it is 2.5 percent. To be 
certain that this is understood, it will be stated another way. The jump will 
no longer operate properly when the tail water depth approaches 0.98D9 for a 
Froude number of 2 or 0.94D2 for a number of 6 or 0.96D2 for a number of 
10, or 0.975D2 for a number of 16. The margin of safety is largest in the 
middle range. For the two extremes of the curve it is advisable to provide a 
tail water greater than conjugate depth to be safe. For these reasons the 
Type II basin should never be designed for less than conjugate depth and a 
minimum safety factor of 5 percent of Dg is recommended. 

Several precautions should be taken when determining tail water elevations. 
First, tail water curves are usually extrapolated for the discharges encoun- 
tered in design, so they can be in error. Secondly, the actual tail water depth 
usually lags, in a temporal sense, that of the tail water curve for rising flow 
and leads the curve for a falling discharge. Extra tail water should therefore 
be provided if reasonable increasing increments of discharge limit the per- 
formance of the structure because of a lag in building up tail water depth. 
Thirdly, a tail water curve may be such that the most adverse condition oc- 
curs at less than the maximum designed discharge; and fourthly, temporary 


or permanent retrogression of the riverbed downstream may be a factor need- 


ing consideration. These factors, some of which are difficult to evaluate, are 
all important in stilling basin design, and suggest that an adequate factor of 
safety is essential. It is advisable to construct a jump height curve, superim- 
posed on the tail water curve, for each basin to determine the most adverse 
operating condition. This procedure will be illustrated later. 

The verification tests repeatedly demonstrated that there is no simple 
remedy for a deficiency in tail water depth. Increasing the length of basin, 
which is the remedy often attempted in the field, will not compensate for 
deficiency in tail water depth. For these reasons, care should be taken to 
consider all factors that may affect the tail water at a future date. 


Length of Basin 


The necessary length of Basin II, determined by the verification tests, is 
shown as the intermediate curve on Figure 12. The squares indicate the test 
points (Columns 10 and 12 of Table 3). The black dots represent existing 
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basins (Columns 11 and 17, Table 2). Conjugate depth was used in the ordi- 
nate ratio rather than actual tail water depth since it could be computed for 
each case. 

The dots scatter considerably but an average curve drawn through these 
points would be lower than the Basin II curve. In Figure 12, therefore, it ap- 
pears that in practice a basin about 3 times the conjugate depth is actually 
used when a basin about 4 times the conjugate is recommended from the veri- 
fication tests. It should be remembered, however, that the shorter basins 
were all model tested and every opportunity was taken to reduce the basin 
length. The extent and depth of bed erosion, wave heights, favorable flood 
frequencies, flood duration and other factors were all used to justify reducing 
the basin length. Lacking definite knowledge of this type in designing a basin 
for field construction without model tests, the longer basins indicated by the 
verification tests curve are recommended. 

The Type II basin curve has been arbitrarily terminated at Froude num- 
ber 4, as the jump may be unstable at lower numbers. The chute blocks have 
a tendency to stabilize the jump and reduce the 4.5 limit discussed for Basin I. 
For basins having Froude numbers below 4.5 see Paper 1404. 


Water-surface Profiles 


Water-surface profiles were measured during the tests to aid in computing 
uplift pressures under the basin apron. As the water surface in the stilling 
basin tests fluctuated rapidly it was felt that a high degree of accuracy in mea- 
surement was not necessary. This was found to be true when the approximate 
water-surface profiles obtained were plotted, then generalized. It was found 
that the profile in the basin could be closely approximated by a straight line 
making an angle o@ with the horizontal. This line can also be considered to be 
a pressure profile. 

The angle a@ (Column 24, Table 3) observed in each of the verification tests 
has been plotted with respect to the Froude number on Figure 13. The slope 
increases with the Froude number. To use the curve in Figure 13, a horizon- 
tal line is drawn at conjugate depth on a scale drawing of the basin. A verti- 
cal line is also drawn from the upstream face of the dentated sill. Beginning 
at the point of intersection, a sloping line is constructed as shown. The above 
procedure gives the approximate water surface and pressure profile for con- 
jugate tail water depth. Should the tail water depth be greater than Dg, the 
profile will resemble the uppermost line on Figure 13; the angle remains un- 
changed. This information applies only for the Type II basin, constructed as 
recommended in this section. 


CONCLUSIONS 


The following rules are recommended for generalization of Basin II, 
Figure 14: 


1. Set apron elevation to utilize full conjugate tail water depth, plus an 
added factor of safety if needed. An additional factor of safety is advisable 
for both low and high values of the Froude number (see Figure 11). A mini- 
mum margin of safety of 5 percent of D2 is recommended. 


2. Basin II may be effective down to a Froude number of 4 but the lower 


values should not be taken for granted (see Paper 1404 for values less than 
4.5). 
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3. The length of basin can be obtained from the intermediate curve on 
Figure 12. 


4. The height of chute blocks is equal to the depth of flow entering the ba- 
sin, or Dj, Figure 14. The width and spacing should be equal to approxi- 
mately D1; however, this may be varied to eliminate fractional blocks. A 


space equal to + is preferable along each wall to reduce spray and maintain 


desirable pressures. 


5. The height of the dentated sill is equal to 0.2D9, while the maximum 
width and spacing recommended is approximately 0.15Dg. In this case a block 
is recommended adjacent to each side wall, Figure 14. The slope of the con- 
tinuous portion of the end sill is 2:1. In the case of narrow basins, which 
would involve only a few dentates according to the above rule, it is advisable 
to reduce the width and the spacing so long as this is done proportionately. 
Reducing the width and spacing actually improves the performance in narrow 
basins, thus, the minimum width and spacing of the dentates is governed only 
by structural considerations. 


6. It is not necessary to stagger the chute blocks and the sill dentates. In 
fact this practice is usually inadvisable from a construction standpoint. 


7. The verification tests on Basin II indicated no perceptible change in the 
stilling basin action with respect to the slope of the chute preceding the basin. 
The slope of chute varied from 0.6:1 to 2:1 in these tests, Column 25, Table 3. 
Actually, the slope of the chute does have an effect on the hydraulic jump in 
some cases. This subject will be discussed in more detail in Paper 1405 with 
regard to sloping aprons. It is recommended that the sharp intersection be- 
tween chute and basin apron, Figure 14, be replaced with a curve of reason- 
able radius (R $ 4D) when the slope of the chute is 1:1 or greater. Chute 
blocks can be incorporated on the curved face as readily as on the plane sur- 
faces. On steep chutes the length of top surface on the chute blocks should be 
made sufficiently long to deflect the jet. 

Following the above rules will result in a safe, conservative stilling basin 
for spillways with fall up to 200 feet and for flows up to 500 cfs per foot of 
basin width, provided the jet entering the basin is reasonably uniform both as 
to velocity and depth. For greater falls, larger unit discharges, or possible 
asymmetry, a model study of the specific design is recommended. 


Aids in Computation 


Before presenting an example illustrating the method of proportioning Ba- 
sin II, a chart will be presented which should be of special value for prelimi- 
nary computations. The chart makes it possible to determine Vj and Dj with 
a fair degree of accuracy, for chutes having slopes of 0.8:1 or steeper, where 
computation is a difficult and arduous procedure. The chart presented as 
Figure 15 represents a composite of experience, computation, and a limited 
amount of experimental information obtained from prototype tests on Shasta 
and Grand Coulee Dams. There is much to be desired in the way of experi- 
mental confirmation; however, it is felt that this chart is sufficiently accurate 
for preliminary design. 

The ordinate on Figure 15 is the fall from reservoir level to stilling basin 
floor, while the abscissa is the ratio of actual to theoretical velocity at 
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entrance to the stilling basin. The theoretical velocity V7 = VY 2g(Z-H/2) (see 
Figure 15). The actual velocity is the term desired. The curves represent 
different heads, H, on the crest of the spillway. As is reasonable, the larger 
the head on the crest, the more nearly the actual velocity a‘ the base of the 
spillway will approach the theoretical. For example, with H = 40 feet and 

Z = 230 feet, the actual velocity at the base of the dam would be 0.95 of the 
computed theoretical velocity; while with a head of 10 feet on the crest, the 
actual velocity would be 0.75 VjT. The value of Dj may then be computed by 
dividing the unit discharge by the actual velocity obtained from Figure 15. 

The chart is not applicable for chutes flatter than 9.6:1 as frictional re- 
sistance assumes added importance in this range. Therefore, it will be ne- 
cessary to compute the hydraulic losses starting at the gate section where 
critical depth is known. 

Insufflation, produced by air from the atmosphere mixing with the sheet of 
water during the fall, need not be considered in the hydraulic jump computa- 
tions. Insufflation need be considered principally in determining the height of 
chute and stilling basin walls. It is usually not possible to construct walls 
sufficiently high to confine all spray and splash; thus, the best that can be 
hoped for is a height that is reasonable and commensurate with the material 
and terrain to be protected. 


Application of Results (Example 2) 


The crest of an overfall dam, having a downstream slope of 0.7:1, is 200 
feet above the horizontal floor of the stilling basin. The head on the crest is 
30 feet and the maximum discharge is 480 cfs per foot of stilling basin width. 
Proportion a Type ‘I stilling basin for these conditions. 


Entering Figure 15 with a head of 30 feet over the crest and a total fall of 
230 feet, 


The theoretical velocity Vr = 2g(230-5>) = 117.6 ft/sec 


The actual velocity V, = V 


r = 117.6 x 0.92 = 108.2 ft/sec 


1 


==> = 4.44 feet 


The Froude number 


Vy 108.2 


Entering Figure 11 with a Froude number of 9.04, the solid line gives 


9.04 


TW _ 123 
D, 


As TW and Dp are synonomous in this case, the conjugate tail water depth 
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BRADLEY - PETERKA 
Dy = 12.3 x 4.44 = 54.6 feet 


The minimum tail water line for the Type II basin on Figure 11 shows that 


a factor of safety of about 4 percent can be expected for the above Froude 
number. 


Should it be desired to provide a margin of safety of 7 percent, the follow- 
ing procedure may be followed: Consulting the line for minimum TW depth 
for the Type II basin, Figure 11, 


* = 11.85 for a Froude number of 9.04 
1 


The tail water depth at which sweepout is incipient: 
T= 11.85 x 4.44 = 52.6 feet 
so 


Adding 7 percent to this figure, the stilling basin apron should be posi- 
tioned for a tail water depth of 


52.6 + 3.7 = 56.3 feet or 1.03D, 


The length of basin can be obtained by entering the intermediate curve on 
Figure 12 with the Froude number of 9.04 


Ly = 4.28 x 54.6 = 234 feet (see Figure 14) 


The height, width, and spacing of the chute blocks as recommended is Dj, 
thus the dimension can be 4 feet 6 inches. 

The height of the dentated sill is 0.2Do or 11 feet, while the width and 
spacing of the dentates can be 0.15D9 or 8 feet 3 inches. 
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HYDRAULIC DESIGN OF STILLING BASINS: SHORT STILLING BASIN 
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SMALL SPILLWAYS (BASIN III) 


J. N. Bradley! M. ASCE, and A. J. Peterka2 M. ASCE 
(Proc. Paper 1403) 


FOREWORD 


This is one of a group of six papers on the hydraulic design of stilling ba- 
sins and their associated appurtenances. Although original data from hydrau- 
lic models are presented in tabular form, the papers emphasize practical de- 
sign procedures. Basin sizes and dimensions are given in easy-to-use 
dimensionless forms so that dependable stilling basins may be designed with- 
out the need for exceptional judgment or extensive experience. Sample prob- 
lems are included. 

The “Introduction” and “Experimental Equipment” are included in Paper 
1401 but apply to all papers in this group. Paper 1401 is an academic study 
of the hydraulic jump on a flat floor. Useful information is presented con- 
cerning energy losses, applicability of the hydraulic jump formula, length of 
jump, and a new classification of the various types of jumps. 

Paper 1402 covers the design of a short hydraulic jump stilling basin hav- 
ing an end sill and chute blocks. Paper 1403 describes a shorter stilling ba- 
sin which utilizes baffle piers. In both papers information is given to deter- 
mine the critical dimensions of the stilling basin, the tail water range, and 
the water surface profile in the basin. 

Paper 1404 describes a special type of hydraulic jump basin for use when 
the Froude number of the incoming flow is low, and the jump produces waves 
in the downstream channel. An alternate design and two types of wave sup- 
pressors are also developed. 

Paper 1405 describes the design of a stilling basin having a sloping apron. 
The extra tail water depth required and the merits of a sloping apron are 
evaluated, as is the basin length for a range of apron slopes. 

Paper 1406 develops the design of impact basins for discharges up to 340 


Note: Discussion open until March 1, 1958. Paper 1403 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 5, October, 1957. 


1. Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly of 
U. S. Bureau of Reclamation, Denver, Colo. 
2. Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 
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second feet and incoming velocities up to 30 feet per second. No tail water is 
required. The bibliography in this paper applies to all papers in this group. 


ABSTRACT 


Continuing the development in Paper 1402, this paper describes tests to 
generalize the design of a shorter stilling basin containing chute blocks, end 
sill, and one row of baffle piers. General design rules and sample problems 
are included. Critical basin dimensions, water surface profiles, and tail 

water ranges are presented in dimensionless forms. 


INTRODUCTION 


Basin II (paper 1402) often is considered too conservative and consequently 
overcostly for structures carrying relatively small discharges at moderate 
velocities. In this Paper a minimum basin is developed for a class of smaller 
structures in which the velocity at the entrance to the basin is moderate or 
low (5 to 60 feet per second, corresponding to an overall head of about 100 
feet). Further economies in basin length are accomplished with baffle piers. 
Verification test results are used to determine the range over which Basin II 


will operate satisfactorily, and to obtain data necessary to generalize the 
design. 


Development 


The most effective way to shorten a stilling basin is to modify the jump by 
the addition of appurtenances in the basin. One restriction imposed on these 
appurtenances, however, is that they must be self-cleaning or nonclogging. 
This restriction thus limits the appurtenances to blocks or sills which can be 
incorporated on the stilling basin apron. 

Numerous experiments were therefore performed using various types and 
arrangements of baffle piers and sills on the apron in an effort to obtain the 
best possible solution. Some of the arrangements tested are shown on 
Figure 16. The blocks were positioned in both single and double rows with 
the second row staggered with respect to the first. Arrangement “a” on 
Figure 16 consisted of a solid curved sill which was tried in several positions 
on the apron. This sill required an excessive tail water depth to be effective. 
The solid sill was then replaced with baffle piers. For certain heights, 
widths, and spacing, block “b” performed quite well, resulting in a water sur- 
face similar to that shown on Figure 19. Block “c” was ineffective for any 
height. The velocity passed over the block at about a 45° angle, thus was not 
impeded, and the water surface downstream was very turbulent with waves. 
Stepped block “d,” both for single and double rows, was much the same as 
“c.” The cube “e” was effective when the best height, width, spacing, and 
position on the apron were found. The front of the jump was almost vertical 
and the water surface downstream was quite flat and smooth, much like the 
water surface shown on Figure 19. Block “f” performed identically with the 
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cubical block “e.” The important feature as to shape appeared to be the ver- 
tical upstream face. The foregoing blocks and others not mentioned here 
were all tested in single and double rows. The second row in each case was 
of little value, sketch “h,” Figure 16. 

Block “g” is the same as block “f” with the corners rounded. It was found 
that rounding the corners greatly reduced the effectiveness of the blocks. In 
fact a double row of blocks with rounded corners did not perform as well as a 
single row of blocks “b,” “e,” or “f.” As block “f” is usually preferable from 
a construction standpoint, it was used throughout the remaining tests to de- 
termine a general design with respect to height, width, spacing, and position 
on the apron. 

In addition to experimenting with the baffle piers, variations in the size 
and shape of the chute blocks and the end sill were also tested. It was found 
that the chute blocks should be kept small, no larger than Dj if possible, to 
prevent the chute blocks from directing the flow over the baffle piers. The 
end sill had little or no effect on the jump proper when haffle piers are placed 
as recommended. Thus, there is no need for a dentated end sill and most any 
type of solid end sill will suffice. The basin as finally developed is shown on 
Figure 17. This basin is principally an impact dissipation device whereby 
the baffle piers are called upon to do most of the work. The chute blocks aid 
in stabilizing the jump and the solid type end sill is for scour control. 


Verification Tests 


At the conclusion of the development work, a set of verification tests was 
made to examine and record the performance of this basin, which will be 
designated as Basin III, over the entire range of operating conditions that may 
be met in practice. The tests were made on a total of 14 basins constructed 
in Flumes B, C, D, and E. The conditions under which the tests were run, 
the dimensions of the basin, and the results are recorded in Table 4. The 
headings are identical with those of Table 3 except for the dimensions of the 


baffle piers and end sills. The additional symbols can be identified from 
Figure 18. 


Stilling Basin Performance and Design 


Stilling basin action was quite stable for this design; in fact, more so than 
for either Basins I (paper 1401) or II (paper 1402), The front of the jump was 
steep and there was less wave action to contend with downstream than in 
either of the former basins. In addition, Basin III has a large factor of safety 
against jump sweepout and operates equally well for all values of the Froude 
number above 4.0. The verification tests served to show that Basin III was 
very satisfactory. 

Basin III should not be used where baffle piers will be exposed to veloci- 
ties above the 50 to 60 feet per second range without the full realization that 
cavitation and resulting damage may occur. For velocities above 50 feet per 
second Basin II or hydraulic model studies are recommended. 


Chute Blocks 


The recommended proportions for Basin III are shown on Figure 18. The 
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Table 3 


VERIFICATION TESTS ON TYPE II STILLING BASIN 


Q :Width: per: 
Flume:Test: cfs : of : ft : : : : : > D5 ‘sweep 
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height, width, and spacing of the chute blocks are equal to Dj, the same as 
was recommended for Basin II (paper 1402). Larger heights were tried, as 
can be observed from Column 18, Table 4, but are not recommended. The 
larger chute blocks tend to throw a portion of the high-velocity jet over the 
baffle piers. Some cases will be encountered in design, however, where Dj 

is less than 8 inches. In such cases the blocks may be made 8 inches high, 
which is considered by some designers to be the minimum size possible from 
a construction standpoint. The width and spacing are the same as the height, 
but this may be varied so long as the aggregate width of spaces approximately 
equals the total width of the blocks. 


Baffle Piers 


The height of the baffle piers increases with the Froude number as can be 
observed from Columns 22 and 10, Table 4. The height, in terms of Dj, can 
be obtained from the upper line on Figure 18. The width and spacing may be 
varied but the total of the spaces should equal the total width of blocks. The 
most satisfactory width and spacing was found to be three-fourths of the 
height. It is not necessary to stagger the baffle piers with the chute blocks 
as this is often difficult and there is little to be gained from a hydraulic 
standpoint. 

The most effective position of the baffle piers is 0.8D9 downstream from 
the chute blocks as shown in Figure 17. The actual positions used in the veri- 
fication tests are shown in Column 25, Table 4. The recommended position, 
height and spacing of the baffle piers on the apron should be adhered to care- 
fully, as these dimensions are important. For example, if the blocks are set 
appreciably upstream from the position shown they will produce a cascade 
with resulting wave action. On the other hand, if the baffles are set farther 
downstream than shown, a longer basin will be required. Likewise, if the 
baffles are too high they can produce a cascade, while if too low a rough wa- 
ter surface will result. It is not the intention to give the impression that the 
position or height of the baffle blocks are critical. Their position or height 
are not critical so long as the above proportions are followed. There exists 
a reasonable amount of leeway in all directions; however, one cannot place 
the baffle piers on the pool floor at random and expect anything like the excel- 
lent action associated with the Type III basin. 

The baffle piers may be in the form shown on Figure 17, or they may be 
cubes; either shape is effective. The corners of the baffle blocks should not 
be rounded, as the sharp edges are effective in producing eddies which in 


turn aid in the dissipation of energy. Small chambers on the pier edges are 
permissible. 


End Sill 


The height of the solid end sill is also shown to vary with the Froude num- 
ber although there is nothing critical about this dimension. The heights of the 
sills used in the verification tests are shown in Columns 27 and 28 of Table 4. 
The height of the end sill in terms of Dj is plotted with respect to the Froude 
number and shown as the lower line on Figure 18. A slope of 2:1 was used 

throughout the tests since previous sill experiments had showed that mini- 

mum wave heights and erosion could be expected with this slope. 
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As in the case of Basin II (paper 1402), full conjugate depth, measured 
above the apron, is also recommended for Basin III. There are several rea- 
sons for this: First, the best operation for this stilling basin occurs at full 
conjugate tail water depth; secondly, if less than the conjugate depth is used, 
the surface velocities leaving the pool are high, the jump action is impaired, 
and there is a greater chance for scour downstream; and thirdly, if the baffle 
blocks erode with time, the additional tail water depth will serve to lengthen 
the interval between repairs. On the other hand, there is no hydraulic advan- 
tage in using greater than the conjugate depth, as the action in the pool will 
show little or no improvement. The same precautions should be considered 
when determining the tail water for Basin III that were discussed for Basin II, 
(paper 1402). 

The margin of safety for Basin III varies from 15 to 18 percent depending 
on the value of the Froude number, as can be observed by the dashed line 
labeled “Minimum Tail Water Depth—Basin III,” on Figure 11. The points, 
from which the line was drawn, were obtained from the verification tests, 
Columns 10 and 14, Table 4. Again, this line does not represent complete 
sweepout, but the point at which the front of the jump moves away from the 
chute blocks and the basin no longer functions properly. In special cases it 
may be necessary to encroach on this wide margin of safety, however, it is 
not advisable as a general rule for the reasons stated above. 


Length of Basin 


The length of Basin II, which is related to the Froude number, can be ob- 
tained by consulting the lower curve of Figure 12, (paper 1402). The points, 
indicated by circles, were obtained from Columns 10 and 12, Table 4, and in- 
dicate the extent of the verification tests. The length is measured from the 
downstream end of the chute blocks to the downstream end of the end sill, 
Figure 17. Although this curve was determined conservatively, it will be 
found that the length of Basin III is less than one-half the length needed for a 
basin without appurtenances. Basin II, as was true of Basin II, may be ef- 
fective for values of the Froude number as low as 4.5, thus the length curve 
was terminated at this value. 


Water Surface and Pressure Profiles 


Approximate water-surface profiles were obtained for Basin II during the 
verification tests. The front of the jump was so steep, Figure 19, that only 
two measurements were necessary; the tail water depth and the depth up- 
stream from the baffle piers. The tail water depth is shown in Column 6 and 
the upstream depth is recorded in Column 29 of Table 4. The ratio of the up- 
stream depth to conjugate depth is shown in Column 30. As can be observed, 
the ratio is much the same regardless of the value of the Froude number. 
The average of the ratios in Column 30 is 0.52. Thus it will be assumed that 
the depth upstream from the baffle blocks is one-half the tail water depth. 

The profile represented by the crosshatched area, Figure 19, is for conju- 
gate tail water depth. For a greater tail water depth, Dz, the upstream depth 

D 
would be > . For a tail water depth less than conjugate, Dy, the upstream 
D 

depth would be approximately >. There appears to be no particular 
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Significance in the fact that this ratio is one-half. 

The information on Figure 19 applies only to Basin II, proportioned ac- 
cording to the rules set forth. It can be assumed that for all practical pur- 
poses the pressure and water-surface profiles are the same. There will be 
a localized increase in pressure on the apron immediately upstream from 
each baffle block but this has been taken into account, more or less, by ex- 


tending the diagram to full tail water depth beginning at the upstream face of 
the baffle blocks. 


Recommendations 


The following rules pertain to the design of the Type III basin, Figure 17: 


1. The stilling basin operates best at full conjugate tail water depth, Do. 
A reasonable factor of safety is inherent in the conjugate depth for all values 
of the Froude number (Figure 11, paper 1402) and it is recommended that this 
margin of safety not be encroached upon. 


2. The length of pool, which is less than one-half the length of the natural 
jump, can be obtained by consulting the curve for Basin III in Figure 12, (pa- 
per 1402). As a reminder, an excess of tail water depth does not substitute 
for pool length or vice versa. 


3. Stilling Basin III may be effective for values of the Froude number as 
low as 4.0 but this cannot be stated for certain (consult Paper 1404 for values 
under 4.5). 


4. Height, width, and spacing of chute blocks should equal the average depth 
of flow entering the basin, or Dj. Width of blocks may be decreased, pro- 
vided spacing is reduced a like amount. Should Dj prove to be less than 8 
inches, make the blocks 8 inches high. 


5. The height of the baffle piers varies with the Froude number and is 
given on Figure 18. The blocks may be cubes or they may be constructed as 
shown on Figure 17; the upstream face should be vertical and in one plane. 
The vertical face is important. The width and spacing of baffle piers are also 
shown on Figure 17. In narrow structures where the specified width and spac- 
ing of blocks do not appear practical, block width and spacing may be reduced, 


provided they are reduced a like amount. A half space is recommended adja- 
cent to the walls. 


6. The upstream face of the baffle piers should be set at a distance of 


0.8D2 from the downstream face of the chute blocks (Figure 17). This dimen- 
sion is also important. 


7. The height of the solid sill at the end of the basin can be obtained from 
Figure 18. The slope is 2:1 upward in the direction of flow. 


8. It is undesirable to round or streamline the edges of the chute blocks, 
end sill or baffle piers. Streamlining of baffle piers may result in loss of half 
of their effectiveness. Small chamfers to prevent chipping of the edges is 
permissible. 


9. It is recommended that a radius of reasonable length (R 5 4D ) be used 
at the intersection of the chute and basin apron for slopes of 45° or greater. 
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10. As a general rule the slope of the chute has little effect on the jump 
unless long flat slopes are involved. This phase will be considered in Paper 
1405 on sloping aprons. 


Since Basin III is close-coupled, the above rules should be followed closely 
for its proportioning. If the proportioning is to be varied from that recom- 
mended, or if the limits given below are exceeded, a model study is advisable. 
Arbitrary limits for the Type III basin are set at 200 cfs per foot of basin 
width and 50 to 60 feet per second entrance velocity until experience demon- 
strates otherwise. 


Application of Results (Example 3) 


Given the following computed values for a small overflow dam 


cfs 


3,900 78.0 

3,090 61.8 

2,022 40.45 
662 13.25 


and the tail water curve for the river, identified by the solid line on Figure 20, 
proportion Basin III for the most adverse condition utilizing full conjugate tail 
water depth. The flow is symmetrical and the width of the basin is 50 feet. 
(The purpose of this example is to demonstrate the use of the jump elevation 
curve.) 

The ffrst step is to compute the jump elevation curve which in this case 
is D2 plus the elevation of the basin floor. As V, and Dj are given, the 
Froude number is computed and tabulated in Column 2, Table 5, below: 


Table 5 


ae Junp elevation 
ft Curve a Curve a' 


wy ~ 


1.130 617.5 615.0 
0.936 615.2 612.7 
0.€42 612.0 609.5 
0.260 606.1 603.6 


Entering Figure 11 (paper 1402) with these values of the Froude number values 


of ™ are obtained for conjugate tail water depth from the solid line. These 
D 


values are also . and are shown listed in Column 3. The conjugate tail water 


D, 


depths for the various discharges, Column 5, were obtained by multiplying the 
values in Column 3 by those in Column 4. 

If it were assumed that the most adverse operating condition occurs at the 
maximum discharge of 3,900 cfs, the stilling basin apron should be placed at 
elevation 617.5 - 17.8 or elevation 599.7 


| 
Q 4 V1 Di fo 
cfs = it/sec ft_ 
69 1.130 
66 0.936 
63 0.642 
51 0.260 
| 
| 
cfs Fi Dy 
3,900 11.42 15.75 
3,090 12.02 16.60 
2,022 13.85 19.20 
662 17.62 24.5 
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With the apron at elevation 599.7 the tail water required for conjugate tail 
water depth for each discharge would follow the elevations listed in Column 6. 
Plotting Columns 1 and 6 on Figure 20 results in Curve a, which shows that 
the tail water depth is inadequate for all but the maximum discharge. 

The tail water curve is unusual in that the most adverse tail water condi- 
tion occurs at a discharge of approximately 2,850 cfs rather than maximum. 
As full conjugate tail water depth is desired for the most adverse tail water 
condition, it is necessary to shift the jump elevation curve downward to match 
the tail water curve for a discharge of 2,850 cfs (see Curve a' , Figure 20). 
The coordinates for Curve a' are given in Columns 1 and 7, Table 5. This 
will place the basin floor 2.5 feet lower, or elevation 597.2 feet, as shown in 
the sketch on Figure 20. 

Although the position of the basin floor was set for a discharge of 2,850 
cfs, the remaining details are proportioned for the maximum discharge 
3,900 cfs. 

Entering Figure 12 with a Froude number of 11.42, 

Lin 


—— = 2.75, and the length of 
D2 


basin required Lyy = 2.75 x 17.80 = 48.95 feet. 


(Notice that conjugate depth was used, not tail water depth.) 


The height, width, and spacing of chute blocks are equal to Dy or 1.130 feet 
(use 13 or 14 inches). 


The height of the baffle piers for a Froude number of 11.42 (Figure 19) is 


2.5D 1. 
hg = 2.5 x 1.130 = 2.825 feet (use 34 inches). 


The width and spacing of the baffle piers are preferably three-fourths of 
the height or 


0.75 x 34 = 25.5 inches. 


From Figure 17, the upstream face of the baffle piers should be 0.8D9 
from the downstream face of the chute blocks, or 


0.8 x 17.80 = 14.24 feet. 
The height of the solid end sill, Figure 18, is 1.60D , or 
hy = 1.60 x 1.130 = 1.81 feet (use 22 inches). 


The final dimensions of the basin are shown on Figure 20. 
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FOREWORD 


This is one of a group of six papers on the hydraulic design of stilling ba- 
sins and their associated appurtenances. Although original data from hydrau- 
lic models are presented in tabular form, the papers emphasize practical de- 
sign procedures. Basin sizes and dimensions are given in easy-to-use 
dimensionless forms so that dependable stilling basins may be designed with- 
out the need for exceptional judgment or extensive experience. Sample prob- 
lems are included. 

The “Introduction” and “Experimental Equipment” are included in Paper 
1401 but apply to all papers in this group. Paper 1401 is an academic study 
of the hydraulic jump on a flat floor. Useful information is presented con- 
cerning energy losses, applicability of the hydraulic jump formula, length of 
jump, and a new classification of the various types of jumps. 

Paper 1402 covers the design of a short hydraulic jump stilling basin hav- 
ing an end sill and chute blocks. Paper 1403 describes a shorter stilling ba- 
sin which utilizes baffle piers. In both papers information is given to deter- 
mine the critical dimersions of the stilling basin, the tail water range, and 
the water surface profile in the basin. 

Paper 1404 describes a special type of hydraulic jump basin for use when 
the Froude number of the incoming flow is low, and the jump produces waves 
in the downstream channel. An alternate design and two types of wave sup- 
pressors are also developed. 

Paper 1405 describes the design of a stilling basin having a sloping apron. 
The extra tail water depth required and the merits of a sloping apron are 
evaluated, as is the basin length for a range of apron slopes. 

Paper 1406 develops the design of impact basins for discharges up to 340 


Note: Discussion open until March 1, 1958. Paper 1404 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 5, October, 1957. 


1. Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly of 
U. S. Bureau of Reclamation, Denver, Colo. 
2. Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 
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second feet and incoming velocities up to 30 feet per second. No tail water is 
required. The bibliography in this paper applies to all papers in this group. 


ABSTRACT 


Where the Froude number of the incoming flow is low, the hydraulic jump 
produces objectionable waves downstream. This study generalizes the design 
of two stilling basins and two types of wave suppressors for stabilizing the 
water surface. Sample problems are used to illustrate the design procedures 
recommended. 


INTRODUCTION 


In this Paper the characteristics of the hydraulic jump and the design of 
an adequate stilling basin for Froude numbers between 2.5 and 4.5 are dis- 
cussed. This range is encountered principally in the design of canal struc- 
tures, but occasionally !sw dams and outlet works fall in this category. In 
the 2.5 to 4.5 Froude number range, the jump is not fully developed and the 
previously discussed methods of design do not apply (papers 1401-1403). The 
main problem concerns the waves created in the hydraulic jump, making the 
design of a suitable wave suppressor a part of the stilling basin problem. 

Four means of reducing wave heights are discussed. The first is an in- 
tegral part of the stilling basin design and should be used only in the 2.5 to 
4.5 Froude number range. The second may be considered to be an alternate 
design and may be used over a greater range of Froude numbers. These 
types are discussed as a part of the stilling basin design. The third and 
fourth devices are considered as appurtenances which may be included in an 
original design or added to an existing structure. Also, they may be used in 
any open channel flow-way without consideration of the Froude number. 
These latter devices are described under the heading Wave Suppressors. 


Jump Characteristics—Froude Number 2.5 to 4.5 


For low values of the Froude number, 2.5 to 4.5, the entering jet oscillates 
intermittently from bottom to surface, as indicated in Figure 9B, (paper 1401) 
with no particular period. Each oscillation generates a wave which is difficult 
to dampen. In narrow structures, such as canals, waves may persist to some 
degree for miles. As they encounter obstructions in the canal, such as bridge 
piers, turnouts, checks, and transitions, reflected waves may be generated 
which tend to dampen, modify, or intensify the original wave. Waves are 
destructive to earth-lined canals and riprap and produce undesirable surges 
at gaging stations and in measuring devices. Structures in this range of 
Froude numbers are the ones which have been found to require the most 
maintenance. 

On wide structures, such as diversion dams, wave action is not as pro- 
nounced when the waves can travel laterally as well as parallel to the direc- 
tion of flow. The combined action produces some dampening effect but also 
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results in a choppy water surface. These waves may or may not be dissipat- 
ed in a short distance. Where outlet works, operating under heads of 50 feet 
or greater, fall within the range of Froude numbers between 2.5 and 4.5, a 
model study of the stilling basin is imperative. A model study is the only 
means of including preventive or corrective devices in the structure so that 
proper performance can be assured. 


Stilling Basin Design—Froude Number 2.5 to 4.5 


Development Tests 


The best way to combat a wave problem is to eliminate the wave at its 
source; in other words, concentrate on altering the condition which generates 
the wave. In the case of the stilling basin preceded by an overfall or chute, 
two schemes were apparent for eliminating waves at their source. The first 
was to break up the entering jet by opposing it with directional jets deflected 
from baffle piers or sills. The second was to bolster or intensify the roller, 
shown in the upper portion of Figure 9B, (paper 1401) by directional jets de- 
flected from large chute blocks. 

The first method was unsuccessful in that the number and size of appur- 
tenances necessary to break up the jet occupied so much volume that the de- 
vices themselves posed an obstruction to the flow. This conclusion was based 
on tests in which various shaped baffle and guide blocks were systematically 
placed in a stilling basin in combination with numerous types of spreader 
teeth and deflectors in the chute. The program involved dozens of tests, and 
not until all conceivable ideas were tried was this approach abandoned. A 
few of the basic ideas tested are shown on Figure 21, a, b, c, f, g, and h. 


Final Tests 


Deflector Blocks 

The second approach, that of attempting to intensify the roller, yielded bet- 
ter results. In this case, large blocks were placed well up on the chute, while 
nothing was installed in the stilling basin proper. The object was to direct a 
jet into the base of the roller in an attempt to strengthen it. After a number 
of trials, the roller was actually intensified which did improve the stability of 
the jump. Sketches d and e on Figure 21 indicate the only schemes that 
showed promise, although many variations were tried. After finding an ar- 
rangement that was effective, approximations were tested to make the field 
construction as simple as possible. The dimensions and proportions of the 
adopted deflector blocks are shown on Figure 22. 

The object in the latter scheme was to place as few appurtenances as pos- 
sible in the path of the flow, as volume occupied by appurtenances helps to 
create a backwater problem, thus requiring higher training walls. The num- 
ber of deflector blocks shown on Figure 22 is a minimum requirement to ac- 
complish the purpose set forth. The width of the blocks is shown equal to Dj 
and this is the maximum width recommended. From a hydraulic standpoint it 
is desirable that the blocks be constructed narrower than indicated, prefera- 
bly 0.75D,. The ratio of block width to spacing should be maintained as 1:2.5. 
The extreme tops of the blocks are 2D, above the floor of the stilling basin. 
The blocks may appear to be rather high and, in some cases, extremely long, 
but this is essential as the jet must play at the base of the roller to be 
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effective. To accommodate the various slopes of chutes and ogee shapes en- 
countered, a rule has been established that the horizontal length of the blocks 
should be at least 2D,. The upper surface of each block is sloped at 5° ina 
downstream direction as it was found that this feature resulted in better oper- 
ation, especially at the discharges lower than the design flow. 


Tail Water Depth 

A tail water depth 5 to 10 percent greater than the conjugate depth is 
strongly recommended for the above basin. Since the jump is very sensitive 
to tail water depth at these low values of the Froude number, a slight de- 
ficiency in tail water depth may allow the jump to sweep completely out of the 
basin. The jump performs much better and wave action is diminished if the 
tail water depth is increased to approximately 1.1D9. 


Basin Length and End Sill 

The length of this basin, which is relatively short, can be obtained from 
the upper curve on Figure 12 (paper 1402). No additional blocks or appur- 
tenances are needed in the basin, as these will prove a greater detriment than 
aid. The addition of a small triangular sill placed at the end of the apron for 
scour control is desirable. A sill of the type used on Basin III may also be 
used, Figure 18 (paper 1403). 


Performance 

If designed for the maximum discharge, this stilling basin will perform 
satisfactorily for ali flows. Waves below the stilling basin will still be in 
evidence but will be of the ordinary variety usually encountered with jumps of 


a higher Froude number. This design is applicable to rectangular cross sec- 
tions only. 


Alternate Stilling Basin Design—Small Drops 


Performance 


An alternate basin for reducing wave action at the source, for values of the 
Froude number between 2.5 and 4.5, is particularly applicable to small drops 
in canals. The Froude number in this case is computed the same as though 
the drop were an overflow crest. A series of steel rails, channel irons or 
timbers in the form of a grizzly are installed at the drop, as shown in Figure 
23. The overfalling jet is separated into a number of long, thin sheets of wa- 
ter which fall nearly vertically into the canal below. Energy dissipation is 
excellent and the usual wave problem is avoided. If the rails are tilted down- 
ward at an angle of 3° or more, the grid is self-cleaning. 

The use of this device is particularly justified when the Froude number is 
below 3.0. If use of a jump were possible the maximum energy loss would be 
less than 27 percent at best, Figure 8, (paper 1401). The suggested device 


accomplishes nearly as much energy loss and provides a smooth water sur- 
face in addition. 


Design 


Two spacing arrangements of the beams were tested in the laboratory: in 
the first, the spacing was equal to the width of the beams, and in the second, 
the spacing was two-thirds of the beam width. The latter was the more effec- 
tive. In the first, the length of beams required was about 2.9 times the depth 
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of flow (y) in the canal upstream, while in the second, it was necessary to in- 
crease the length to approximately 3.6y. The following expression can be 
used for computing the length of beams: 


L= 


(4) 
CSN ,/2gy 


where Q is the total discharge in cfs, C is an experimental coefficient, S is 
the width of a space in feet, N is the number of spaces, g is the acceleration 
of gravity and y is the depth of flow in the canal upstream (see Figure 23). 
The value of C for the two arrangements tested was 0.245. 

Should it be desired to maintain a certain level in the canal upstream, the 
grid may be made adjustable and tilted upward to act as a check; however, 
this arrangement may pose a cleaning problem. 


Wave Suppressors 


The two stilling basins described above may be considered to be wave sup- 
pressors, although the suppressor effect is obtained from the necessary fea- 
tures of the stilling basin. If greater wave reduction is required on a pro- 
posed structure, or if a wave suppressor is required to be added to an existing m 
flow-way, the two types discussed below may prove useful. Both are ap- ba 
plicable to most open channel flow-ways having rectangular, trapezoidal, or 
other cross-sectional shapes. The first or raft type may prove more eco- 
nomical than the second or underpass type, but rafts do not provide the degree 
of wave reduction obtainable with the underpass type. Both types may be used 
without regard to the Froude number. 


Raft Type Wave Suppressor 


In a structure of the type shown in Figure 24, there are no means for 
eliminating waves at their source. Tests showed that appurtenances in the 
stilling basin merely produced severe splashing and created a backwater ef- 
fect, resulting in submerged flow at the gate for the larger flows. Submerged 
flow reduced the effective head on the structure, and in turn, the capacity. 
Tests on several suggested devices showed that rafts provided the best an- 
swer to the wave problem when additional submergence could not be tolerated. 
The general arrangement of the tested structure is shown in Figure 24. The 
Froude number varied from 3 to 7, depending on the head behind the gate and 
the gate opening. Velocities in the canal ranged from 5 to 10 feet per second. 
Waves were 1.5 feet high, measured from trough to crest. 

During the course of the experiments a number of rafts were tested; thick 
rafts with longitudinal slots, thin rafts made of perforated steel plate, and 
others, both floating and fixed. Rigid and articulated rafts were tested in 
various arrangements. 

The most effective raft arrangement consisted of two rigid stationary rafts 
20 feet long by 8 feet wide, made from 6- by 8-inch timbers, placed in the 
canal downstream from the stilling basin (Figure 24). A space was left be- 
tween timbers and lighter cross pieces were placed on the rafts parallel to the 
flow, giving the appearance of many rectangular holes. Several essential re- 
quirements for the raft were apparent: (1) that the rafts be perforated ina 
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regular pattern; (2) that there be some depth to these holes; (3) that at least 
two rafts be used; and (4) that the rafts be rigid and held stationary. 

It was found that the ratio of hole area to total area of raft could be from 
1:6 to 1:8. The 8-foot width, W on Figure 24, is a minimum dimension. The 
rafts must have sufficient thickness so that the troughs of the waves do not 
break free from the underside. The top surfaces of the rafts are set at the 
mean water surface in a fixed position so that they cannot move. Spacing be- 
tween rafts should be at least three times the raft dimension, measured paral- 
lel to the flow. The first raft decreases the wave height about 50 percent, 
while the second raft effects a similar reduction. Surges over the raft dissi- 
pate themselves by flow downward through the holes. For this specific case 
the waves were reduced from 18 to 3 inches in height. 

Under certain conditions wave action is of concern only at the maximum 
discharge when freeboard is endangered, so the rafts can be a permanent in- 
Stallation. Should it be desired to suppress the waves at partial flows, the 
rafts may be made adjustable, or, in the case of trapezoidal channels, a 
second set of rafts may be placed under the first set for partial flows. The 
rafts should perform equally well in trapezoidal as well as rectangular 
channels. 

The recommended raft arrangement is also applicable for suppressing 
waves with a regular period such as wind waves, waves produced by the start- 
ing and stopping of pumps, etc. In this case, the position of the downstream 
raft is important. The second raft should be positioned downstream at some 
fraction of the wave length. Placing it at a full wave length could cause both 
rafts to be ineffective. Thus, for narrow canals it may be advisable to make 
the second raft portable. However, if it becomes necessary to make the rafts 
adjustable or portable, or if a moderate increase in depth in the stilling basin 
can be tolerated, consideration should be given to the type of wave suppressor 
discussed below. 


Underpass Type Wave Suppressor 


General Description 
By far the most effective wave dissipator is the short-tube type of under- 
pass suppressor. The name “short-tube” is used because the structure has 


many of the characteristics of the short-tube discussed in hydraulic textbooks. 


This wave suppressor may be added to an existing structure or included in the 
original construction. In either case it provides a sightly structure, perman- 
ent in nature, which is economical to construct and effective in operation. 

Essentially, the structure consists of a horizontal roof placed in the flow 
channel with a headwall sufficiently high to cause all flow to pass beneath the 
roof. The height of the roof above the channel floor may be set to effectively 
reduce wave heights for a considerable range of flows or channel stages. 

The length of the roof, however, determines the amount of wave suppression 
obtained for any particular roof setting. 

The recommendations for this structure are based on three separate model 
investigations, each having different flow conditions and wave reduction re- 
quirements. The design is then generalized and design procedures given, in- 
cluding a sample problem. 


Performance 
The effectiveness of the underpass wave suppressor is illustrated in 
Figures 25 and 26. Figure 25 shows one of the hydraulic models used to 
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develop the wave suppressor and the effect of the suppressor on the waves in 

the canal. Figure 26 shows before and after photographs of the prototype in- 
stallation, indicating that the prototype performance was as good as predicted 

by the model. In this instance it was desired to reduce wave heights entering { 
a lined canal to prevent overtopping of the canal lining at near maximum dis- t 
charges. Below 3,000 second-feet, waves were in evidence but did not overtop 
the lining. For larger discharges, however, the stilling basin produced 
moderate waves which were actually intensified by the short transition be- * { 
tween the basin and the canal. These intensified waves overtopped the lining 
at 4,000 second-feet and became a real problem at 4,500 second-feet. Tests 
were made with a suppressor 21 feet long using discharges from 2,000 to 
5,000 second-feet. The suppressor was located downstream from the stilling 
basin. 

Figure 27, Test 1, shows the results of tests to determine the optimum 
opening between the roof and the channel floor using the maximum discharge 
5,000 second-feet. With a 14-foot opening, waves were reduced from about 
8 feet to about 3 feet. Waves were reduced to less than 2 feet with an opening 
of 11 feet. Smaller openings produced less wave height reductions, due to the 
turbulence created at the underpass exit. Thus, it may be seen that an open- 
ing of from 10 to 12 feet produced optimum results. 

With the opening set at 11 feet the suppressor effect was then determined 
for other discharges. These results are shown on Figure 27, Test 2. Wave 
height reduction was about 78 percent at 5,000 second-feet, increasing to § 
about 84 percent at 2,000 second-feet. The device became ineffective at about 7 
1,500 second-feet when the depth of flow became less than the height of the 
roof. 

To determine the effect of suppressor length on the wave reduction, other 
factors were held constant while the length was varied. Tests were made on 
suppressors 10, 21, 30, and 40 feet long for discharges of 2, 3, 4, and 5 thou- 
sand second-feet, Figure 27, Test 3. Roof lengths in terms of the downstream 
depth, Dg, for 5,000 second-feet were 0.62D2, 1.31D2, and 2.5D2, respectively. } 
In terms of a 20-foot-long underpass, halving the roof length almost doubled 
the downstream wave height while doubling the 20-foot length almost halved 
the resulting wave height. 

The same type of wave suppressor was successfully used in an installation 
where it was necessary to obtain optimum wave height reductions, since flow 
from the underpass discharged directly into a measuring flume in which it 
was desired to obtain accurate discharge measurements. The capacity of the 
structure was 625 cubic feet per second but it was necessary for the under- 
pass to function for low flows as well as for the maximum. With an underpass 
3.5D2 long and set as shown in Figure 28, the wave reductions were as shown 
in Table 6. 

Figure 28 shows actual wave traces recorded by an oscillograph. Here it 
may be seen that the maximum wave height, measured from minimum trough 
to maximum crest did not occur on successive waves. Thus, the water sur- 
face will appear smoother to the eye than is indicated by the maximum wave 
heights recorded in Table 6. 


{ 


General Design Procedure 

To design an underpass for a particular structure there are three main 
considerations: (1) how deeply should the roof be submerged, (2) how long an 
underpass should be constructed to accomplish the necessary wave reduction, 
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Table 6 


WAVE HEIGHTS IN FEET--PROTOTYPE 
Maximum Head 
Discharge : 625 : > 4OO : 200 : 
in cfs :Upstream*:Downstream*: U: D: U: : : 


Wave heights: **3.8 0.3 
in feet lus : : : : : : : 


and (3) how much increase in flow depth will occur upstream from the under- 
pass. These considerations are discussed in order. 

Based on the two installations shown in Figures 27 and 28, and on other 
experiments, it has been found that maximum wave reduction occurs when the 
roof is submerged about 33 percent, i.e., when the under side of the underpass 
is set 33 percent of the flow depth below the water surface for maximum dis- 
charge, Figure 29C. Submergences greater than 33 percent (for the cases 
tested) produced undesirable turbulence at the underpass outlet resulting in 
less overall wave reduction. With the usual tail water curve, submergence 
and the percent reduction in wave height will become less, in general, for 
smaller than maximum discharges. This is illustrated by the upper curve in 
Figure 29C. The lower curve shows a near constant value for less submer- 
gence, but it is felt that this is a somewhat special case since the wave heights 
for less than maximum discharge were smaller and of shorter period than in 
the usual case. 

It is known that the wave period greatly affects the performance of a given 
underpass, with the greatest wave reduction occurring for short period waves. 
Since the wave periods to be expected are usually not known in advance, it is 
desirable to eliminate this factor from consideration. Fortunately, wave ac- 
tion below a stilling basin usually has no measurable period but consists of a 
mixture of generated and reflected waves best described as a choppy water 
surface. This fact makes it possible to provide a practical solution from 
limited data and to eliminate the wave period from consideration except in 
this general way: waves must be of the variety ordinarily found downstream 
from hydraulic jumps or energy dissipators. These usually have a period of 
not more than about 5 seconds. Longer period waves may require special 
treatment not covered in this discussion. Fortunately, too, there generally is 
a tendency for the wave period to become less with decreasing discharge. 
Since the suppressor provides a greater percentage reduction on shorter 
period waves, this tends to offset the characteristics of the device to give less 
wave reduction for reduced submergence at lower discharges. It is therefore 
advisable in the usual case to submerge the underpass about 33 percent for 
the maximum discharge. For less submergence, the wave reduction can be 
estimated from Figure 29C. 

The minimum length of underpass required depends on the amount of wave 
reduction considered necessary. If it is sufficient to obtain a nominal reduc- 
tion to prevent overtopping of a canal lining at near maximum discharge or to 
prevent waves from attacking channel banks, a length 1D2 to 1.5D2 will 
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provide from 60 to 75 percent wave height reduction, provided the initial 
waves have periods up to about 5 seconds. The shorter the wave period the 
greater the reduction for a given underpass. For long period waves, little 
wave reduction may occur because of the possibility of the wave length being 
nearly as long or longer than the underpass, with the wave passing ‘un- 
touched” beneath the underpass. 

To obtain greater than 75 percent wave reductions, a longer underpass is 
necessary. Under ideal conditions an underpass 2Dg9 to 2.5Dg in length may 
provide up to 88 percent wave reduction for wave periods up to about 5 
seconds. Ideal conditions include a velocity beneath the underpass of less 
than, say, 10 feet per second and a length of channel 3 to 4 times the length of 
the underpass downstream from the underpass which may be used as a quiet- 
ing pool to still the small turbulence waves created at the underpass exit. 

Wave height reduction up to about 93 percent may be obtained by using an 
underpass 3.5D9 to 4D2 long. Included in this length is a 4:1 sloping roof ex- 
tending from the underpass roof elevation to the tail water surface. The 
sloping portion should not exceed about one-quarter of the total length of un- 
derpass. Since slopes greater than 4:1 do not provide the desired draft tube 
action they should not be used. Slopes flatter than 4:1 provide better draft 
tube action and are therefore desirable. 

Since the greatest wave reduction occurs in the first D2 of underpass 
length, it may appear advantageous to construct two short underpasses rather 
than one long one. In the one case tested, two underpasses each 1D2 long, 
with a length 5D2 between them, gave an added 10 percent wave reduction ad- 
vantage over one underpass 2D9 long. The extra cost of another headwall 
should be considered, however. 

Table 7 summarizes the amount of wave reduction obtainable for various 
underpass lengths. 

To determine the backwater effect of placing the underpass in the channel, 
Figure 29B will prove helpful. Data from four different underpasses were 


Table 7 
EFFECT OF UNDERPASS LENGTH ON WAVE REDUCTION 


For Underpass Submergence 33 Percent and 
Maximum Velocity Beneath 14 ft/sec 


Underpass length | Percent wave reduction* 


1Dp to 1.5D. 
to 2.5Do 
3.5 to 4.0D5 


60 to 75 


80 to 88 


90 to 93** 


*For wave periods up to 5 seconds. 
**Upper limit only with draft tube type 


outlet. 
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used to obtain the two curves shown. Although the test points from which the 
curves were drawn showed minor inconsistencies, probably because factors 
other than those considered also affected the depth of water upstream from 
the underpass, the submitted curves are sufficiently accurate for design pur- 
poses. Figure 29B shows two curves of the discharge coefficient “C” versus 
average velocity beneath the underpass, one for underpass lengths of 1D9 to 
2D2 and the other for lengths 3D2 to 4D2. Intermediate values may be inter- 
polated although accuracy of this order is not usually required. 

Pressures on the underpass were measured by means of piezometers to 
determine the direction and magnitude of the forces acting. Average pres- 
sures on the headwall were found to be about static, with distribution in a 
straight line from zero at the water surface to static pressure at the bottom. 
Pressures along the underside of the roof were found to be 1 to 2 feet below 
atmospheric; for design purposes they may be considered to be atmospheric. 
Pressures on the downstream vertical wall were equal to static pressures. 

In other words there is only a slight tendency (except for the force of breaking 
waves which was not measured) to move the underpass downstream, and there 
is a slight resultant force tending to hold the underpass down. 


Sample Problem, Example 4 

To illustrate the use of the preceeding data in designing an underpass, a 
sample problem will be helpful. 

A rectangular channel 30 feet wide and 14 feet deep flows 10 feet deep at 
maximum discharge, 2,400 cfs. It is estimated that waves will be 5 feet high 
and of the ordinary variety having a period less than 5 seconds. It is desired 
to reduce the height of the waves to approximately 1 foot at maximum dis- 
charge by installing an underpass-type wave suppressor without increasing 
the depth of water upstream from the underpass more than 15 inches. 

To obtain maximum wave reduction at maximum discharge, the underpass 
should be submerged 33 percent. Therefore, the depth beneath the underpass 
is 6.67 feet with a corresponding velocity of 12 ft/sec, 2,400 To 

(V = 30x 6.67" 
reduce the height of the waves from 5 feet to 1 foot, an 80 percent reduction 
in wave height is indicated, and, from Table 7, requires an underpass approxi- 
mately 2D9 in length. 

From Figure 29B, C = 1.07 for 2D2 and a velocity of 12 ft/sec. 

From the equation given on Figure 29B: 


Total head, h + hy = on. = (Gore 2 = 1.95 feet 
CA y 2g 


h + hy is the total head required to pass the flow and h represents the back- 
water effect of increase in depth of water upstream from the underpass. The 
determination of values for h and hy is done by trial and error. As a first 
determination, assume that h + hy represents the increase in head. 


Then, channel approach velocity, V, = g = 


2,400 
To + 1.95730 * 6.7 ft/sec 


| 
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and h = 1.95 - 0.70 #2 1.25 feet 


To refine the calculation, the above computation is repeated using the new 
head 


2,400 


vy, = Tio + 1.25)30 * 7.1 ft/sec 


= 0.72 foot 


and h = 1.17 feet 


Further refinement is unnecessary. 

Thus, the average water surface upstream from the underpass is 1.2 feet 
higher than the tail water which satisfies the assumed design requirement of 
a maximum backwater of 15 inches. The length of the underpass is 2D9 or 20 
feet, and the waves are reduced 80 percent to a maximum height of approxi- 
mately 1 foot. 

If it is desired to reduce the wave heights still further, a longer underpass 
is required. Using Table 7 and Figure 29B as in the above problem, an under- 
pass 3.5 to 4.0D9 or 35 to 40 feet in length reduces the waves 90 to 93 percent, 
making the downstream waves approximately 0.5 foot high and creating a 
backwater, h, of 1.61 feet. 

In providing freeboard for the computed backwater, h, allowance should be 
made for waves and surges which, in effect, are above the computed water 
surface. One-half the wave height or more, measured from crest to trough, 
should be allowed above the computed surface. Full wave height would pro- 
vide a more conservative design for the usual short period waves encountered 
in flow channels. 

The headwall of the underpass should be extended to this same height and a 
seawall overhang placed at the top to turn wave spray back into the basin. An 
alternate method would be to place a cover, say 2D2 long, upstream from the 
underpass headwall. 

To insure obtaining the maximum wave reduction for a given length of un- 
derpass, a 4:1 sloping roof should be provided at the downstream end of the 
underpass, as indicated on Figure 28. This slope may be considered as part 
of the overall length. The sloping roof will help reduce the maximum wave 
height and will also reduce the frequency with which it occurs, providing in 
all respects a better appearing water surface. If the flow entering the under- 
pass contains entrained air in the form of rising air bubbles, a few small 
vents in the underpass roof will relieve the possibility of air spurts and re- 
sulting surface turbulence at the underpass exit. 

If the underpass is to be used downstream from a stilling basin the under- 
pass must be placed sufficiently downstream to prevent turbulent flow, such 
as occurs at the end of a basin, from entering and passing through the wave 
suppressor. In highly turbulent flow the underpass is only partly effective. 

A close inspection of the submitted data will reveal that slightly better re- 
sults were obtained in the tests than are claimed in the example. This was 
done to illustrate the degree of conservatism required, since it should be 


understood that the problem of wave reduction can be very complex if unusual 
conditions prevail. 
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FIGURE 
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FIGURE 22 
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FIGURE 24 


HYDRAULIC JUMP STUDIES 
RAFT WAVE SUPPRESSOR (TYPE IZ) 
FOR FROUDE NUMBERS 2.5- 4.5 
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Figure 25 


Suppressor in place - Length 1.3D2, submerged 
30 percent 


Performance of Underpass Wave Suppressor 
1:32 Scale Model 
Discharge 5, 000 Second-feet 
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Friant-Kern Canal - Q = 3900 c.f.s. 
Before wave suppressor was installed 


Friant-Kern Canal - Q = 3900 c.f.s. 
After wave suppressor was installed 
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FIGURE 27 
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FIGURE 29 
A = Flow area beneath underpass 
h = Flow-producing head 
» hy= Velocity head in approach 
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HYDRAULIC DESIGN OF STILLING BASINS: STILLING BASIN WITH 
SLOPING APRON (BASIN V) 
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FOREWORD 


This is one of a group of six papers on the hydraulic design of stilling ba- 
sins and their associated appurtenances. Although original data from hydrau- 
lic models are presented in tabular form, the papers emphasize practical de- 
sign procedures. Basin sizes and dimensions are given in easy-to-use 
dimensionless forms so that dependable stilling basins may be designed with- 
out the need for exceptional judgment or extensive experience. Sample prob- 
lems are included. 

The “Introduction” and “Experimental Equipment” are included in Paper 
1401 but apply to all papers in this group. Paper 1401 is an academic study 
of the hydraulic jump on a flat floor. Useful information is presented con- 
cerning energy losses, applicability of the hydraulic jump formula, length of 
jump, and a new classification of the various types of jumps. 

Paper 1402 covers the design of a short hydraulic jump stilling basin hav- 
ing an end sill and chute blocks. Paper 1403 describes a shorter stilling ba- 
sin which utilizes baffle piers. In both papers information is given to deter- 
mine the critical dimensions of the stilling basin, the tail water range, and 
the water surface profile in the basin. 

Paper 1404 describes a special type of hydraulic jump basin for use when 
the Froude number of the incoming flow is low, and the jump produces waves 
in the downstream channel. An alternate design and two types of wave sup- 
pressors are also developed. 

Paper 1405 describes the design of a stilling basin having a sloping apron. 
The extra tail water depth required and the merits of a sloping apron are 
evaluated, as is the basin length for a range of apron slopes. 

Paper 1406 develops the design of impact basins for discharges up to 340 


Note: Discussion open until March 1, 1958. Paper 1405 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers , 
Vol. 83, No. HY 5, October, 1957. 


1. Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly of 
U. S. Bureau of Reclamation, Denver, Colo. 


2. Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 
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second feet and incoming velocities up to 30 feet per second. No tail water is 
required. The bibliography in this paper applies to all papers in this group. 


ABSTRACT 


Procedures and rules for the design of a stilling basin with a sloping apron 
are presented along with a discussion of the relative merits of sloping and 
horizontal aprons. Critical basin dimensions and tail water requirements are 
presented in dimensionless form. Recommendations are verified from an 

analysis of 13 existing basins having sloping aprons. 


INTRODUCTION 


Much has been said, pro and con, concerning the advantages and disadvan- 
tages of stilling basins with sloping aprons. Discussions have continued in- 
definitely because there were not sufficient supporting data available from 
which to draw conclusions. In this study, therefore, the sloping apron was 
investigated sufficiently to answer many of the debatable questions and also 
to provide more definite design data. 

Four flumes, A, B, D, and F, Figures 1, 2, and 3, (paper 1401) were used 
to obtain the range of Froude numbers desired for the tests. In Flumes A, B, 
and D, floors were installed to the slope desired, while Flume F could be 
tilted to obtain slopes from 0° to 12°. The slope in this discussion is the 
tangent of the angle between the floor and the horizontal and is designated as 
“tan @.” Five principal measurements were made in these tests, namely: 
the discharge, the average depth of flow entering the jump, the length of the 
jump, the tail water depth, and the slope of the apron. The tail water was ad- 
justed so that the front of the jump formed either at the intersection of the 
spillway face and the sloping apron or, in the case of the tilting flume, at a 
selected point. 

The jump on a sloping apron takes many forms depending on the slope and 
arrangement of the apron, the value of the Froude number, and the concentra- 
tion of flow (discharge per foot of width), but the dissipation is as effective as 
occurs in the true hydraulic jump on a horizontal apron. 


Previous Experimental Work 


Previous experimental work on the sloping apron has been carried on by 
several experimenters. In 1934, the late C. L. Yarnell of the United States 
Department of Agriculture supervised a series of experiments on the hydrau- 
lic jump on sloping aprons. Carl Kindsvater(5) later compiled these data and 
presented a rather complete picture, both experimentally and theoretically, 
for one slope, namely: 1:6 (tan ¢ = 0.167). G. H. Hickox(5) presented data for 
a series of experiments on a slope of 1:3 (tan ¢ = 0.333). Bakhmeteff(1) and 
Matzke(6) performed experiments on slopes of 0 to 0.07 in a flume 6 inches 
wide. 


From an academic standpoint, the jump may occur in several ways ona 
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sloping apron, as outlined by Kinsdvater, presenting separate and distinct 
problems, Figure 30. Case A has the jump on a horizontal apron. In Case B, 
the toe of the jump forms on the slope, while the end occurs over the horizon- 
tal apron. In Case C, the toe of the jump is on the slope, and the end is at the 
junction of the slope and the horizontal apron; in Case D, the entire jump 
forms on the slope. With so many possibilities, it is easily understood why 
experimental data have been lacking on the sloping apron. Messrs. Yarnell, 
Kindsvater, Bakhmeteff, and Matzke limited their experiments to Case D. 

B. D. Rindlaub(7) of the University of California concentrated on the solution 


of Case B, but his experimental results are complete for only one slope, that 
of 12.330 (tan @ = 0.217). 


Sloping Apron Tests 


From a practical standpoint, the scope of the present test program need 
not be so broad as outlined in Figure 30. For example, the action in Cases C 
and D is for all practical purposes the same, if it is assumed that a horizon- 
tal floor begins at the end of the jump for Case D. Sufficient tests were made 
on Case C to verify the above statement that Cases C and D can be considered 
as one. The first of the current experiments to be described in this chapter 
involves Case D. The second set of tests deals with Case B. Case B is vir- 
tually Case A operating with excessive tail water depth. As the tail water 
depth is further increased, Case B approaches Case C. The results of Case 
A have already been discussed in the preceding chapters, and Cases D and B 
will be considered here in order. 


Tail Water Depth (Case D) 


Data obtained from the four flumes used in the sloping apron tests (Case D 
experiments) are tabulated in Table 8. The headings are much the same as 
those in previous tables, but need some explanation. Column 2 lists the tan- 
gents of the angles of the slopes tested. The depth of flow entering the jump, 
D1, Column 8, was measured at the beginning of the jump in each case, cor- 
responding to Section 1, Figure 30. It represents the average of a generous 
number of point gage measurements. The velocity at this same point, Vj, 
Column 7, was computed by dividing the unit discharge, q (Column 5), by Dj. 
The length of jump, Column 11, was measured in the flume, bearing in mind 
that the object of the test was to obtain practical data for stilling basin de- 
sign. The end of the jump was chosen as the point where the high velocity jet 
began to lift from the floor, or a point on the level tail water surface immedi- 
ately downstream from the surface roller, whichever occurred farthest down- 
stream. The length of the jump, as tabulated in Column 11, is the horizontal 
distance from Sections 1 to 2, Figure 30. The tail water depth, tabulated in 
Column 6, is the depth measured at the end of the jump, corresponding to the 
depth at Section 2 on Figure 30. 


The ratio x (Column 9, Table 8) is plotted with respect to the Froude 


number (Column 10) for sloping aprons having tangents 0.05 to 0.30 on Figure 
31. The plot for the horizontal apron (tan @ = 0) is the same as shown in 
Figure 5, (paper 1401). Superimposed on Figure 31 are data from 
Kindsvater,(5) Hickox,(5) Bahkmeteff,(1) and Matzke.(6) The agreement is 
within experimental error. 
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The small chart on Figure 31 was constructed using data from the larger 
chart, and shows, for a range of apron slopes, the ratio of tail water depth | 
for a continuous sloping apron, to conjugate depth for a horizontal apron. As 
indicated on the small sketch in Figure 31, Dg and TW are identical for a 
horizontal apron. The conjugate depth, Dg, listed in Column 14, Table 6, is 


the depth necessary for a jump to form on an imaginary horizontal floor be- 
ginning at Section 1, Figure 31. | 
The small chart, therefore, shows the extra depth, required for a jump of | 


a given Froude number to form on a sloping apron, rather than on a horizon- 
talapron. For example, if the tangent of the slope is 0.10, a tail water depth 
equal to 1.4 times the conjugate depth (D2 for a horizontal apron) will occur 
at the end of the jump; while if the slope is 0.30, the tail water depth at the f 
end of the jump will be 2.8 times the conjugate depth, Dg. The conjugate 
depth, Dg, used in connection with a sloping apron is merely a convenient 
reference figure which has no other meaning. It will be used throughout this 
discussion on sloping aprons. 


Length of Jump (Case D) 


The length of jump for the Case D experiments has been presented in two 
ways. First, the ratio length of jump to tail water depth, Column 12, was 
plotted with respect to the Froude number on Figure 32 for sloping aprons 
having tangents from 0 to 0.25. Secondly, the ratio of length of jump to the | : 
conjugate tail water depth, Column 16, Table 8, has been plotted with respect 
to the Froude number for the same range of slopes on Figure 33. Although not iS 
evident on Figure 32, it can be seen from Figure 33 that the length of jump on | 
a sloping apron is longer than that which occurs on a horizontal floor. For 
example, for a Froude number of 8, the ratio a varies from 6.1, for a hori- 
zontal apron, to 7.0, for an apron with a slope of 0.25. Length determinations 
from Kindsvater(5) for a slope of 0.167 are also plotted on Figure 32. The 
points show a wide spread. 


Expression for Jump on Sloping Apron (Case D) 


Several mathematicians and experimenters have developed expressions for 
the hydraulic jump on sloping aprons,(2,5,6,13) so there is no need to repeat 
any of these derivations here. An expression presented by Kindsvater(5) is 
the more common and perhaps the more practical to use. 


De 1 “cos? 5 


All symbols have been referred to previously, except for the coefficient K, a 
dimensionless parameter called the shape factor, which varies with the 
Froude number and the slope of the apron. Kindsvater and Hickox evaluated 
this coefficient from the profile of the jump and the measured floor pressures. 


| 
Surface profiles and pressures were not measured in the current tests but, as | 
a matter of interest, K was computed from Expression 5 by substituting ex- 

perimental values and solving for K. The resulting values of K are listed in I 


Column 17 of Table 8, and are shown plotted with respect to the Froude 
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number for the various slopes on Figure 34A. Superimposed on Figure 34A 
are data from Kindsvater for a slope of 0.167, and data from Hickox on a slope 
of 0.333. The agreement is not particularly striking nor do the points plot 
well, but it should be remembered that the value K is dependent on the method 
used for determining the length of jump. The current experiments indicate 
that the Froude number has little effect on the value of K. Assuming this to 
be the case, values of individual points for each slope were averaged and K is 
shown plotted with respect to tan g on Figure 34B. This phase is incidental 
to the study at hand and has been discussed only to complete the data analysis. 


Jump Characteristics (Case B) 


Case B is the one usually encountered in sloping apron design where the 
jump forms both on the slope and over the horizontal portion of the apron. Al- 
though this form of jump may appear quite complicated, it can be readily 
analyzed when approached from a practical standpoint. The primary concern 
in sloping apron design is the tail water depth required to move the front of 
the jump up the slope to Section 1, Figure 30B. There is little to be gained 
with a sloping apron unless the entire length of the sloping portion is utilized. 

Referring to the sketches on Figure 35A, it can be observed that for tail 
water equal to the conjugate depth, Dg, the front of the jump will occur ata 
point 0, a short distance up the slope. This distance is noted as lo and varies 
with the degree of slope. If the tail water depth is increased a vertical incre- 
ment, A Yj, it would be reasonable to assume that the front of the jump would 
raise a corresponding increment. This is not true; the jump profile undergoes 
an immediate change as the slope becomes part of the stilling basin. Thus, 
for an increase in tail water depth, AY}, the front of the jump moves up the 
Slope to Point 1, or moves a vertical distance AY; , which is several times 
AY. Increasing the tail water depth a second increment, say AY2, produces 
the same effect to a lesser degree, moving the front of the jump to Point 2. 
Additional increments of tail water depth produce the same effect but toa 
still lesser degree, and this continues until the tail water depth approaches 
1.3D9. For greater tail water depths, the relation is geometric; an increase 
in tail water depth, AY4, moves the front of the jump up the slope an equal 
vertical distance AY4, from Point 3 to 4. Should the slope be very flat, as in 
Figure 35B, the horizontal movement of the front of the jump is even more 
pronounced. 

The following studies were made to definitely tabulate the characteristics 
described above for conditions encountered in design since it has been neces- 
Sary in the past to check practically all sloping apron designs by model 
studies to be certain that the entire sloping portion of the apron was utilized. 


Experimental Results (Case B) 


The experiments for determining the magnitude of the profile characteris- 
tics were carried out on a large scale in Flume D, and the results are re- 
corded in Table 9. A sloping floor was placed in the flume as in Figure 30B. 
A discharge was established (Column 3, Table 9) and the depth of flow, Dy 
(Column 6), was measured immediately upstream from the front of the jump 
in each instance. The velocity entering the jump, Vj (Column 7), and the 
Froude number (Column 8) were computed. Entering Figure 31 with the com- 


puted values of Fy, the ratio D (Column 9) was obtained from the line labeled 
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“Horizontal apron.” Multiplying this ratio by Dy gives the conjugate depth for 
a horizontal apron which is listed in Column 10 of Table 9. The tail water 
was then set at conjugate depth (Point 0, Figure 35) and the distance, 1p), mea- 
sured and tabulated. The distance, lo, gives the position of the front of the 
jump on the slope, measured from the break in slope, for conjugate depth. 

The tail water was then increased, moving the front of the jump up to Point 1, 
Figure 35. Both the distance, 11, and the tail water depth were measured, and 
these are recorded in Columns 11 and 12, respectively, of Table 9. The tail 
water was then raised, moving the front of the jump to Point 2 while the 
length, 19, and the tail water depth were recorded. The same procedure was 
repeated until the entire apron was utilized by the jump. In each case, Dj 
was measured immediately upstream from the front of the jump, thus com- 
pensating for frictional resistance on the slope. The velocity, Vj, and the 
Froude number were computed at the same location. The tests were made 
for slopes with tangents varying from 0.05 to 0.30, and in some cases, several 
lengths of floor were used for each slope, as indicated in Column 15 of 

Table 9. 

The resulting lengths and tail water depths, divided by the conjugate depth, 
are shown in Columns 13 and 14 of Table 9, and these values have been 
plotted on Figure 36. The horizontal length has been used rather than the 
vertical distance, AY, as the former dimension is more convenient to use. 
Figure 36 shows that the straight lines for the geometric portion of the graph 
= and = 0.92, indicated by the 
circle on the graph. The change in the profile of the jump as it moves froma 
horizontal floor to the slope is evidenced by the curved portion of the lines. 

Case C, Figure 30, is the upper extreme of Case B; and as there is prac- 
tically no difference in the performance for Cases D and C, data for Case D 
(Table 4) can again be utilized. By assuming that a horizontal floor begins at 
the end of the jump in Case D, Columns 15 and 16 of Table 8 can be plotted 
on Figure 36. In addition, data from experiments by D. D. Hindlaub of the 
University of California, for a slope of 0.217, have been plotted on Figure 36. 
The agreement of the information from the three sources is very satisfactory. 


tend to intersect at a common point 


Length of Jump (Case B) 


It is suggested that the length of jump for Case B be obtained from Figure 
33. Actually, Figure 33 is for continuous sloping aprons, but these lengths 
can be applied to Case B with but negligible error. In some cases the length 
of jump is not of particular concern because it may not be economically pos- 
sible to design the basin to confine the entire jump. This is especially true 
when sloping aprons are used in conjunction with medium or high overfall 
spillways where the rock in the riverbed is in fairly good condition. When 
sloping aprons are designed shorter than the length indicated on Figure 33, 
the rock in the river downstream must act as part of the stilling basin. On 
the other hand, when the quality of foundation material is questionable, it is 
advisable to make the apron sufficiently long to confine the entire jump, 
Figure 33. 
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Existing Structures 


To determine the practical value of the methods given for the design of 
sloping aprons, existing basins employing sloping aprons were, in effect re- 
designed using the current experimental information. Pertinent data for 13 
existing spillways are tabulated in Table 10. The slope of the spillway face 
is listed in Column 3; the tangent of the sloping stilling basin apron is listed 
in Column 4; the elevation of the upstream end of the apron, or front of the 
jump, is listed in Column 7; the elevation of the end of the apron is listed in 
Column 8; the fall from headwater to upstream end of apron is tabulated in 
Column 9; and the total discharge is shown in Column 11. Where outlets dis- 
charge into the spillway stilling basin, that discharge was also been included 
in the total. The length of the sloping portion of the apron is given in Column 
14; the length of the horizontal portion of the apron is given in Column 15; 
and the overall length is given in Column 16. Columns 17 through 27 show 
computed values similar to those in the previous table. 

The lower portions of the curves of Figure 36 have been reproduced to a 
larger scale on Figure 37. The coordinates from Columns 26 and 27 of 
Table 10 have been plotted on Figure 37 for each of the 13 spillways. Longi- 
tudinal sections through the basins are shown on Figures 38 and 39. 

Each point in Figure 37 has been connected with an arrow to the tan¢d 
curve corresponding to the apron slope. Points which lie to the right and 
below the corresponding tan ¢ curve indicate that if the tail water depth is 
correct the sloping portion of the apron is excessively long; if the length of 
the slope is correct the tail water is insufficient to move the jump upstream 
to Section 1 on the slope. Only the points for Capilano and Madden Dams 
show an excess of tail water depth for the length of slope used. In both these 
cases the jump will occur upstream from Section 1 as shown in Figures 38 
and 39. Friant and Dickinson show almost perfect agreement with the derived 
curves while Bhakra (final) and Norris show agreement within practical 
limits. All other points indicate that the tail water depth is insufficient to 
move the toe of the jump upstream to Section 1. The rather large chute 
blocks on Keswick Dam may compensate for the discrepancy indicated by the 
point in the margin of Figure 37. 

All of the structures listed in Table 10 and shown on Figures 38 and 39 
were designed with the aid of model studies. The degree of conservatism 
used in each case was dependent on local conditions and the individual de- 
signer. 

The overall lengths of apron provided for the above 13 existing structures 
are shown in Column 16 of Table 10. The length of jump for the maximum 
discharge condition for each case is tabulated in Column 29 of the same table. 
The ratio of total length of apron to length of jump is shown in Column 30. 
The total apron length ranges from 39 to 83 percent of the length of jump; or 
considering the 13 structures collectively, the average total length of apron 
is 60 percent of the length of the jump. Considering all aspects of the model 
tests on the individual structures and the current basin tests it is believed 
that 60 percent is sufficient for most installations. Longer basins are needed 
only when the downstream riverbed is in very poor condition. Shorter basins 
may be used where a solid bed exists. 
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Evaluation of Sloping Aprons 


Many sloping aprons have been designed so that the jump height curve 
matches the tail water curve for all discharge conditions. This procedure 
results in what has been designated a “tailor-made” basin. Some of the exist- 
ing basins shown on Figures 38 and 39 were designed in this manner. Asa 
result of the current experiments it was discovered that this course is not the 
most desirable. Matching of the jump height curve with the tail water curve 
should be a secondary consideration, except for the maximum discharge con- 
dition. 

The first consideration in design should be to determine the apron slope 
that will require the minimum amount of excavation, the minimum amount of 
concrete, or both, for the maximum discharge and tail water condition. This 
is the prime consideration. Only then is the jump height checked to determine 
whether the tail water depth is adequate for the intermediate discharges. It 
will be found that the tail water depth usually exceeds the required jump 
height for the intermediate discharges resulting in a slightly submerged con- 
dition for intermediate discharges, but performance will be very acceptable. 
The extra depth will provide a smoother water surface in and downstream 
from the basin and greater stability at the toe of the jump. Should the tail 
water depth be insufficient for intermediate flows, it will be necessary to 
increase the depth by increasing the slope, or reverting to a horizontal apron. 
It is not necessary that the front of the jump form at the upstream end of the 
sloping apron for low or intermediate discharges provided the tail water depth 
and the length of basin available for energy dissipation are considered ade- 
quate. With this method, the designer is free to chose the slope he desires, 
since the current tests showed, beyond a doubt, that the slope itself had little 
effect on the performance of the stilling basin. 

It is not possible to standardize design procedures for sloping aprons to 
the degree shown for the horizontal aprons; greater individual judgment is re- 
quired. The slope and overall shape of the apron must be determined from 
economic reasoning, while the length must be judged by the type and sound- 
ness of the riverbed downstream. The existing structures shown on Figures 


38 and 39 should serve as a guide in proportioning future sloping apron 
designs. 


Sloping Apron Versus Horizontal Apron 


The Bureau of Reclamation has constructed very few stilling basins with 
horizontal aprons for its larger dams. It has been the consensus that the hy- 
draulic jump on a horizontal apron is very sensitive to slight changes in tail 
water depth. The current tests demonstrate this to be true for the larger 
values of the Froude number, but this characteristic can be remedied. Sup- 
pose a horizontal apron is designed for a Froude number of 10. The basin 
will operate satisfactorily for conjugate tail water depth, but as the tail water 
is lowered to 0.98D9 the front of the jump will begin to move. By the time the 
tail water is dropped to 0.96D9, the jump will probably be completely out of 
the basin. Thus, to design a stilling basin in this range the tail water depth 
must be known with certainty or a factor of safety provided in the design. To 
guard against deficiency in tail water depth, the same procedure used for 
Basins I and II (papers 1401 and 1402) is suggested here. Referring to the 
minimum tail water curve for Basins I and II on Figure 11, (paper 1402) the 
margin of safety can be observed for any value of the Froude number. It is 
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recommended that the tail water depth for maximum discharge be at least 5 


Froude number greater than 9, a 10 percent factor of safety may be advisable 
as this will not only stabilize the jump but will improve the performance of 


the basin. With the additional tail water depth, the horizontal apron will per- 


form on a par with the sloping apron. Thus, the primary consideration in de- 
sign need not be hydraulic but structural. The basin, with either horizontal 
or sloping apron, which can be constructed at the least cost is the most 
desirable. 


Effect of Slope of Chute 


A factor which occasionally affects stilling basin operation is the slope of 
the chute upstream from the basin. The foregoing experimentation was suffi- 
ciently extensive to shed some light on this factor. The tests showed that the 
slope of chute upstream from the stilling basin was unimportant, as far as 
jump performance was concerned, so long as the velocity distribution in the 
jet entering the jump was reasonably uniform. In the case of steep chutes or 
short flat chutes, the velocity distribution can be considered normal. Diffi- 
culty is experienced, however, with long flat chutes where frictional resist- 
ance on the bottom and side walls is sufficient to produce a center velocity 
greatly exceeding that on the bottom or sides. When this occurs, greater 
activity results in the center of the stilling basin than at the sides, producing 
an asymmetrical jump with strong side eddies. This same effect is also wit- 
nessed when the angle of divergence of a chute is too great for the water to 
follow properly. In either case the surface of the jump is unusually rough and 
choppy and the position of the front of the jump is not always predictable. 

When long chutes precede a stilling basin the practice has been to make 
the upstream portion unusually flat, then increase the slope to 2:1, or that 
corresponding to the natural trajectory of the jet, immediately preceding the 
stilling basin. Figure 1A, (paper 1401) which shows the model spillway for 
Trenton Dam, illustrates this practice. Bringing an asymmetrical jet into the 
stilling basin at a steep angle usually helps to redistribute the flow to stabilize 
the jump. This is not effective, however, where very long flat slopes have 
caused the velocity distribution to be completely out of balance. 

The most adverse condition has been observed where long canal chutes 
terminate in stilling basins. A typical example is the chute and basin at Sta- 
tion 25419 on the South Canal, Uncompahgre Project, Colorado, Figure 40. 
The operation of this stilling basin is not particularly objectionable, but it 
will serve as an illustration. The above chute is approximately 700 feet long 
with a slope of 0.0392. The stilling basin at the end is also shown on 
Figure 40. A photograph of the prototype basin operating at normal capacity 
is shown on Figure 41. The action is of the surging type; the jump is unusual- 
ly rough, with a great amount of splash and spray. Two factors contribute to 
the rough operation: the unbalanced velocity distribution in the entering jet, 
and excessive divergence of the chute in the steepest portion. 

A definite improvement can be accomplished in future designs where long 
flat chutes are involved by utilizing the Type III basin described in Paper 1403. 
The baffle piers on the floor tend to alter the asymmetrical jet, resulting in 
an overall improvement in operation. 
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FIGURE 40 
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Recommendations 


The following rules have been devised for the design of the sloping aprons 
developed from the foregoing experiments: 


1. Determine an apron arrangement which will give the greatest economy 
for the maximum discharge condition. This is the governing factor and the 
only justification for using a sloping apron. 


2. Position the apron so that the front of the jump will form at the up- 
stream end of the slope for the maximum discharge and tail water condition 
by means of the information on Figure 37. Several trials will usually be re- 
quired before the slope and location of the apron are compatible with the hy- 
draulic requirement. It may be necessary to raise or lower the apron, or 
change the original slope entirely. 


3. The length of the jump for maximum or partial flows can be obtained 
from Figure 33. The portion of the jump to be confined on the stilling basin 
apron is a decision for the designer. In making this decision, Figures 38 and 
39 may be helpful. The average overall apron in Figures 38 and 39 averages 
60 percent of the length of jump for the maximum discharge condition. The 
apron may be lengthened or shortened, depending upon the quality of the rock 
in the riverbed and other local conditions. If the apron is set on loose mate- 
rial and the downstream channel is in poor condition, it may be advisable to 
make the total length of apron the same as the length of jump. 


4. With the apron designed properly for the maximum discharge condition, 
it should then be determined that the tail water depth and length of basin 
available for energy dissipation are sufficient for, say, 1/4, 1/2, and 3/4 
capacity. If the tail water depth is sufficient or in excess of the jump height 
for the intermediate discharges, the design is acceptable. If the tail water 
depth is deficient, it may then be necessary to try a flatter slope or reposi- 
tion the sloping portion of the apron. It is not necessary that the front of the 
jump form at the upstream end of the sloping apron for partial flows. In 
other words, the front of the jump may remain at Section 1 (Figure 30B), move 
upstream from Section 1, or move down the slope for partial flows, providing 


the tail water depth and length of apron are considered sufficient for these 
flows. 


5. Horizontal and sloping aprons will perform equally well for high values 
of the Froude number if the proper tail water depth is provided. 


6. The slope of the chute upstream from a stilling basin has little effect on 
the hydraulic jump when the velocity distribution and depth of flow are rea- 
sonably uniform on entering the jump. 


7. A small solid triangular sill, placed at the end of the apron, is the only 
appurtenance needed in conjunction with the sloping apron. It serves to lift 
the flow as it leaves the apron and thus acts to control scour. Its dimensions 
are not critical; the most effective height is between 0.05D9 and 0.10D9 and a 
slope of 3:1 to 2:1 (see Figures 38 and 39). 


8. The spillway should be designed to operate with as nearly symmetrical 
flow in the stilling basin as possible. (This applies to all stilling basins.) 
Asymmetry produces large horizontal eddies that can carry riverbed material 
onto the apron. This material, circulated by the eddies, can abrade the apron 


| 
| ASCE 1405-31 

| 

| 

| 
| 

eee 


1405-32 HY 5 October, 1957 


and appurtenances in the basin at a very surprising rate. Eddies can also 
undermine wing walls and riprap. Asymmetrical operation is expensive oper- 
ation, and operating personnel should be continually reminded of this fact. 


9. Where the discharge over high spillways exceeds 500 cfs per foot of 
apron width, where there is any form of asymmetry involved, and for the 
higher values of the Froude number where stilling basins become increasing- 
ly costly and the performance less acceptable a model study is advisable. 
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HYDRAULIC DESIGN OF STILLING BASINS: SMALL BASINS FOR PIPE 
OR OPEN CHANNEL OUTLETS—NO TAIL WATER REQUIRED (BASIN VI) 
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FOREWORD 


This is one of a group of six papers on the hydraulic design of stilling ba- 
sins and their associated appurtenances. Although original data from hydrau- 
lic models are presented in tabular form, the papers emphasize practical de- 
sign procedures. Basin sizes and dimensions are given in easy-to-use 
dimensionless forms so that dependable stilling basins may be designed with- 
out the need for exceptional judgment or extensive experience. Sample prob- 
lems are included. 

The “Introduction” and “Experimental Equipment” are included in Paper 
1401 but apply to all papers in this group. Paper 1401 is an academic study 
of the hydraulic jump on a flat floor. Useful information is presented con- 
cerning energy losses, applicability of the hydraulic jump formula, length of 
jump, and a new classification of the various types of jumps. 

Paper 1402 covers the design of a short hydraulic jump stilling basin hav- 
ing an end sill and chute blocks. Paper 1403 describes a shorter stilling ba- 
Sin which utilizes baffle piers. In both papers information is given to deter- 
mine the critical dimensions of the stilling basin, the tail water range, and 
the water surface profile in the basin. 

Paper 1404 describes a special type of hydraulic jump basin for use when 
the Froude number of the incoming flow is low, and the jump produces waves 
in the downstream channel. An alternate design and two types of wave sup- 
pressors are also developed. 

Paper 1405 describes the design of a stilling basin having a sloping apron. 
The extra tail water depth required and the merits of a sloping apron are 
evaluated, as is the basin length for a range of apron slopes. 

Paper 1406 develops the design of impact basins for discharges up to 340 


Note: Discussion open until March 1, 1958. Paper 1406 is part of the copyrighted 
Journal of the Hydraulics Division of the American Society of Civil Engineers, 
Vol. 83, No. HY 5, October, 1957. 


1. Hydr. Engr., U. S. Bureau of Public Roads, Washington, D. C., formerly of 
U. S. Bureau of Reclamation, Denver, Colo. 
. Hydr. Engr., U. S. Bureau of Reclamation, Denver, Colo. 
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second feet and incoming velocities up to 30 feet per second. No tail water is 
required. The bibliography in this paper applies to all papers in this group. 


ABSTRACT 


An impact-type stilling basin is developed which will handle 21 to 340 
cubic feet per second. No tail water is required and incoming velocities may 
be as great as 30 feet per second. General design rules and a table listing 


critical dimensions and proper riprap sizes for 9 different basin sizes are 
presented. 


INTRODUCTION 


The stilling basin developed in these tests is an impact-type energy dis- 
sipator, contained in a relatively small boxlike structure, which requires no 
tail water for successful performance. Although the emphasis in this discus- 
sion is placed on use with pipe outlets, the entrance structure may be modi- 
fied for use with an open channel entrance. 

Generalized design rules and procedures are presented to allow determin- 
ing the proper basin size and all critical dimensions for a range of discharges 
up to 339 feet per second and velocities up to about 30 feet per second. 
Greater discharges may be handled by constructing multiple units side by 
side. The efficiency of the basin in accomplishing energy losses is greater 
than a hydraulic jump of the same Froude number. 

The development of this short impact-type basin was initiated by the need 
for some 50 or more stilling structures on a single irrigation project. The 
need was for relatively small basins providing energy dissipation independent 
of a tail water curve or tail water of any kind. 

Since individual model studies on 50 small stilling structures were too 
costly a procedure, tests were made on a single setup which was modified as 
necessary to generalize the design for the range of expected operations. 


Test Procedure 


Hydraulic Models 


Hydraulic models were used to develop the stilling basin, determine the 
discharge limitations, and obtain dimensions for the various parts of the ba- 
sin. Basins 1.6 to 2.0 feet wide were used in the tests. The inlet pipe was 
6 3/8 inches, inside diameter, and was equipped with a slide gate located 
well upstream from the basin entrance so that the desired relations between 
head, depth, and velocity could be obtained. The pipe was transparent so that 
backwater effects in the pipe could be studied. Discharges of over 3 cubic 
feet per second and velocities up to 15 feet per second could be obtained dur- 
ing the tests. Hydraulic model-prototype relations were used to scale up the 
results to predict performance for discharges up to 339 second-feet and 
velocities up to 30 feet per second. 


The basin was tested in a tail box containing gravel formed into a 
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trapezoidal channel. The size of the gravel was changed several times during 
the tests. The outlet channel bottom was slightly wider than the basin and 

had 1:1 side slopes. A tail gate was provided at the downstream end to evalu- 
ate the effects of tail water. 


Development of Basin 


The shape of the basin evolved from the development tests was the result 
of extensive investigations on many different arrangements. These tests are 
discussed briefly to show the need for the various parts of the adopted design. 

With the many combinations of discharge, velocity, and depth possible for 
the incoming flow, it became apparent during the early tests that some device 
was needed at the stilling basin entrance to convert the many possible flow 
patterns into a common pattern. The vertical hanging baffle proved to be this 
device, Figure 42. Regardless of the depth or velocity of the incoming flow 
(within the prescribed limits) the flow after striking the baffle acted the same 
as any other combination of depth and velocity. Thus, some of the variables 
were eliminated from the problem. 

The effect of velocity alone was then investigated, and it was found that for 
velocities 30 feet per second and below the performance of the structure was 
primarily dependent on the discharge. Actually, the velocity of the incoming 
flow does affect the performance of the basin, but from a practical point of 
view it could be eliminated from consideration. Had this not been done, a 
prohibitive amount of testing would have been required to evaluate and ex- 
press the effect of velocity. 

For velocities of 30 feet per second or less the basin width W was found to 
be a function of the discharge, with other basin dimensions being related to 
the width, Figure 42. To determine the necessary width, erosion test results, 
judgment, and operating experiences were all used, and the advice of labora- 
tory and design personnel was used to obtain the finally determined limits. 
Since no definite line of demarcation between a “too wide” or “too narrow” 
basin exists, it was necessary to work between two more definite lines, 
shown on Figure 42 as the upper and lower limits. These lines required far 
less judgment to determine than a single intermediate line. 

Various basin sizes, discharges, and velocities were tested taking note of 
the erosion, wave heights, energy losses, and general performance. When the 
upper and lower limit lines had been established a line about midway between 


the two was used to establish the proper width of basin for various discharges. 


The exact line is not shown because strict adherence to a single curve would 
result in difficult to use fractional dimensions. Accuracy of this degree is 
not justifiable. Figure 43 shows typical performance of the recommended 
stilling basin for three limits. It is evident that the center photograph rep- 
resents a compromise between the upper limit operation which is very mild 
and the lower limit operation which is approaching the unsafe range. 

Using the middle range of basin widths, other basin dimensions were de- 
termined, modified, and made minimum by means of trial and error tests on 
the several models. Dimensions for nine different basins are shown in Table 
11. These should not be arbitrarily reduced since in the interests of economy 
the dimensions have been reduced as far as is safely possible. 


Performance of Basin 


Energy dissipation is initiated by flow striking the vertical hanging baffle 
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and being turned upstream by the horizontal portion of the baffle and by the 
floor, in vertical eddies. The structure, therefore, requires no tail water for 
energy dissipation as is necessary for a hydraulic jump basin. Tail water as 


high as d + 4 Figure 42, however, will improve the performance by reducing 


outlet velocities, providing a smooth water surface, and reducing tendencies 
toward erosion. Excessive tail water, on the other hand, will cause some flow 
to pass over the top of the baffle. This should be avoided if possible. 

The effectiveness of the basin is best illustrated by comparing the energy 
losses within the structure to those which occur in a hydraulic jump. Based 
on depth and velocity measurements made in the approach pipe and in the 
downstream channel (no tail water), the change in momentum was computed 
as explained in Paper 1401 for the hydraulic jump. The Froude number of the 
incoming flow was computed using the D1 obtained by converting the flow area 
in the partly full pipe into an equivalent rectangle as wide as the pipe diame- 
ter. Compared to the losses in the hydraulic jump, Figure 44, the impact 
basin shows greater efficiency in performance. Inasmuch as the basin would 
have performed just as efficiently had the flow been introduced in a rectangu- 
lar cross section, the above conclusion is valid. 


Basin Design 


Table 11 and the key drawing, Figure 42, may be used to obtain dimensions 
for the usual structure operating within usual ranges. However, a further 
understanding of the design limitations may help the designer to modify these 
dimensions when necessary for special operating conditions. 

The basin dimensions, Columns 4 to 13, are a function of the maximum dis- 
charge to be expected, Column 3. Velocity at the stilling basin entrance need 
not be considered except that it should not greatly exceed 30 feet per second. 

Columns 1 and 2 give the pipe sizes which have been used in field installa- 
tions. These may be changed as necessary, however. The suggested sizes 
were obtained by assuming the velocity of flow to be 12 feet per second. The 
pipes shown would then flow full at maximum discharge or they would flow 
half full at 24 feet per second. The basin operates as well whether a small 
pipe flowing full or a larger pipe flowing partially full is used. The pipe size 
may therefore be modified to fit existing conditions, but the relation between 
structure size and discharge should be maintained as given in the table. In 
fact, a pipe need not be used at all; an open channel having a width less than 
the basin width will perform equally as well. 

The invert of the entrance pipe, or open channel, should be held at the ele- 
vation shown on the drawing of Figure 42, in line with the bottom of the baffle 
and the top of the end sill, regardless of the size of the pipe selected. The 
entrance pipe may be tilted downward somewhat without affecting performance 
adversely. A limit of 15° is a suggested maximum although the loss in effi- 
ciency at 20° may not cause excessive erosion. For greater slopes use a hori- 
zontal or sloping pipe (up to 15°) 2 or more diameters long just upstream 
from the stilling basin. 

For submerged conditions a hydraulic jump may be expected to form in 
the downstream end of the pipe, sealing the exit end. If the upper end of the 
pipe is also sealed by incoming flow, a vent may be necessary to prevent 
pressure fluctuation in the system. A vent to the atmosphere, say one-sixth 
the pipe diameter, should be installed upstream from the jump. 
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The notches shown in the baffle are provided to aid in cleaning out the ba- 
sin after prolonged nonuse of the structure. When the basin has silted level 
full of sediment before the start of the spill, the notches provide concentrated 
jets of water to clean the basin. If cleaning action is not considered neces- 
sary the notches need not be constructed. The basin is designed, however, to 
carry the full discharge, shown in Table 11, over the top of the baffle if for 
any reason the space beneath the baffle becomes clogged, Figure 45C. Per- 
formance is not as good, naturally, but acceptable. With the basin operating 
normally, the notches provide some concentration of flow passing over the end 
sill, resulting in some tendency to scour, Figure 45A. Riprap as shown on 
the drawing will provide ample protection in the usual installation, but if the 
best possible performance is desired, it is recommended that the alternate 
end sill and 45° end walls be used, Figures 45B and 42. The extra sill length 
reduces flow concentration, scour tendencies, and the height of waves in the 
downstream channel. 

Figure 46 shows the performance of a prototype structure designed from 
Table 11. The basin, designed for a maximum discharge of 165 second-feet, 
is shown discharging 130 second-feet with a higher than recommended en- 
trance velocity of about 39 feet per second. Performance is entirely 
satisfactory. 


CONCLUSIONS AND RECOMMENDATIONS 


The following procedures and rules pertain to the design of Basin VI: 


1. Use of Basin VI is limited to cases where the velocity at the entrance 
to the stilling basin does not greatly exceed 30 feet per second. 


2. From the maximum expected discharge, determine the stilling basin 
dimensions, using Table 11, Columns 3 to 13. The use of multiple units side 
by side may prove economical in some cases. 


3. Compute the necessary pipe area from the velocity and discharge. The 
values in Table 11, Columns 1 and 2, are suggested sizes based on a velocity 
of 12 feet per second and the desire that the pipe run full at the discharge 
given in Column 3. Regardless of the pipe size chosen, maintain the relation 
between discharge and basin size given in the table. An open channel entrance 
may be used in place of a pipe. The approach channel should be narrower 
than the basin with invert elevation the same as the pipe. 


4. Although tail water is not necessary for successful operation, a moder- 
ate depth of tail water will improve the performance. For best performance 
set the basin so that maximum tail water does not exceed d + é, Figure 42. 

5. Suggested thicknesses of various parts of the basin are given in 
Columns 14 to 18, Table 11. 


6. The suggested sizes for the riprap protective blanket, given in Column 
19 of Table 11, show the minimum size of individual stones which will resist 
movement when critical velocity occurs over the end sill. Since little is 
known regarding the effect of interlocking rock pieces, most of the riprap 
should consist of the sizes given or larger. An equation(34,35) for determin- 
ing minimum stone sizes, which appears from a limited number of experi- 
ments and observations to be accurate, is given below. 
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Vp = bottom velocity in feet per second 
d = diameter of rock in inches 


The rock is assumed to be of the ordinary variety having a specific gravity of 
about 2.65. The accuracy of the equation is not known for velocities above 16 
feet per second. 


7. The entrance pipe or channel may be tilted downward about 15° without 
affecting performance adversely. For greater slopes use a horizontal or 
sloping pipe (up to 15°) 2 or more diameters long just upstream from the 
stilling basin. Maintain proper elevation of invert at entrance as shown on 
the drawing. 


8. If a hydraulic jump is expected to form in the downstream end of the 
pipe and the pipe entrance is sealed by incoming flow, install a vent about one- 
sixth the pipe diameter at any convenient location upstream from the jump. 


9. For best possible operation of basin use alternate end sill and 45° wall 
design shown on Figure 42. Erosion tendencies will be reduced as shown on 
Figure 45. 
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FIGURE 42 
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Figure 46 
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Discussion of 
“FLOOD PROTECTION OF CANALS BY LATERAL SPILLWAYS” 


by Harald Tults 
(Proc. Paper 1077) 


HARALD TULTS,! A.M. ASCE.—The discussion of this paper by Professor 
Ven Te Chow describes the application of the momentum and energy princi- 
ples to the flow with increasing and decreasing discharges. 

Professor Chow states that the water level variation in spacially varied 
flow may be computed either by the energy or momentum principles. 

The energy principle is generally applicable when the change of the energy 
in the flow is computable. That is the case with the flow along the lateral 
spillway crest, which is subject to friction loss only. However, the kinetic 
energy of spilled water is almost entirely lost at the drop into a side channel 
as demonstrated by the laboratory tests. Thus, the statement by Professor 
Chow that “the energy due to the added discharge dQ per elementary distance 
dx should be added to the total energy along the length of the channel during 
the time interval dt” is not substantiated. The flow in the side channels is in- 
vestigated by Messrs. W. H. R. Nimmo,!1 H. Favre,17 J. Hinds,!2 and T. R. 
Camp,18 but none of them attempts to compute the water level in the side 
channels by the Bernouilli energy equation. 

Unfortunately, Professor Chow does not show how his Eq. (D-1) is derived 
and how it can be transformed to Eq. (D-2). A superficial analysis indicates, 
that Eq. (D-1) may be applicable only for rectangular canals because the term 
“D” is defined as hydraulic depth equal to the cross sectional area, A, divided 
by its top width. 

Furthermore, Eq. (D-3) for lateral spillway flow represents only the super- 
critical flow (Cases d and e in Fig. 1) as shown by the discussion of the dif- 
ferential equation (Eq. 18) for lateral spillway flow. This differential equa- 
tion for the lateral spillway was omitted in the writers paper because the 
application of the (y,Q)-curves is much simpler in practical design. However, 
in regard to the discussion by Professor Chow, it wouic be expedient to in- 
clude the development of the differential equation of lateral spillway flow into 
the closure of the paper. 

For the differentiation of the Bernouilli equation, the friction term, Sp, as 
a function of the velocity, of the cross section, and of the wetted perimeter is 
eliminated by equating it to the bottom slope, S. 

The first derivative of the simplified version of Eq. 10, 


H 


is 
Ss 


v 2gke ’ 
1. Hydr. Engr., Bechtel Corp., Vernon, Calif. 
17. “Contributions a 1’Etude des Eaux Currents,” Editor Rascher & Cie., 

Zurich, Leipzig and Stuttgart, 1933. 


18. “Lateral Spillway Channels,” ASCE Transactions, Vol. No. 105, 1940, 
pp. 636-637. 
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dA d 15 
* (15) 


is zero because of the constancy of the energy head, Hg, of lateral spillway 
flow. 


This general equation applied to a rectangular uniform lateral spillway 


canal is transformed by replacing a“ by ve to 


= ga? 


in which the spillway discharge per increment of length, Cis e 


3/ 
= - 2)? (17) 


As the value of aa is always negative, the water level variation, , may 
be positive or negative depending on the sign of the factor CyQb - gA3. 

When the value of CyQ2b - gA3 is negative, the water level along the spill- 
way canal is rising in the downstream direction as in the subcritical flow, 
shown in Fig. 1 -a, b, 

In case, when C, yo? b- gA3 is positive, the water level along the spillway 
crest is subsiding as in the case of supercritical flow, Fig. 1 -dande. In 
this case, Eq. 16 is fully identical with Eq. (D-3) by Professor Chow, when S 
is made equal to So. 


For the integration, Eq. 16 is transformed using Eqs. 11 and 17 into the 
following form: 


2c y) 3/ 
= -2) ° 
b \2g 2H, - 39 


(18) 


Integration of Eq. 18 between the limits x and xg (Xo - the abscissa at which y 
is equal to Hg, i.e. Q is equal to zero according to Eq. 11) results in the fol- 
lowing equation as shown by Professor A. Schocklitsch13 and modified by 


blz 2H |#,-¥ 
x= 


Ho - Zz 
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The required length of the lateral spillway, replacing 


s ifs 
b ly 
C ple (20) 


It is obvious from the mathematical work required that the application of 
Eq. 11 for the computation of the dimensions of a lateral spillway is much 
simpler than the integration of the fundamental differential Eq. 18 derived for 
the lateral spillway flow in canals of rectangular cross-section. Besides that, 
integration of Eq. 15, adapted to a more complex spillway cross section than 
to the rectangular one, is very complicated, perhaps even impossible. How- 
ever, it is possible for any kind of canal cross-section to establish an equa- 
tion Q = f(A,Hs,y) of the same type as Eq. 11. 

The writer is obliged to Professor Ven Te Chow for corrections, particu- 
larly in connection with the derivation of Eq. 6. Eq. 6 itself is correct, but 
the sign error during typing of Eq. 3 caused the writer to change “sub- 
tracted” to “added” at last check. 

Actually, the writer expected more discussions, particularly field observa- 
tions about the performances of lateral spillways and about the priming of the 
spillway at the filling of the canal. 

It is hoped that the writers paper and closure will contribute to the under- 


Standing of lateral spillway flow and of its application for flood protection of 
open canals. 
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Discussion of 
“GRAPHICAL DETERMINATION OF WATER-SURFACE PROFILES” 


by Francis F. Escoffier 
(Proc. Paper 1114) 


FRANCIS F. ESCOFFIER, | A.M. ASCE.—Mr. Ven Te Chow has pointed out 
how the accuracy of the graphical method is affected by the variation of the 
critical slope and proposes a different graphical construction in which he 

plots as the abscissa the product wB rather than the quantity B. The differen- 


tial equation which serves as the starting point for both methods can be 
written. 


Sodx=dy yowdB (22) 


When this is integrated from x, to xg there is obtained 


The writer followed Bakhmeteff in assuming, as an approximation, that w re- 
mains constant. The solution proposed by Mr. Chow replaces eq. (23) by 
So (2-1) -Y) “Yo (WB)z + Yo (WB) 


To compare eq.’s (23) and (24) we consider the identity 


By a 


~ 


which can be recognized as a variation of the formula for integration by parts. 
On appropriate substitution from eq. (25) into eq. (24) there is obtained 


So («2 


wdB-y, [2 Bdw (26) 


B, 


On comparing eq. (26) with eq. (23) it is seen that eq. (26) contains the 


extraneous term 
w 
2 Bdw 
Wy 


1. Hydr. Engr., Corps of Engrs., U. S. Dept. of the Army, Mobile Dist., 


Mobile, Ala. 
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and that therefore eq.’s (24) and (26) are in error to that extent. 

Mr. Chow has outlined a graphical method in which he uses the varied- 
flow functions published by him. Another method similar to that shown in 
fig. 7 can also be used with these functions. In this method the abscissa is 


1 B, the ordinate in the lower diagram is B - 7) and that in the upper diagram 


is 7. The slope w is 


w $e) 


The two methods are mathematically equivalent and will therefore give the 
same numerical results. However, Mr. Chow’s method requires a replotting 


of the curves for each change in discharge whereas that of the writer does not. 


Mr. Chow’s tables of varied-flow functions for adverse slopes represent a 
distinct advance over Matzke’s table not only because of the wider range of 
exponents that are included but also because of the wider range in the depth 
variable 7 and because of the smaller intervals in both the exponent and the 
depth variable. 

Mr. Chow refers to a chart in his own papers similar to the chart in fig. 8. 
His chart (fig. 2 in his paper) includes a larger number of side slopes for 
trapezoidal channels than fig. 8 and can be read more accurately. Attention 
should also be directed at this point to figs. 3, 4, and 5 in his paper. Fig. 3 
is similar to fig. 2 but gives the exponent m instead of the exponent n. Fig. 4 
is useful in determining the normal depth and fig. 5, in determining the criti- 
cal depth. These four charts will prove invaluable to anyone engaged in com- 
puting water-surface profiles in trapezoidal or circular channels. 

The writer wishes to thank Mr. Chow for his most interesting discussion. 
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Discussion of 
“MECHANICS OF SEDIMENT-RIPPLE FORMATION” 


by H. K. Liu 
(Proc. Paper 1197) 


THOMAS MADDOCK, JR.,1 M. ASCE.—Dr. Liu has presented an important 
addition to our slowly growing understanding of the sediment movement prob- 
lem. From Dr. Liu’s work one would conclude that ripple formation is 
limited to those conditions during which the flow around the sediment grains 
is neither laminar or viscous nor turbulent, but is in that area which Horton 
called “mixed.” 

In conventional hydraulics “mixed” flow is relatively unimportant and has 
been the subject of little study. It appears likely that an understanding of this 
type of flow is a prerequisite to a resolution of sediment problems since it 
may be that ripple forms are simply the visual evidence of unstable flow con- 
ditions around the grains. The lengthening of the ripple form as the flow 
changes from laminar to turbulent in terms of the Reynolds number Ud/v as 
noted by Shields,(26) might be explained in an oversimplified fashion, as being 
due to changes in the cycle of formation of vortices in the boundary layer. 

Thus, we are not interested in the values of a Reynolds number at which 
flow around the grains becomes fully turbulent nor fully laminar but in learn- 
ing what happens in between these stages. It is not clear to the writer that 
those authorities cited by the author fully discuss what happens during the 
transition stage. 

One of the most interesting presentations by Dr. Liu is his fig. 10 showing 
the relation between his values for the beginning of ripples and the Shields 
curve for the beginning of bed-load movement. The two curves appear to coin- 
cide in the turbulent range but should not coincide in the laminar range. 

There is a considerable amount of sediment in movement at the time ripples 
form—enough to affect the hydraulics of the channel. Nevertheless, it might 
be assumed that ripples are found as soon as enough sediment becomes avail- 
able toform them. Thus, it is easy to presume that the flow conditions re- 
lated to the ripple existed at some smaller value of Ups d Yi v but, lacking ade- 
quate movement of sediment, there was no way to distinguish them. However, 
ripple formation might also be related to velocity distributions in the cross 
section of the channel. 

The shape of the curves on figs. 9 and 10 bears a considerable similarity 
to the usual Drag Coefficient (Cp) - Reynolds Number plot of particles falling 
in water. It is to be noted that U,/w becomes constant at high Reynolds num- 
bers. So does the drag coefficient, hence, it may be possible to presume that 


Us /w 0c Co . Now, if the Liu and er” curves on fig. 10 are extended 
in the laminar range, we find that /wee or, 


of the equation are functions of Vv, ot w and awh 
or Ur /ww oc Vuh Ate At small Reynolds numbers wed? so 
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Ur oc But inthis size Cpeoe Soit canbe 
concluded that Ux Jw 
Since from equation (8) and since Cz- Cp, 


must also be proportional to Cp in spite of the way in which it is derived. 
More importantly, however, this expression indicates a strong relationship 
between U, and d. 


where wW is generally called the flow intensity. 
The data from the Liu and Shields curves indicate, therefore, that 


Ux /gd or /y- is roughly constant over a wide range of particle sizes at 
the time sediment movement begins or ceases. It should be noted that the 


ratios of Ox/ ah? found in Table I are 1.23, 1.14, 1.07, 1.10, 0.93, and 
.975. 


This relationship is confirmed to some extent by inspection of the curves 
of the Einstein bed-load function and the Meyer-Peter formula. Both show 
constant values of w for low values of ¢, the intensity of bed-load transport. 
Furthermore, fig. 11 shows that there is no direct relationship between 


Whe -r)d and Ur d/y . While there is considerable scatter, there is 
no regular variation apparent in tH cs -rd . This latter value is equivalent 
to RS a or Vy 
There is enough scatter in fig. 11, however, to indicate that R5 (3:6) 


cannot be assumed to be precisely constant at the time movement begins. 
Part of the scatter may be caused by variation in channel shape or the un- 
recognized effect of particle shape on settling velocity. Certainly sediment 
load is responsible for part. There are also indications that another parame- 


ter may be involved since Shields(26) concluded that irs -r)d was not 
a constant at the beginning of sediment movement. The writer believes that 


Dr. Liu has answered his own question and that Te/(-1)d. is not a 
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satisfactory parameter to define the beginning of sediment movement. It is 
too constant to be used as a variable and too variable to be used as a constant. 
The evidence presented by Dr. Liu indicates that Ux/us is a better 


parameter when related to Us d/v. The discussion here indicates that 


(232) ad could be used as a substitute. However, it would be effective 


only if precise information were available to describe conditions at the begin- 
ning of the movement. 

It could be assumed that if w is reasonably constant at the beginning of 
movement of a wide range of sediment sizes, it would be a satisfactory 
parameter to describe the rate of sediment movement. As a matter of fact, 
Ning Chien (1956) showed that 


“ 0/88) 


which is in accord with such an assumption. 

However, the same Journal carrying Dr. Liu’s paper contains Brook’s 
closure of his very fine paper, “Mechanics of Streams with Movable Beds of 
Fine Sand.” In it he shows a poor relation between y and sediment movement 
and concludes “Because of enormous changes in roughness have been ob- 
served to occur with almost negligible changes in bed shear, it appears that 
the latter is a poor choice of independent variable.” Thus, there is strong 


evidence to indicate that R5/ 42a is a poor independent parameter to 


be used to describe either the beginning or the rate of sediment movement. 

It is worthy of note that Durand (1953) who has presented the only really 
definitive description of the movement of sediment, in this case as it flows in 
pipes, did not use the shear velocity as an independent parameter. This pa- 
per, which does not seem to have attracted the attention it deserves, is note- 
worthy for its use of the drag coefficient of the sediment with respect to its 
fall velocity as an independent parameter. The purpose of this discussion is 
to show that Dr. Liu’s paper indicates that this parameter should have an 
equal application to the movement of sediment in open channels. 

In his equation (5) Dr. Liu presents the conventional relationship presumed 
to exist in a flow near a smooth boundary. At the time sediment just begins to 
move, such an assumption may be satisfactory. However, there is a growing 


mass of evidence that once sediment is moving these conventional relations no 
longer hold. 
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Discussion of 
“CHARACTERISTICS OF A LARGE THROATED SIPHON” 


by J. C. Stevens 
(Proc. Paper 1198) 


E. F. RICE, | M. ASCE.—The Sullivan siphon is of boldly novel design and 
unprecedented dimensions. As a solution to a difficult engineering problem, 
the design could hardly have been improved upon for economy, performance, 
or dependability. 

Some question might be raised as to the use of the terms “efficiency” and 
“discharge coefficient.” “Efficiency” can be defined as the ratio of estab- 
lished merit to unattainable ideal. The word should not be applied to siphon 
spillways, for few authors would agree as to what should be their established 
merit or unattainable ideal. If “efficiency” is to be regarded as the ratio of 
energy used in passing the designed flow to the energy available for such use, 
one could, with some justification, conclude that the “efficiency” is always 
100 per cent: All the energy available for producing flow is dissipated either 
in friction within the spillway or in eddies downstream from the structure. 

If efficiency is somehow to be regarded as a measure of the extent of the re- 
duction of unwanted friction losses, a similar conclusion could be reached. 
For each obstruction to flow that is removed or reduced, there will result a 
corresponding increase in the flow that may be passed through the siphon for 
a given head. This increased flow is a product of higher velocities which also 
result in increased friction losses until the total head loss again becomes the 
total operating head; the entire head is used up and no rational expression for 
efficiency results. If efficiency were to be defined as the ratio of energy re- 
covered to that expended, then the efficiency would be zero, for the siphon 
spillway cannot recover energy from falling water, but can only use it to ac- 
celerate the water along its way. Several writers have chosen to abandon any 
definition of efficiency in terms of energy in favor of more general criteria 
describing the utility of various designs. Frequently the “efficiency” is taken 
to be the coefficient of discharge, although in many cases designers have not 
agreed upon which dimensions should be used to determine this value. Asa 
result, coefficients of discharge have been reported from about 0.33 to over 
1.5, i.e., 33 per cent ‘efficiency” to 150 per cent for various designs of 
siphon spillways. The wide range of values and the extremely low apparent 
efficiencies of some of the best designs have been the sources of consider- 
able controversy. The difficulty arises in the selection of the area to be used 


in the equation Q= CAY2gH , where C is the disputed discharge co- 
efficient and A the area upon which the dispute hinges. If the cross-sectional 
area is constant there is no problem, but for the designer who hopes to in- 
crease the flow by making the lower leg divergent, as a draft tube, the choice 
becomes more difficult. Since the majority of designers have chosen the 
throat area as the basis for computing C it will be well to examine here the 


reason for the presence in occasional siphons of an “efficiency” exceeding 
100 per cent. 


1. Prof. of Civ. Eng., Univ. of Alaska, College, Alaska. 
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A siphon with a divergent lower leg fitted with a normally converging en- 
trance is somewhat like a bent Venturi tube. Equating energy at the forebay é 


to that of any point in the interior of the siphon using the outlet elevation as 
datum, 


2 2 


* * | 


Where Cx is a coefficient embodying all losses upstream from the point con- 
sidered, z is the elevation of the point considered, and V is the mean velocity 
of the liquid passing the cross-section in which the point is located. Writing 
the same equation for the cross-section at the exit 


z 2 
Rg 
| 
| 
y | 


where | 


H is the full operating head. C is the coefficient of discharge of the siphon 
and is frequently expressed as a percentage and called the efficiency. If 
there were no losses, C would be unity, corresponding to 100 per cent 
efficiency or 100 per cent conversion of the available energy into velocity | 


head. (If efficiency were to be defined as the ratio of output energy, % , to 


input energy, H, then efficiency, N, would be C2.) 

No coefficient greater than unity is possible if the area, A, used in the 
equations above is the outlet area of the siphon tube. However, most authors 
have used, instead, the area of greatest constriction (usually the throat area 
at the spillway crest) as the area in the equation 


Q=CA 


Thus, the arbitrarily reduced A demands a proportionally increased C, for un- 
changed flow and head. This new coefficient may be much greater than unity 
and “efficiencies” greater than 100 per cent may be reported. 

A coefficient based on any area other than the area of the issueing cross- 
section of liquid has no merit as a measure of the extent to which losses have 


Q= AV = CA 
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been eliminated. The erroneous use of the throat area is common, however, 
and is justified wherever the outlet area is not well defined, as in an incom- 
pletely submerged draft tube of trumpet shape which discharges horizontally. 

Mr. Stevens has long held that the efficiency of a siphon can be defined as 
the ratio between the average velocity across the throat of the siphon and the 
velocity theoretically obtainable under a head of 1 atmosphere. This measure 
of efficiency has the advantage that it accurately portrays the extent to which 
losses have been eliminated in the design of a siphon. Its legitimacy is based 
on two criteria: (1) the throat of a siphon is the highest, frequently most con- 
stricted, part, and (2) the flow through it cannot be greater than the flow due 
to a head of 1 atmosphere. To be completely correct the “1 atmosphere” 
should perhaps be reduced by the partial pressures of gases and the water 
vapor pressure. Should the term efficiency be used at all with regard to 
siphons, this definition will be found the least objectionable, even though it 
varies for varying heads on the structure. 

Since the true discharge coefficients (i.e., coefficients based on outlet 
area) of most siphons vary between 0.4 and 0.8, siphons compare with any 
other short tubes. It is something of a paradox that divergent tubes have low 
coefficients while straight or convergent tubes have high ones, although the 
reasons for this are fairly apparent upon reflection. Speaking of divergent 


tubes, King, Wisler and Woodburn, in a standard hydraulics text, have this to 
Say: 


“Experiments indicate that even under favorable conditions the value of 
Cy based on velocity at outlet is small. Venturi and Eytelwein, experi- 
menting with the flow of water in a tube eight inches long, one inch in 
diameter at the throat and 1.8 inches at the outer end, obtained a value of 
Cy of about 0.46. The lost head was therefore approximately 0.79 H. 
Even with this large loss of head the discharge was about 2 1/2 times the 
discharge from a sharp-edged orifice having the same diameter as the 
throat of the tube.” 


This last statement illustrates the fallacy of assuming that a low discharge 
coefficient is evidence of a low degree of hydraulic effectiveness. Unless 
based on outlet area, a coefficient of discharge is not even a valid basis for 
comparison of various designs. The tube cited above, though 2 1/2 times as 
effective as an orifice, would show less “efficiency,” even by Mr. Stevens’ 
criterion. 

From a practical standpoint, very little attention should be paid to reduc- 
ing the hydraulic losses by extreme refinements in design, because a slight 
increase in the cross-sectional area of the siphon will increase the flow 
enough to compensate for a profile of high loss. By this means any desired 
capacity may be achieved, up to a flow of nearly 200 cubic feet per second 
per foot of crest attained in the Sullivan siphon. 

The author’s correction of an obviously erroneous energy gradient is 
questionable. A cursory study of the literature would seem to Show that the 
velocity-head correction, a , rarely exceeds 2.0 even in natural streams, 
whose uneven distributions of velocity might imply the largest possible values. 
It remains to determine @ for tortuous, man-made conduits such as the 
Sullivan siphon. 

Prior to undertaking the model tests on the Sullivan siphon for Mr. 
Stevens, the writer made some velocity traverses in a siphon model of sim- 
ple design with the results shown in Figure 1. At section (2) the velocity 
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distribution closely obeyed the so-called free-vortex law, wherein the pro- 
duct of velocity and radius is constant for all stream-tubes. (While the dia- 
gram shown does not take into account the deviation of the streamlines from 
true radial lines, it nevertheless shows the relative velocities obtaining at 
each of the measuring points. The value of @ for this region (assuming free- 
vortex distribution) is calculated to be 1.5. The values of @ for sections (3) 
and (4) are hardly greater, and at sections (1) and (5), approaches unity. 
Even if the velocity distribution in the Sullivan siphon was equivalent to a uni- 
from variation from zero on one side to maximum at the other, would be no 
greater than 2. 

To postulate that a is much higher than 2 at any section would be to imply 
eddies so strong as to reverse the flow in some portions of the siphon. No 
such reversals were present either in the low head or high head models, and 
it is hardly possible that they existed in the prototype—certainly not at 
piezometer set No. 1. 

Therefore, uneven velocity distribution is not likely to be the true explana- 
tion for the misaligned energy gradient as measured with the piezometer 
tubes. 

The 11 per cent error of prediction of design flow of the Sullivan siphon by 
model studies warrants some explanation, especially since the prediction errs 
on the unsafe side (which is unusual in hydraulic models). Such an error is 
often acceptable since it may be well within other factors of uncertainty. In 
this case, the 11 per cent error includes not only the “model error,” but also 
considerable variation due to design changes, for Mr. Stevens wisely used the 
model as a basis for modifications rather than as a mere “checking” device. 

Figure 2 shows the changes made in the design. Both the addition of the 
air relief valves and the redesign of the outlet were dictated by factors deter- 
mined by the model’s priming characteristics. Under conditions of high tail- 
water, the sheet of falling water could not exhaust the air, for the low veloci- 
ties prevailing during free fall over the crest were unable to sweep the 
bubbles out through the descending leg. Also, the rising forebay displaced 
air within the siphon, causing a rise in the air pressure inside, inhibiting 
priming by depressing that part of the forebay within the siphon inlet. Air 
check valves proved able to correct the priming; the change in the descending 
leg was made to assure higher velocities for carrying air through the siphon. 
While they affected the discharge coefficient adversely, these modifications 


were expected to improve priming and to make possible much better low-flow 
characteristics. 


| 
| 
| 
is 
| 
| 
| 
| 


October, 1957 


nd SV d31S31 Sv 


ML \- 


O15 GE 

6 


AVMT1dS NOHdIG JLLAaWV TIM 


1417-18 HY 5 
‘N | 
& \e 3\ | 
AN 
“a - \\\ \\ \\ \\ 
\\\\\\ \ \ | 
| 
Q\ | 
y 
f | 
EM 
| 
| 
VA 
x 
| 
| 
| 
|: bw | 


Discussion of 
“OCEAN WAVE FORCES ON CIRCULAR CYLINDRICAL PILES” 


by R. L. Wiegel, K. E. Beebe and James Moon 
(Proc. Paper 1199) 


T..E. STELSON, | J.M. ASCE.--Large scale model tests such as the ones 
described in this paper are of great value in obtaining immediate design data 
as well as providing a basis for evaluating both theoretical relationships and 
the results of smaller scale model studies. When so much time, money and 
skill result in data that show such tremendous variation, non-uniformity and 
scatter, however, something is lacking, either in experimental method or 
theoretical understanding. Measured values for apparently identical experi- 
mental conditions vary by factors of more than ten as shown in Figures 11, 
12, and 13. A paper such as this provides the stimulus to evaluate our 
present knowledge in this field considering both basic theory and supporting 
experimental data. 

One example of such an area of deficient understanding is that of inertia 
drag. The added mass coefficient, Cyy, as defined in equation (3) would have 
a theoretical value of one* for a circular cylinder of infinite length, accelerat- 
ed in an ideal fluid of infinite extent in irrotational motion at rest at infinity. 
This theoretical coefficient was erroneously given as two on pages 22 and 23. 

The experimental data summarized in Figure 13 would then show that 91 
per cent of the measured values for the mass coefficients were in excess of 
the theoretical value. The mean experimental value of 2.5 is 150 per cent in 
excess of the theoretical value. Furthermore, the data presented in Table 3 
shows the maximum coefficient to be more than seven times the minimum 
coefficient. 

Some deviations from the theoretical value would be expected because of 
divergence in physical conditions such as (1) limited length of the test cylin- 
der, (2) wake formation behind the pile, (3) non-uniform movement within the 
wave and incorrect determination of acceleration, and (4) presence of the free 
surface as well as the variation in turbulence as mentioned by the authors. 
None of these items, however, offer any reasonable explanation for such great 
diversity and scatter of data unless it could possibly be item (3) above. 

The coefficient of drag as shown in Figure 11 also showed tremendous 
variation with the larger values being more than twenty times the smaller 
values for the same pile and Reynolds’ number. 

The forces computed using mean values of drag and mass coefficients and 
plotted in Figures 16 and 17 still show considerable scatter indicating that 
(1) measured forces are incorrect, (2) mean coefficients are incorrect, or 
(3) understanding of the phenomena is inadequate. The tendency for computed 
forces to be higher than measured forces when forces are low and for 


1. Asst. Prof. of Civ. Eng. and Acting Dept. Head, Carnegie Inst. of Tech- 
nology, Pittsburgh, Pa. 

* Milne-Thomson, L. M. Theoretical Hydrodynamics, 3rd ed., The Mac- 

millan Company, New York (1955) P. 233. 
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computed forces to be lower than measured forces when forces are high (a 
dangerous tendency), may point to the cause of some of the irregularity. 
Force due to fluid inertia varies as the height of the waves whereas force due 
to drag varies as the square of the height of the waves. Thus, acccleration 
dependent forces seem to be too large and velocity dependent forces too small 
since the larger forces are associated with the higher waves. This again 
would indicate that mean coefficients of mass are too large and that mean 
coefficients of drag are too small. Lower values for the mass coefficients 
would make them more consistent with theory, and the possibility of higher 
drag forces has been stated by the authors on page 22. 

A possible conclusion from this paper is that if the experimental measure- 
ments are correct, a re-evaluation of all concepts of wave movement, fluid 
inertia and viscous drag should be made. They certainly do not begin to pro- 
vide an adequate basis for reasonable correlation for the reported data. So 
much divergence cannot be tolerated in most engineering situations unless 
large safety factors are consistently used. Thus designs tend to be conserva- 
tive and expensive and fewer projects are economically feasible. 
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Discussion of 
“THE ESTIMATION OF THE FREQUENCY OF RARE FLOODS” 


by Benjamin A. Whisler, M. ASCE, and Charles J. Smith, A. M. ASCE 
(Proc. Paper 1200) 


LEO R. BEARD.!—The equation given by the authors for relating the 
probability of monthly flows to annual flows is apparently based on the as- 
sumption that flow probabilities during each month of the year are equal for 
equal flows, which condition is rarely, if ever, attained under actual condi- 
tions. As a consequence of the inequity of the monthly probabilities in most 
regions of this country, the relationship is complex, and it is believed that 

the authors’ conclusions should be reexamined in the light of empirical rela- 
tionships among the frequency curves for the various calendar months. 

These relationships will be different for different geographical regions. 

The authors’ findings that frequency curves of annual flows curve upward 
on logarithmic-probability paper and that those of monthly flows are straight 
are not borne out by the examples given and might well be tested using data 
from other parts of the United States. The Washington District, Corps of En- 
gineers, has made such a test for the annual-event curve and has found2 that 
frequency curves of annual maximum flows for 108 of the streams throughout 
the United States having the longest records of unregulated runoff average 
very nearly a straight line. While frequency curves of annual flows on 
streams in some regions, particularly in the northeastern United States, do 
have a tendency to curve upward, those on streams in some other regions 
show the opposite tendency. Whether these regional irregularities are due to 
regional characteristics or to accidental excess or deficiency of floods or 
droughts in the period of record is a question. 

One of the advantages claimed for the authors’ method is that it is based 
on more numerous data than the annual-event method, but this advantage may 
not be as great as appears at first sight. The writer questions whether a 
curve based on monthly maximum flows can actually be extrapolated much 
more accurately than a curve based on annual maximum flows. Extrapolation 
must start at the same magnitude, since the largest observed flow is common 
to both series. Also, by the authors’ table it can be seen that while the 
monthly series has 12 events (1 percent) in 100 years above a given magni- 
tude, the annual series has 11.4 on the average above that same magnitude. 
Therefore, in the region where the trend of the frequency curve is critical 
for extrapolation purposes, the numerical advantage of the monthly-event 
curve is small. 

It is hoped that the advantages of statistical techniques over graphical 
techniques are greater than the authors infer. Certainly the authors have not 
exaggerated the unreliability of estimates based entirely on a single record. 
No method, no matter how ingenious, can make up for an accidental deficiency 


1. Hydr. Engr., Corps of Engrs., U.S. Dept. of the Army, Sacramento, Calif. 
2. Stream Flow Volume-Duration-Frequency Studies, June 1955, Office of the 
District Engineer, Washington District, Washington 25, D. C. 
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or excess of large floods during the period of record at one location, without 
comparing the experience at that location with experience at other locations. 
Statistical methods provide an effective vehicle for comparing frequency 
characteristics and for differentiating between accidental and characteristic 
variations. The writer believes that statistical techniques, improved from 
time to time by new ideas, such as those suggested by the authors, and tested 
in the light of experience, as is being done under the Corps of Engineers 
Civil Works Investigations Program, will pave the way to much more precise 
frequency estimates in the future. 


MANUEL A. BENSON,! A.M. ASCE.—The authors’ ideas are ingenious and 
certainly well worth investigation. Their discussion of some present-day 
methods of using mathematical statistics in flood-frequency analysis is dis- 
cerning and points out the weaknesses and assumptions involved. The writer 
would like to express complete agreement with and reemphasize their obser- 
vations that 


*, .. many investigators utilized mathematical techniques involving the 


coefficient of variation and the coefficient of skew to fit skew curves to 
peak annual flow data that plotted as a curve. This technique becomes 
little more than a mathematical method of fitting an empirical curve to 
the data and is little, if any, better than drawing such a curve graphical- 
ly by eye through the plotted points, and may actually be worse.” 


The authors’ contention that the series of annual flood peaks is not a 
proper series for statistical analysis is debatable. They have said, “It is 
faulty statistical technique to use selected data in studying any phenomenon 
which is presumed to follow the basic laws of probability.” There is nothing 
inherently incorrect in selecting data to be subjected to statistical analysis, 
particularly such selection as is represented by the highest event in a time 
series. Some types of selection are proper, others are not. There may be 
something inherently incorrect, however, about assuming that any phenome- 
non follows any of the basic laws of probability if basic law is intended to 
mean a frequency distribution. This applies equally to phenomena for which 
“all” the data are considered. 

It must be recognized that a “flood peak” per se is not as clear cut an item 
as most of the phenomena ordinarily subjected to statistical analysis. A flood 
peak needs to be defined before it can be determined. By “flood peak” most 
people would mean the single momentary maximum or daily maximum occur- 
ring during a rise caused by a flood of damaging or notable magnitude. The 
annual flood comes as close to fitting this definition as is possible. The con- 
cept of “floods above a base” attempts to include all the high or damaging 
floods regardless of the period of time within which they occur, but in doing 
so introduces uncertainty because flood crests often occur in rapid succession 
and may be part of the same general rise. The authors’ definition of “all” 
peaks will mean the use of many events occurring at low flows, most of which 
are not ordinarily considered as flood peaks at all. It also involves consider- 


ing each of the separate crests during a flood rise as being independent events. 


There does not seem to be any reason for expecting that a series of such 
events, defined arbitrarily, should follow a log-normal or any other statistical 
distribution, despite the fact that “all” such events are being considered. As 


1. Hydr. Engr., U. S. Geological Survey, Wash. 25, D. C. 
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the authors realize, any such contention must be shown to conform to actual 
streamflow data in order to justify the use of a fitted distribution which aver- 
ages the data within the defined range. Whether the distribution is or is not 

a good fit within the range of data, extrapolation is risky and not justified un- 
less as a minimum condition the evidence is numerous and shows a definite 
and undeniable tendency. Even then, the possibility exists of a change in 
relationship beyond the range of data. 

Langbein (1949) has demonstrated the relation between the annual flood 
series and “floods above a base.” The same relation holds true between the 
annual series and the series of “all” peaks, since “all” peaks represent the 
floods above a zero base. As Langbein shows, there is virtually no difference 
between the frequency curves derived from the two above a recurrence inter- 
val of about 10 years. The additional floods included in the partial-duration 
(or complete-duration) series serve only to add information at the lower end, 
where it is not badly needed. For this reason there is no advantage in con- 
sidering “all” floods or monthly floods, unless the relationship shown by such 
a series is to be extrapolated beyond the range of data, as the authors have 
done. 

The plotted data for the three stations analyzed by the authors do not show 
clear-cut evidence of having produced better methods of fitting the data than 
are now available. Their final computed curves of annual peak flow do not 
appear to give any better fit of the data than do the computed straight lines 
of best fit in the case of Lehigh River and Iowa City (authors’ fig. 1, 3). 

Their curve is not as good a fit for Tunkhannock Creek (authors’ fig. 2). 

The curves for Lehigh River (authors’ fig. 1) have been cited as demon- 
strating “the great errors in the determination of the recurrence interval of 
rare floods which are made in assuming that annual peak flows follow the 
logarithmically normal laws of probability.” It would be equally as valid to 
show by means of the curves for Tunkhannock Creek (authors’ fig. 2) “the 
great errors which are made in assuming that monthly peak flows follow the 
log-normal laws.” Actually, neither the monthly nor annual peaks are demon- 
strated by the fit of the data as being representative of the log-normal law. 

Data for all the peaks on Lehigh River show a slight positive curvature 
(concave upward) for Lehigh, a definite negative curvature for Tunkhannock. 
The monthly peaks for Iowa River exhibit a slight negative curvature. It is 
therefore not at all certain, particularly in view of only three gaging-station 
records, that ‘all” peaks or monthly peaks exhibit a log-normal distribution. 

The method has been tested by the writer using unregulated records at 8 
stations in California (unpublished records of monthly peak discharges avail- 
able at District office of U. S. Geological Survey, Menlo Park, Calif.) and 6 
stations in Great Britain (1937-45, 1945-53, 1953-54). For these stations 
records of momentary monthly and annual peaks are available. Daily peaks 
have but a limited use in engineering design involving storage, whereas 
momentary peaks are the important factor in the design of most hydraulic 
structures and in problems of flood prediction and flood damage. For this 
reason no attempt was made to analyze “all” peaks and investigation of the 
authors’ methods was confined to momentary peaks. 

The longest record is that for Kern River near Bakersfield, Calif., for 
which monthly peaks have been published in water-supply papers from October 
1912 to December 1950. The other California stations have records from 10 
to 18 years in length. The records in Great Britain are 17 years long. 

The curves for Kern River near Bakersfield, Calif., are shown in figure 1 
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and those for River Avon at Evesham, England in figure 2. The data for 
River Avon are typical of those for all the English stations and the remainder 
of the California stations. In all cases the actual monthly peaks appear to de- 
fine an S-curve. In all cases (except Kern River) the top points are well to 
the right of the least squares line through the monthly peaks. At every station 
but Kern River, the actual annual peaks are well to the right of the trans- 
formed annual peak curve. (Note that in figures 1 and 2 only one in ten of the 
monthly peaks are shown within the middle range, for clarity). 

The following table shows the recurrence interval from the theoretical 
(transformed) annual peak curve, of 10-, 50-, and 100-year peaks on the an- 
nual peak least-squares line. 

The table on page 25 shows the extremely large difference between the an- 
nual peak least-squares lines and the authors’ derived annual peak curves. 
For example, at River Avon, the discharge with 100-year recurrence interval 
on the least-squares line has a 14.3-year recurrence interval on the authors’ 
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annual peak curve. The closer agreement, shown by the table, between the 
two curves for Kern River should not hide the fact that the actual annual 
peaks define a quite different relation than either of the two curves (see 
fig. 1). The writer considers the least-squares line through the peaks to be 
a better means of predicting future annual peaks than the authors’ trans- 
formed curve. However, a graphically-fitted line through the upper points 
would be even better. 

Some explanation seems necessary to account for the large discrepancies 
between the actual annual peak floods experienced and the authors’ derived 
annual flood curves. There appears to be a basic accumption in making the 
transformation from monthly to annual peaks that is not borne out by the 
facts. The statement reads, “If P,, is the probability of a given event occur- 
ring in any one month, the probability, Py, that it will occur in a year is 
much greater and can be expressed as Py = 1-(1 - Pm)!2.” In actuality, the 
probability of a peak of a given magnitude being a monthly peak is not the 
same in any two consecutive months and is far from the same in 12 consecu- 
tive months. This is well illustrated by the record on the Kern River near 
Bakersfield, Calif., where all 39 annual peaks occurred during the 8-month 
period from November through June. Monthly peaks of the magnitudes oc- 
curring in the 4-month period from July through October have only an ex- 
tremely small chance (if any) of being the annual peaks. Or, working back- 
wards, annual peaks which occur with high frequency (high probability) in 
November-June have practically no probability of occurring during July- 
October. It might be thought that this situation could be handled by computing 
the probability of, say, a 1,000 cfs peak occurring in each of the 12 months 
and multiplying these. However, even this is not correct, since the principle 
of independence of probabilities which is assumed in the combining process 
does not actually exist. For example, if all annual peaks occurred in May, 
the probability of May peaks and annual peaks would be identical, and all the 
other monthly peaks would have no effect on the probability of the annual 


peak. The situation at Bakersfield and any other river is analogous, though 
not as extreme. 


In summary: 


1. The authors’ method of extrapolation beyond the range of data is de- 
pendent on the validity of the log-normal distribution for “all” peaks 
(approximated by distribution of monthly peaks). Conformance to this 
distribution cannot be predicated by theory and is not borne out by 
actual streamflow records. 


- The method of transforming from monthly to annual probabilities is 
based on assumptions which are not in accord with the hydrologic facts. 
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Discussion of 
“HIGH HEAD CAVITATION TEST STAND FOR HYDRAULIC TURBINES” 


by W. G. Whippen and G. D. Johnson 
(Proc. Paper 1201) 


LESLIE J. HOOPER, ! M. ASCE.—The authors are to be complimented on 
their clear presentation of a carefully worked out design which seems to be 
admirably adapted to the purpose intended. The provision in the design for 
mounting the water wheel casings and dynamometer in one frame work seems 
to be particular noteworthy. In this way it is hard to see how deflections 
caused by changes in pressure or changes in the loading on the pipe can pos- 
sibly affect the alignment of the shaft and friction readings of the rotating 
parts. This is an important feature bearing upon the repeatability of tests. 

With respect to repeatability of tests it is noted that the repeatable 
accuracy of any test point appears to be within + 1/4 of 1%. It is presumed 
that this refers to a single test set-up where the gate and speed are read- 
justed to the same operating conditions. It would be interesting to know what 
the discrepancies might be when a unit has been tested, completely torn down 
and then reinstalled at a later date and tested again. 

There are one or two questions that come to mind that probably can be 
very easily answered. It is noted at one point that “since equipment is avail- 
able to cool both fresh and de-aerated water, the effect of air content upon 
performance and cavitation may be studied.” Inasmuch as no provision has 
been made in the design to re-absorb air that might be liberated in a passage 
through the test circuit it is difficult to see how the studies on the effect of 
air content can be made when the air content is not controlled. 

A second question that arises is concerned with the prediction of full scale 
prototype efficiency of model test results. It is noted in Fig. 4 that model 
test efficiencies of the order of 91.5% are being obtained. Assuming a 12" 
diameter model, a 240"' diameter prototype and the Moody step-up formula 
with the 1/4 power of the diameter ratio there results a predicted efficiency 
of 96%. If the 1/5 power of the diameter ratio is assumed the resulting 
efficiency is still 95.4%. Either of these results is very complimentary to the 
design of the water wheel but do seem a little difficult to obtain by field test- 
ing methods. 

There is no question that laboratory tests can be made more accurately 
than field tests because the control and possibilities of calibration are much 
better. However, there is still a place for field tests inasmuch as a model 
test can never duplicate all of the conditions that exist in the field and it 


would seem that the step-up formulas now in use still leave a little to be 
desired. 


B. L. VANDER BOEGH.2—The S. Morgan Smith high head test stand offers 


. Prof. of Hydraulic Eng., Worcester Polytechnic Inst., Worcester, Mass. 


2. Hydr. Engr., Newport News Shipbuilding and Dry Dock Co., Newport News, 
Va. 
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a means heretofore unavailable by which the answer to several baffling prob- 
lems may well be determined. 

Mention ‘= made of the lower head flume. A comparison of performance 
test results of homologous models between the low head and the new high 
stand would certainly be of interest. A later comparison with a prototype 
test might indicate the need or, conversely, the lack of a requirement that all 
model testing be done at field heads if possible. 

A visual comparison of cavitation formation at high and low heads is not 
possible inasmuch as this feature is not incorporated in the new stand. It 
would seem that this comparison would be of definite value in outlining the 
surfaces to be overlaid with stainless steel, and in the observance of tip 
cavitation as well as trailing edge flow patterns. 

It is assumed that provision for the study of the hydraulic moment and 
forces on the vanes of a Kaplan model has been made in the low head stand 
inasmuch as no mention is made of this feature in the high head stand. This 
might well have been considered in the high stand as the moments are known 
to vary with sigma, and this stand should give more accurate sigma values. 

It would seem that a filling of the system with city water for each test 
would make it difficult to maintain a uniform air content over the extent of 
the test. At low sigmas, it has been suggested that the air content may 
definitely affect the performance to varying degrees. A storage tank for the 
reuse of the water might seem desirable. 

The cooling system for the stand is definitely interesting. It would be 
helpful to know the effect of temperature change of the water on the cavitation 
results, (inasmuch as the water vapor pressure is included in the sigma 
formula) with a suggestion or recommendation as to the reliable temperature 
range for accurate cavitation testing. 

Many of the newer test stands use head and discharge tanks. The question 
is raised as to whether the flow conditions into the case inlet of the model 
are affected without the use of an inlet section or penstock. 

In general, the new stand seems to be well designed, incorporating pro- 
visions for exceptionally accurate testing, and the test results should be both 
interesting and reliable. 
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Discussion of 
“IS THE WRITING OF FLOOD INSURANCE FEASIBLE?” 


by John F. Neville 
(Proc. Paper 1202) 


E. M. LAURSEN, ! A.M. ASCE, and A. TOCH,2 J.M. ASCE.—The three 
papers (Proc. Papers 1164, 1165, and 1202) of this series describe very 
clearly the difficulties involved in the various aspects of writing flood in- 
surance—the first two papers dealing primarily with the technical problems, 
this paper with the practical problem of selling insurance. The conclusion 
that appears to be inescapable is that such insurance would be possible tech- 
nically but impossible practically. It would seem, however, that the conclu- 
sion “insurance against the peril of flood applicable to fixed property cannot 
successfully be written” is based on a possibly shaky premise as to the type 
of insurance program to be evolved. That flood insurance must be written on 
the basis of the average annual loss incurred by each property is not pre- 
ordained. However, if this assumption is made it does become necessary for 
the insurance company to make hydrologic and damage studies (as outlined by 
Foster in the second paper of the series) which will cover each individual 
property to be insured. Not only is the average annual loss high for any 
property on the flood plain, but there is then a sizable “expense loading” 
added to the premium. 

Average annual losses should be thought of, quite properly, as part of the 
price paid for the benefits derived from occupation of the flood plain. The 
cost of repairing damage caused by flood stages which occur frequently can 
and should be paid just like taxes or any other cost. 

The following hypothetical problem of determining the location of steam- 
power plant will serve to illustrate this assertion. Given two alternate loca- 
tions, one on the flood plain next to the river and another on high ground out 
of danger of flooding, certain advantages will be inherent in each site. Con- 
struction and operating costs for the condenser-water supply will be less at 
the flood-plain location, as may be storage costs for the coal supply. Costs 
attributable to flood damage will be least at the high-ground location. Other 
costs, such as land acquisition, foundations, and spur trackage, will be less 
for one or the other site depending on local conditions. An economic analysis 
will determine the best site, with average annual flood damage taken into ac- 
count either as an operating cost or by inclusion of proper structural safe- 
guards. 

Thus there is no advantage in having insurance against floods of frequent 
occurrence, since the premiums must necessarily be greater than the 
“losses” due to those floods. Insurance against the unusual damage caused 
by a very large—and therefore very infrequent—flood may be desirable, how- 
ever, because a considerable reserve fund would otherwise be needed in 


1. Research Engr., Iowa Inst. of Hydr. Research, Iowa City, Iowa. 


2. Research Engr. and Research Associate, Iowa Inst. of Hydr. Research, 
Iowa City, Iowa. 
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order to overcome such a catastrophic loss. Insurance can be substituted for 
all or part of this reserve. The amount of insurance does not need to be re- 
lated precisely to the damage which actually occurs, but can be an arbitrary 
amount which will cushion the effect of the flood damage. Since only floods of 
rare frequency would be involved, the premium should be much smaller than 
would be required to cover the average annual loss. Stage-probability and 
stage-damage studies for the individual property do not need to be as precise, 
and the burden of these studies can be placed on the insured. The policies 
could be written simply on the occurrence of floods of various magnitudes as 
determined at the nearest gaging station, and stage-probability relationships 
for these stations would form the basis for the determination of risk and cor- 
responding premium per dollar of coverage. 

The argument for considering flood losses as a cost of doing business is 
readily apparent in the case of a steam-power plant. It is equally true for 
any other business. If the benefits derived from a flood-prone location do 
not outweigh the costs resulting from flood damages, a move to high ground 
is indicated. If the benefits do outweigh the average annual “losses,” it would 
be utter folly to pay an insurance premium greater than these ‘losses” except 
as a cushion against extreme cases of catastrophic losses due to very rare 
floods. Farming would differ from other business enterprises only in that the 
timing and duration of a flood might assume, to the insured, as much impor- 
tance as stage. It is difficult to see any benefits of a flood-prone location 
accruing to dwellings except that of cheap land. If the damages exceed this 
benefit, moving the house to high ground or putting it on stilts would seem 
preferable to buying flood insurance for frequent floods. 

In order to protect the insurance company from opportunists, two safe- 
guards could easily be used: a time lapse of several months between the date 
of payment of premium and the date the policy goes into effect, to eliminate 
policies bought after a definite flood prediction has been made; and a “house 
cut” of sufficient magnitude to discourage mere gambling on the probable 
river stage. Long-term policies could be encouraged, on the other hand, by 
reducing the premium every year for several years on policies kept in force 
continuously. 

It would be interesting to know whether the type of insurance proposed 
here would be more feasible and practicable from the underwriters’ view- 
point than the average-annual-loss type assumed by the author. 
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Discussion of 
“THE EFFICACY OF FLOOR SILLS UNDER DROWNED HYDRAULIC JUMPS” 


by Ahmed Shukry 
(Proc. Paper 1260) 


A. RYLANDS THOMAS,* M. ASCE.—The provision of a sill, lip wall or de- 
flector, at the end of the floor to reduce scour downstream of low-head struc- 
tures is regarded as standard practice in many countries. The height of sill 
is generally about one-tenth of maximum tailwater depth. It is believed that 
this device has been used on many of the barrages and other hydraulic works 
in Egypt, so the author’s preference for a sill at an intermediate position is 
all the more surprising. In this he appears to have ignored the evidence of 
Fig. (2), which shows the beneficial effect of the end sill in reducing scour 
and preferred that of the velocity distribution of Figs. (8) and (11), and cer- 
tain earlier tests not quoted. 

It would, however, appear that the rigid floor in Figs. (8) and (11) was con- 
tinuous beyond section V and therefore the flow pattern at that section was not 
representative of conditions where scour occurs immediately downstream of 
the sill. There is no reason why the effect of an end sill in the case where the 
hydraulic jump is drowned should differ from that where the jump is clear, 
provided that in both cases measures are taken to normalize the flow up- 
stream of the floor end. These measures may include the provision of a sill, 
baffle blocks or baffle piers at an intermediate position and their function is 
separate and distinct from that of the end sill which is required in addition. 

Rehbock® showed that the effect of an end sill was to create a roller with 
horizontal axis which swept bed material back towards the end of the floor, 
thus reducing the depth of scour at that point, and this action has been con- 
firmed in numerous experiments with models, including those of Butcher and 
Atkinson.§ Experience in India and Pakistan with full size structures has 
been favourable, though in the case of larger works the effect is generally 
masked by the presence of concrete blocks laid downstream as an additional 
safety measure. It is only when these blocks have sunk in the sand or have 
been washed downstream, which should not occur under design conditions, 
that the effect of the end sill in reducing scour becomes apparent. 

The author’s contention that an end sill is not beneficial would carry more 
weight if supported by observations of scour in models and, if available, on 
full size structures. It would also be useful if details were given of the gate 


regulation, discharge and tailwater depth in the tests to which Figs. (8) to (11) 
relate. 


- Consultant, Barclay’s Bank, Abbey House, London, England. 

. “Bekampfung der Sohlenauskolkung bei Wehren durch Zahnschwellen” by 
Th. Rehbock. Zeitschrift des Vereins Deutscher Ingenieure, No. 69, 
p. 1382, 1925. 

- Reference quoted on page 9 of Paper. 
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Discussion of 
“A STUDY OF BUCKET-TYPE ENERGY DISSIPATOR CHARACTERISTICS” 


by M. B. McPherson and M. H. Karr 
(Proc. Paper 1266) 


E. A. ELEVATORSKI,* J.M. ASCE.—The authors have presented useful 
information on the characteristics of solid roller bucket-type energy dissi- 
pators. There is a regretable lack of information on this subject and cer- 
tainly more is needed. In part, the lack of information can be attributed to 
the difficulty of relating the many variables which affect the performance of 
the bucket. It is the purpose of this discussion to list the pertinent factors 
and relate them in a convenient dimensionless form. 

Probably the two most important factors in the design of a solid roller 
bucket are the size of the bucket radius and the tailwater depth. A combina- 
tion of the minimum allowable bucket radius and minimum tailwater depth 
exists, which will permit satisfactory bucket action. Below these limits, the 
bucket will not perform satisfactorily. Buckets larger than the minimum 
allowable bucket radius will only slightly improve the bucket performance 
while being far more costly to construct. It is therefore evident that there 
is a minimum bucket size which will still perform properly and be the most 
economical design. 

This writer proposes that the minimum allowable bucket radius is a func- 
tion of the following variables. 


(1) 

In equation (1), vy and dj, represent the computed velocity and depth, 
respectively, for a point at the intersection of the spillway face and tailwater 
surface. The symbols s, g, and y represent the slope of the spillway face, 
acceleration due to gravity, and vertical distance between the bucket lip and 
channel bed, respectively. 


With the aid of dimensional analysis, the variables in equation (1) can be 
reduced to dimensionless form 


) (2) 


s 
— 
y 


gd, 


Studies conducted by the Bureau of Reclamation(1) indicate that both s and 
y do not, within ordinary limits, affect the bucket size. It was pointed out, 
however, that frictional resistance on the spillway face affects the minimum 
allowable bucket radius, as a result of the velocity reduction. Eliminating, 


= from equation (2) we find 


a. Hydr. Engr., Bureau of Reclamation, U. S. Dept. of the Interior, Salt Lake 
City, Utah. 


| ASCE 1417-33 

| 

| 


1417-34 HY 5 October, 1957 


“mine (3) 


where 


gd, 


This is contrary to recommendation (d) of the authors, which states that 
the size of the bucket exerts no influence on bucket action. Investigations by 
Okada and Ishibashi(2) substantiate the functional relationship given by equa- 
tion (3) and further state that “the discharge per unit width of spillway is the 
predominant influence on the bucket size.” 

When the tailwater is insufficient, the high-velocity flow entering the 
bucket will break through the supernatant water blanket, leaving the bucket in 
the form of a jet. The depth at which this action takes place is known as the 
“sweepout” depth. Considerable scour occurs at the point where the jet 
strikes the riverbed, making this type of action undesirable. The tailwater 
depth just safely above the sweepout depth will, hereafter, be called the mini- 
mum tailwater depth. Factors affecting the minimum tailwater depth are 
given in equation (4). 


teWe f\ ,41,5,%,y,¢) (4) 


Again, by dimensional analysis, equation (4) is reduced to dimensionless 
form 


Within ordinary limits, : can be neglected. Consequently, 


teWemin,. = ) (6) 
dy dy 


Equations (3) and (6) express the minimum allowable bucket radius and 
tailwater depth as a function of dimensionless parameters. A relationship 
between the parameters can be determined with the use of laboratory models. 

Some designers may be concerned with the degree of submergence or 
maximum allowable tailwater depth for satisfactory bucket performance. 
Evidence obtained from both laboratory and prototype tests seems to indicate 
that the condition of “diving” flow, for solid roller buckets is impossible ex- 
cept perhaps for rare cases. 

Of interest to the designer will be the dimensions of existing solid roller 
buckets listed in Table 1. As noted in the table, almost all of the buckets 
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have 45-degree lip slopes. Experiments indicate that a 45-degree lip slope 
renders best results.(3) If the lip slope is made steeper than 45-degrees, the 
impact of the flow on the lip wall will disrupt roller action, while for flatter 
slopes the high-velocity flow downstream from the bucket was found to be 
submerged. 
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Discussion of 
“MEASUREMENT OF SEDIMENTATION IN TVA RESERVOIRS” 


by E. H. McCain 
(Proc. Paper 1277) 


LLOYD C. FOWLER,# J.M. ASCE and ROBERT H. LIVESEY.>—There has 
been, in recent years, considerable improvement in the equipment and pro- 
cedures available for hydrographic surveys; however, there has been little 
correlation or standardization in the use of the newer more accurate and 
efficient equipment and procedures. The author’s dissemination of informa- 
tion on the procedures of the TVA surveys should contribute to the general 
advancement. 

The adoption of the sonic sounding to replace the laborious and time- 
consuming sounding by hand line is unquestionably a major forward step in 
sub-aqueous surveys. Instead of a series of point elevations, the sonic 
sounder produces a continuous and permanent record of the bottom configura- 
tion. High frequency sounders will accurately show the surface of low- 
density deposits, and low frequency sounders can penetrate deposited sedi- 
ments, recording the surfaces between superimposed strata of increasing 
density. 

Sounder accuracy of 13 inches between depths of 3 to 50 feet, and +6 inches 
for greater depths, is adequate for most hydrographic surveys over relatively 
flat or gently sloping terrain. However, two sources of error are introduced 
immediately when traversing steep slopes or vertical banks; both originate 
with the transmission of the sound wave in a diverging cone from the trans- 
ducer (antenna) and increase with depth of water. First, an error in depth is 
introduced because the recorded echo was reflected from the shallowest sur- 
face within the radius of the cone, and second, a positioning error due to the 
increasing radius of the cone with depth. Assuming, for example, the cone 
apex is 30° (it can vary from 7 1/2° to 50° depending on the type of trans- 
ducer), then at a depth of 50 feet the recorded echo could have theoretically 
originated from the shallowest depth anywhere within a circle 27 feet in 
diameter, discounting any signal loss due to absorption or radiation. 

The several districts within the Missouri River Division of the Corps of 
Engineers have had considerable experience in the development and improve- 
ment of the equipment and procedures used in hydrographic survey work on 
both rivers and reservoirs. Their experience has indicated that the proper 
craft for work of this nature depends on the conditions imposed by the river 
or reservoir. A skiff with an outboard motor is best used in shallow water 
along the shores of a reservoir, in its headwaters, or in the open river while 
a relatively large inboard cabin cruiser capable of withstanding heavy weather 
and seas may be best for the sounding of long and deep reservoir ranges. 


a. Hydr. Engr., Missouri River Div., Corps of Engrs., U. S. Dept. of the 
Army, Omaha, Nebr. 

b. Head, Reservoir Sedimentation Unit, Omaha Dist., Corps of Engrs., U. S. 
Dept. of the Army, Omaha, Nebr. 
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Some craft may be specially designed to meet specific needs. For instance, 
Figure 7 shows a shallow draft river survey boat designed especially to ob- 
tain data on sediment transport, velocity distribution, and bed form roughness 
in the open river. 

In the boats, launches, or skiffs, used regularly for sounding, the trans- 
ducers should be mounted within the hull. One exception is a scow bow boat 
(see Figure 8). The skiff illustrated in Figure 1, is not generally considered 
as having a scow bow. Transducers mounted inside the hull of scow bow 
boats have proven to be unsatisfactory for sonic sounding because of the pas- 
sage. of entrained air under the hull. The transmission of sound through this 
aerated lamina results in a severe reduction and usually a complete attenua- 
tion of the transmitted signal and returning echo. With a scow bow boat the 
outrigger “shoe” enclosing the antenna is essential (Figure 8). 

Seaworthiness, maneuverability, and draft are relatively constant factors 

inherent in the design of a specific survey boat; however, weight and bulk of 
equipment, plus survey “know how” are variables which directly influence the 
safety, accuracy, and economy of reservoir and river hydrographic surveys. 
Considering the skiff described by the author, removal of the outrigger “shoe” 
with its enclosing well would produce a sizeable reduction in both weight and 
bulk. Early experience with mounting the “shoe” in an open well, similar to 
that illustrated in Figure 2, resulted in a reduction in maximum boat speed of 
almost 50 percent and severely limited the maximum sounding speed by intro- 
ducing aeration below the transducers. All of the boat speed was returned 
and higher sounding speeds were possible when the bottom of the well was 
closed by a smooth steel plate mounted flush with the bottom of the skiff. 
The “shoe” was then placed just above the steel plate and submerged in a few 
inches of water. However, constant attention was required to keep the outside 
surface of the plate smooth to prevent the accumulation of entrained air which 
would obliterate the sounding trace. 

Even greater efficiency and excellent results have been obtained by replac- 
ing the awkward well and “shoe” by separate transmitting and receiving mag- 
neto-restriction transducer units mounted in individual tanks attached to the 
inside of the hull. Figure 9a illustrates the interior transducer mounting ina 
38 foot steel launch and Figure 9b shows the installation in a 16 foot fiber- 
glas skiff. The individual transducer tanks were mounted athwart-ship near 
the keel and just forward of the center of the boat. Each tank was tailored to 
fit the dead rise of the hull so that the transducer reflector cone rested in the 
horizontal plane. After installation of the transducers, the tanks were filled 
with castor oil which doubles as the sounding medium and anti-freeze. 

Utilization of this type of installation requires only the transfer of the re- 
cording unit and battern when changing the survey work from one craft to an- 
other. It also makes possible reconnaissance soundings at speeds of 15 to 20 
miles per hour. During a recent reconnaissance survey, the sounding of the 
thalweg profile of a 100 mile reach of submerged river channel, plus the 
sounding of 20 reservoir ranges, which totalled about 20 miles in length, 
(positioning along the ranges was determined by current meter) was accom- 
plished by two men in 12 hours. 

The accurate location of the boat on the range poses a problem which has 
been ingeniously solved on the TVA surveys by the use of stretched piano wire 
as described by the author. However, the question is raised: can the weight 
and bulk of this equipment be justified to obtain a distance measurement ac- 
curate to 2 feet in 6000 feet when a simultaneous depth measurement of 50 
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feet could have been reflected from any point within a 27 foot diameter 
circle? It is believed that the 6000 foot distance as measured by means of 
the piano wire approaches the maximum feasible distance for this method; 
certainly it would be unsuitable for Missouri River reservoirs where the 
ranges frequently exceed three miles in length, and for use in the open river. 

In the Missouri River Division the distance along the range line is fixed 
by one of three methods: 


a) By direct stadia measurement on short ranges (Figure 10). This method 


is extensively used on the open river where the sight distance is less than 
2000 feet. 


b) By current meter measurement on long ranges in still water. Insofar 
as the writers know, this method of distance measurement was initiated by 
the Southwest Division of the Corps of Engineers. The current meter is 
suspended from a boom near the bow of the boat (Figure 11) and calibrated 
under field conditions. (The distance computed from the standard current 
meter formula will vary from the field distance by an amount depending 

on the field installation). Using the field calibration, check measurements 
have indicated a maximum error of 1 per cent for this method. 


c) By Raydist measurement. This is an electronic device which after cali- 
bration measures the distance between a master station on the boat and a 
relay station on the shore (Figure 12). Over open water it is accurate to 
within 0.02 percent (1 in 5000); however, its accuracy decreases sharply 
when used within about 600 feet of the shore. It has proved to be very ef- 


fective over long ranges but as its initial cost is high the Raydist unit 
could not be justified except for very large reservoirs. 


With both the current meter and Raydist measurements, the distance at the 
ends of the sounding ranges near the shores must be measured by some other 
means such as stadia or a range finder. The sounding boat is kept on the 
range by sighting range markers, and where greater accuracy is required, 
by transit and radio. 

The use of portable radio equipment, as described by the author, has been 
found to be highly beneficial in surveying, both on land and water. They pro- 
vide a real economy of operation as well as convenience. 
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HY 5 October, 1957 


Figure 7a. River survey boat Santee. A shallow draft 
boat especially designed for the study of sediment 
transport in the open river. Note the large amount of 
special measurement equipment on foredeck. 


onl 


Figure 7b. River survey boat Santee in position for 
sediment transport measurement on a 1200 foot tag line 
stretched above the river. Tag line is not used when 
sounding the river bed. 
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Figure 8. River survey boat Terrapin. With the scow 
bow it is necessary to use the outrigger sonic sounder 


antenna. Wote the outrigger mounting at stern of boat; 
the "shoe" is in the water. 


Figure 9a. Installation of transducer wells inside the 
steel hull of the reservoir survey launch Dakota. 
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Figure 9b. Installation of transducer wells inside the 
fiberglas hull of a 16 foot skiff. The wells are 
located out of the way beneath a seat. 


Figure 9c. Transducer wells in hull of fiberglas skiff. 
One transducer has been removed from the well. 
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Figure 10a. Fiberglas skiff sounding a river range. 
Distance measurement to "fixes" on range measured by 
stadia. Fixes marked and reported by portable radio. 


Figure 10b. Fiberglas skiff sounding a river range. 
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Figure 10c. Sonic sounding equipment in fiberglas skiff. 


Figure 10d. Launching fiberglas skiff. Note convenience | 
of hitch mounted on front bumper of power wagon. 
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Figure 11. Current meter used to measure distance along 


reservoir ranges as mounted on bow of reservoir survey 
launch Dakota. 


7 


Figure 12a. Raydist distance measuring equipment 
mounted on Santee. Master control station is 
inside the cabin and the transmitter is mounted on 
the stern. Note large number of antennas. 
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Figure 12b. Raydist phase meter used to indicate 
distance from boat to relay stations on shore. 


Figure l2c. Raydist relay stations on shore. Transit 
is used to keep boat on range and provide "fixes" when 
the boat is less than 600 feet from shore. 
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